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I. INTRODUCTION 


A. SATELLITE DESCRIPTION 
1. General 

The High Temperature Super Conductor Infrared Imaging Satellite 
(HTSCIRIS) is designed to perform the space based infrared imaging 
surveillance mission and is depicted in Figure 1.1. 



Figure 1.1 HTSCIRIS Satellite 
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The design of the satellite follows the ’black box' approach. The payload is 
a stand alone unit, with the spacecraft bus designed to meet the requirements of 
the payload as listed in the statement of work presented in appendix A. 

2. Payload 

The payload weighs 362.9 kg and has the dimensions of 1.52 x 1.52 x 0.51 
meters. The payload attaches to the spacecraft bus on one 1.52 x 1.52 meter 
square surface. On the opposite face is mounted the infrared telescope and 
associated radiators. The infrared focal telescope has a 4 x 4- field of view, and is 
mounted at a 45” angle to the surface of the payload. Ibe telescope has a hmged 
protective cover that opens approximately 50“ after the satellite is deployed. 

The entire cryogenic cooling system is contained within the payload, and is 
designed to maintain the infrared detector at 65 degrees Kelvin during sensing 
operations. When the detector is not sensing, the cryogenic refrigeration system is 
designed to replenish the cold reservoir of liquid nitrogen. During the sensing 
mode, the mechanical portion of the cryogenic cooling system is not operated and 
the cold reservoir is used to maintain the infrared detector at the required 65 

degrees Kelvin. 

The orientation of the payload continuously measured by an inertial 
measurement unit. This IMU is periodically updated through two star sensors 
which are mounted on the front (positive x -axis) face of the payload. Finally, a 
data processor is included within the payload to format the infrared data for 

transmission. 

a. Payload Modes 

The satellite has two modes of operation which include the sensing 
and standby modes. The 3 axis attitude control system is designed to track the sun 
in the standby mode and nadir point during the imaging periods. 


2 


(1) Sensing Mode 

The satellite payload will be earth pointing during sensing. The 
spacecraft body frame will track the moving local vertical reference system, with a 
specified roll axis orientation for the infrared device to view the destred region. 
The system is designed for 15 minutes of sensing time for two consecutive orbtts 
each 24 hour time period. Ibis sensing time limitation is a function of the capac.ty 
of the cryogenic refrigeration system. 150 watts of power is required by the 
payload during this mode. Real time image data transmission to the mobile ground 
site also occurs during sensing. There is no store and dump capability on board the 

spacecraft for the infrared information. 


(2) Standby Mode 

In the standby mode, the spacecraft will track the sun, with the 
solar arrays to charge the batteries. Simultaneously, the payload and telescope 
are pointed away from the sun to optimise thermal control and cooling 
requirements for the payload. 100 watts of power is required by the payload during 

this mode. 

3. Spacecraft Bus 
a. General 

The general configuration of the spacecraft is a 1.524 meter cube with 
the solar panels at the base and the payload located on the top as previously 
depicted in Figure 1.1. The four solar panels are folded for launch and deployed 
following separation from the launch vehicle. The high data rate antenna e 
above the payload and the two pairs of omnidirectional turnstile antennas are 
placed on opposite sides of the the spacecraft. Total weigh, of the spacecraft is 
approximately 1000 kg. The launch vehicle adapter attaches to the central support 

cylinder and thrust cone at the base of the spacecraft. 
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b. Attitude Control 

The attitude control system is designed to continuously stabil.ee t e 
spacecraft while operating in the sensing and the standby modes. The attitude 
control system consists of four reaction wheels for three axis control. Three o e 
wheels are located on orthogonal axes, with a fourth wheel skewed a. 45 degrees 
ou, of plane for redundancy. Momentum dumping is achieved via orthogona 
magnetic torque rods. Eight hydrazine thrusters are used for despm a ter 
deployment from the launch vehicle and to provide orb.tal correct, ons o 
compensate for atmospheric drag. The thrusters are also available for detumble 
mode and as a backup to the torque rods. The attitude control computer rece.ves 
data from the IMU. earth sensor, sun sensors, star backers, and the long term 

gyro assembly. 

c. Propulsion 

The propulsion system is composed of two sets of four 1 Newton 
hydrazine thrusters. The system is designed for reboost and detumble operates 
as we,, as backup to the torque rods for momentum wheel desaturation. Suffice" 
fuel is left over at end of life for deboost. The system is designed to "tamtam u 

capability in the event of single thruster failure. 

d. Electric Power 

The electric power system will provide power to the payload and the 
bus for al, anticipated operating conditions. While operating in the sensing mode, 
the entire satellite is designed to operate completely on battenes. In the Stan y 
mode, the solar arrays provide the power required for housekeepmg an to 
recharge the batteries prior to the next sensing period. Four 10 Ah n.ckle- 
hydrogen batteries provide power a. the bus voltage of 28 ± 4 volts via a full shun, 
regulator. The solar arrays are mounted on four panels that deploy fol.owmg 
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" ZZL* — - * Satellite Control Network for 

i Ttun cets of turnstile antenna pans 
satellite tracking, telemetry and ^ ^ ^ ^ ^ ^ ^ site via 

provide omnidirectional coverag % 

L wtde kand high data rate — r. ^1:^. 

was selected and utilized the binary phase shtft keytng (B ) . _ ^ ^ 

The system was designed with a 3dB safety margtn. ides 

parabolic antenna design. One 0.5 meter antenna on a tnpod s. 

, The 2 1 degree half power band width provtdes a 19-63 km 

full coverage, lne uc & 

footprint. 

f Thermal 

3 v a hnth nassive and active techniques to 

The thermal design utilized both p 

t rp limits Strip heaters with temperature sensors 
maintain the bus with temperature limits. P 

I nsed to maintain the minimum tem.rature during the cold — ~ 
batteries require louvers to maintain the temperature range, 
payload is thermally isolated from the spacecraft bus. 

g ' tI W spacecraft structure was designed to withstand the launch loads 

. . , T . nre dicted launch vehicle loads are 7.2 g axial 
rrf the Delta II launch vehicle. The predicted tauiu, 

of the Delta used for all critic al load carrying 

and 2.5 g laterally. A factor of safe y o . ^ ^ ^ 15 

structures The fundamental launch vehicle natural frequency 1 

n, ™ ^ t : 

. ... TV .... „ ,V ... -f 

conical adapter. The sides 01 onacecraft is 

panels with equipment attached directly to those same pane s. 
designed for a 50 Hz minimum frequency. 
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B. MISSION REQUIREMENTS 

1. General , h SDaC ecraft bus were originated 

Specifications influencing the design 

a, the Naval Research Lab and are summarized below in Table 1.1. 

Table 1.1 General Design Requirements 


Orbit Inclination L 

m A11 bus subsys 7 

■ m „„, The satellite orbit will be circular with an 

were designed to meet this require ■ Launch ve hicle selection will 

altitude of 509.3 km and an incl.natio « ' weigh, 

be discussed in greater detail later in the report and ,s based 

and dimensions, as well as the specified orbit. 

1 Attitude Control Subsystem 

l. muiuut ^ i _ xhp first is to slew 

The attitude control subsystem has three principal tasks. The 

====== 

the sensing mode to ensure proper image r f 

fhra rlpdred geographic area. The Attuuuc v. 
the telescope to view the desired geog y 

requirements are specified in Table 1.2 
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Table 1.2 Attitude Control Specifications 


Method Of Control 

3 Axis Stabilized 

Pointing Accuracy 

± 0.5 Degrees Each Axis (3o) 

Slew Time 

90 Degrees In 15 Minutes 

Settling Time 

1 Minute 

Rate stability 

0.003 Degree Per Second 


Information is also available to the spacecraft bus from the dual star sensors as 
well as the Inertial Measurement Unit that are contained within the payload. 


3. Radio Frequency Communication Subsystem (RFCS) 

a. Omnidirectional Transmit and Receive Capability 

The communications subsystem is designed to provide dual redundant 
communications with the ground control station. The low data rate transmission is 
designed for omnidirectional command and control of the satellite and is capable of 
transferring data at a minimum rate of 16 KBPS. Information to be transmitted 
includes spacecraft operation and housekeeping data as well as pointing 

commands. 

b. Directional High Data Rate Transmitting Capability 

The high data rate downlink is designed to transmit image data in real 
time to a remote mobile ground station. Requirements for the Radio Frequency 
Communications Subsystem (RFCS) are specified in Table 1.3. 
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Table 1.3 Radio Frequency 


Communications Subsystem Requirements 



4. Command and Telemetry Subsystem (CATS) 
a, Command Receive 

The information received by the RFCS from the ground site is used by 
the CATS subsystem to execute commands to the spacecraft subsystems. 
Specifically, the CATS subsystem will decrypt the ground commands, then 
validate, decode, and store them. The information is then distributed to the 
respective subsytem for implementation. 
b. Transmission 

The CATS system will perform the reverse operations as above for 
data transmission via the RFCS. The data to be processed, formatted, and 
encrypted includes the spacecraft operation and housekeeping information, and the 

image data. 

5. Electric Power Subsystem (EPS) 

The electric power subsystem must be sized to support 15 minutes of 
sensing time per orbit for two consecutive orbits during the worst case orbital 
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eclipse period. In the sensing mode 150 watts is required for payload. In the 
standby mode, 100 watts is required for payload operation. 


C LAUNCH VEHICLE DESCRIPTION 
' The Delta II launch vehicle with medium shroud was selected based on 

payload mass and dimensions. The ideal launch site is Vandenberg AFB, 
California, litis site minimizes flight over populated areas during launch and ts more 


suitable for high inclination launches. 
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II. SPACECRAFT 


CONFIGURATION 


A. mission requirements and the 

T- spacecraft is comprised of two " dictat ed t „ 

bus. The specifications and performance cn Research Laboratory 

the design team by the representattves could 

statement of worlc and the goa, of the destgn was to 
support die payload from launch through end-of-ltfe. 

B. description 

1 General meter tall cubic 

The payload is described as a 1 .524 meter square by 0.508 

ot The imaging system and radiators for the 

structure and is depicted in Figure 2.1. The .mag y The 

sensor thermal control system are located on one o thesqm. 
payload-to-bus mount, ng interface is located on die other square ■ 

sensors are mounted on one of the side panels. 

2. Satellite Bus 

TheTufis also a bos-lilce structure that is 1.524 m square by 1.016 m 

( „ The bus is required to provide a stable platform to orient the sensor or 
The 0 aho required t0 p r0 vide electric power, 

imaging and standby mo . transmission of the 

. telemetry downlink as well as reai 

command up, in ^ ^ . cylind rica, aluminum monocoque 

imaging data. The bus is consuu The 

thrust tube tha, is attached to an aluminum monocoque cornea, ad p^ ^ ^ 

cone protrudes through the face opposite ca rry the 

provides the interface to the launch vehicle. The thrust 
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spacecraft launch loads. Four aluminum honeycomb side panels surround and are 
attached to the thrust tube and cone. They act as mounting surfaces for equipment 
and provide radiating area for the thermal control of the bus. An additional panel is 
located along the payload mounting face and provides a mounting surface for the 

payload. 



Figure 2.1 Overall Configuration 
b. Solar Arrays 

Fixed solar arrays are located parallel to the bottom face and are 
mounted on panels that extend outboard from the bus. The arrays are 0.38 m wide 
by 1.47 m long. The silicon solar cells are mounted parallel to the anti-earth face 
and are designed to provide maximum power to recharge the nickek hydrogen 
batteries as well as supply house keeping power while operating in the standby 
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mode. During the imaging mode the entire spacecraft power is provided by 

rechargeable nickel hydrogen batteries. 

c. TT & C and Comunications 

Telemetry and communications are provide by the bus through one 
Wide Band Downlink antenna (WBDL) , two omni-directional transmit antennas 
and two omni directional receive antennas. The omni-directional transmit and 
receive antennas are paired up on the pitch and anti-pitch faces to prov.de 
continuous telemetry and command uplink communications capability. The WBDL 
antenna is gimbal mounted to a boom that extends along one corner of the 
spacecraft, past the payload face, to allow horizon to horizon coverage during 
imaging evolutions. The support equipment, transmit/receivers, and amplifiers are 
mounted interior to the bus. The high heat dissipating components are mounted to 
external panels. Two interior aluminum honeycomb equipment panels (oriented 
parallel to the square faces of the bus) provide mounting surfaces for the non- 

critical equipment. 


d. Attitude Control 

Attitude control for the spacecraft is provided by four reaction wheels 
and eight magnetic torque rods (for momentum desaturation). The reaction wheels 
are mounted near the top of the bus, three inside the adapter cone and the fourth 
(the yaw wheel) inside the bus on the top of the uppermost equipment panel. 
Three of the wheels are oriented on mutually orthogonal axes with the fourth at a 
45 degree angle from the other three. The torque rods are mounted, two per axis, 
on opposite side panels of the bus toward the top surface. Initial orbit correction 
and emergency attitude control (despin) are provided by eight 1-lbf hydrazine 
thrusters, mounted in groups of four on the pitch and anti-pitch faces of the bus. 
The thrusters are utilized as a back-up to the primary attitude control sy 
addition to station keeping. The fuel and pressurant tank are located at the center 
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of the adapter cone with the surface of the tank aligned with the top face of the 
bus. 

e. Thermal Control 

The thermal control system is an active/passive design that maintains 
the bus and components at nominally low temperature, and provides active heating 
to maintain temperature within limits. The system is composed of passive rad.aung 
and insulating material, applied to the surfaces of equipment and structures, and 
electric heaters attached to temperature sensitive components. Louvers are 
attached to the batteries as an aditional measure to dtsapate more heat energy. 

C. MASS SUMMARY 

The mass properties for the principle subsystems are listed in Table 2.1. 

Table 2.1 Satellite Mass Summary 


SUBSYSTEM - 

MASS 

PAYLOAD 

362 kg 

T T & C 

148.5 kg 

FT FfTRlC POWER 

97.0 kg 

PR OPT IT .SION (DRY) 

8.3 kg 

ATTTTTJDE control 

42.4 kg 

[thermal control 

24.8 kg 

structure 

164.9 kg 

» rr'/^TT A VfTr 1 AT TNTFGR ATION 

73.3 kg 

MECHANICAL 1tN 1 CrTJivrk i 

SUB TOTAL 

921.2 kg 

MASS MARGIN (20%) 

184.2 kg 

PROPELLANT 

12.9 kg 

TOTAL 

1118.3 kg 
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D. MOMENT OF INERTIA CALCULATIONS 

Moment of inertia and center of mass calculations for the spacecraft were 
conducted using the Microsoft Excel (ver. 2.2) spreadsheet software. The 
complete spread sheet is included in Appendix B. The reference coordinate 
system used to define the positions in that spread sheet is shown in Figure 2.2. 



Figure 2.2 Moment Of Inertia Reference Frame 
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Calculations of the components were made by idealizing the components rnto 
one of several regular geometric shapes. Component moments of inertia were 
then calculated about die principle axes oriented parallel to the principle axes of the 
spacecraft. Parallel calculations were also made to determine the spacecraft 
center of mass. Smaller components of the spacecraft were modeled as point 
masses and only their displacement from the spacecraft axes were considered^ 
The final calculation of the spacecraft moments of inertia were calcu.ated w„h 
reference to axes that pass through the spacecraft center of mass and are 
perpendicuiar to the faces of the spacecraft. The location of the center of mass ,s 

listed in Table 2.2: 

Table 2.2 Center of Mass Location 


Note that X, Y, and Z pass through the geometric center of the 1.524 m cube and 
are oriented positive along the roll, pitch, and yaw axes. The moments of inertia 

for the system are summarized in Table 2.3: 

Table 2.3 Moment of Inertia Summary 


Z = -0.0107 m 



E. ELECTRICAL POWER REQUIREMENTS 

The spacecraft operates in two distinct modes: sensing and stand-by. During 
the sensing mode the sensor draws it s f*ak power, die dam is being transmuted in 
real rime and the bus is still drawing general maintenance power. Dunng this mode 
the batteries supply the entire load to the spacecraft. In the stand-by mode the 
sensor and payload draw less power, however the solar arrays are required to 
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recharge the batteries for the next sensing period. The peak power requirements 
for the various subsystems are listed in the Table 2.4: 

Table 2.4 Electrical Power Summary 


SUBSYSTEM 

STAND-BY 
LOAD (Watts) 
100 

SENSING LOAD 
(Watts) 

PAYLOAD 

150 

ATTITUDE CONTROL 

50 

50 

THERMAL CONTROL | 

70 

30 

PROPULSION 

1 

1 

TT & C 

25 

25 

DOWNLINK 

0 

140 

TOTAL 

246 

406 


E. EQUIPMENT LAYOUT 

The following sections show the general layout of the equipment on the various 
structural panels. A more complete description of the individual components is 
provided in the following chapters. 

1. Roll Face 

Primary components of the roll face are depicted in Figure 2.3 
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Figure 2.3 Primary Roll Face Components 
2. Anti-Roll Face 

Primary components mounted on the Anti-Roll face of the spacecra 


depicted in figure 2.4. 
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Figure 2.4 Primary Anti-Roll Face Components 
3. Pitch Face 

Primary components mounted on the Pitch face of the spacecraft 
depicted in figure 2.5. 


are 
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Figure 2.5 Primary Pitch Face Components 
4. Anti-Pitch Face 

Primary components mounted on the Anti-Pitch face of the spacecraft are 
depicted in figure 2.6. 
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Figure 2.6 Primary Anti-Pitch Face Components 


5. Anti-Earth Face 

Primary components mounted on the Anti-Earth face of the spacecraft are 
depicted in figure 2.7. 
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HI. ORBITAL DYNAMICS 


The infrared imaging mission of this satellite dictates the orbital parameters 
for this spacecraft. These parameters were defined within the statement of 
work by NRL and included in Appendix A. The orbit was given a 509.3 km 
altitude and a 70“ inclination as summarized in Table 3.1. The ground path of 
this high inclination orbit is depicted in Figure 3.1 and shows that the satellite 
maximizes its time over land masses and passes over the majority of the 
habitable land on the Earth. These characteristics of this orbit make the choice 
ideal to successfully carry out the mission of imaging a wide variety of targets 

with both land and water backgrounds. 

Table 3 . 1 . Summary of Orbital Parameters 


Parameter 

Value 

Parameter 

Value 

Apogee 

6887.43 km 

Perigee 

6887.43 km 

Period 

1.58 hrs 

Inclination 
Longitude of the 
Ascending Node 

70° 

Argument of Perigee 

N/A 

TBD(l) 

Eccentricity 


Altitude 

509.3 km 


Note: (1) Longitude of Ascending Node determined by launch date. 


A. ORBITAL PERTURBATIONS 

The statement of work defined the allowable tolerances and are defined as one nautical 
mile of orbital altitude (1.852 km) and 1° of inclination. For this particular satellite the drag 
force on the spacecraft is the primary factor in maintaining this tolerance over the three year 
lifetime. The inclination is effected by several different forces, these include the Earth s 
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UEHIC1E CROUMBTBftCK 



Figure 3.1 Orbit Ground Track 
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LOWJITBDE Cdejr> 


oblateness, and the attraction of the sun and tnoo. Since the orbit is eclat land there 
oo perigee or rotation of the argnnten, of perigee within the orhitai plane however, 
orbit's ground trade is affect by the rotation of the longitude of the amending node. 

1 . Drag 

The accuracy of the drag force on the satellite and its effect on the or it are as 

accurate as die prediction of die atmospheric density. Two methods, hand C ^ U ^ 01 ^ 
the ASAP program, are used to predict the dug perturbations on the orbit. The 
program estimated the larges, altitude loss and will be used for orbital matntenance 

below in Table 3.2 

Table 3.2 Summary of Altitude Loss for Spacecraft 


Method 

Estimated Altitude Loss 

Hand Calculations 

1.65 km/yr 

ASAP Program 

4.50 km/yr 


2 . Inclination 

The effect of the Earth, Sun, and Moon on the inclination of the orbit is complex 
in nature. These effects not only change as the satellite orbits around the Earth but also 
change as the Moon routes around the Earth and as the Earth routes around the Sun. e 
rotal effect of inclination change on the orbit is based on the maximum presented values 
for effect from the Sun and Moon, it will be shown tha, for the circular orbt, the 
perturbations from the Earth's oblateness has no effect on the inchnanon of the orb,,. 

a Effects from Earth's Oblateness 
’ As mentioned above there is no cumulative effect on the inclination of the 

orbit. The Earth's equatorial bulge does exert an influence on the satellite, but this is small 
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in nature and has a period equal to the orbital period of the satellite. A detailed derivation 

and analysis are include in Appendix C. 
b. Effects from Sun 

The Sun's influence on the inclination of the orbit depends on four 
parameters, the apparent inclination of the sun, the inclination of the orbit, the altitude of 
the orbit and the longitude of the ascending node. Because the altitude and inchnatton of 
,he orbit are specified, the Sun's effect on the inclination will be influenced only by the 
apparent inclination of the sun and tire position of the satellites longitude of the ascendtng 
node. For design purposes interest lies in the maximum value of di/dt. Equation 3.1 shows 
that the maximum value for di/dt will occur when the Sun's inclination is ±23.5" and ts a 

minimum when the Sun’s inclination is 0 . 


luiium ▼-***»*■— 

dL = lb. if sin(Q) cos(Q) sin(i) sin(i s ) 2 + sin(Q) cos(i) sin(i s ) cos(i s )J 
dt 2 h j3 l 

s 


(3.1) 


Figure 3.2 shows the inclination change as a function of Q when the Sun's apparent 
inclination is set to 23.5", from th e graph the maximu m value is found to be .0084 deg/yn 



Figure 3.2. Perturbation Effect Due to Sun's Influence 
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c. Effects from Moon 

The influence of the Moon on the orbital perturbation is similar to the effect 
of the Sun. The magnitude of the change in inclination per year is slightly hig 
Moon's influence than i, is for die Sun's. The governing etpiadon for di/d, is shown to be 

di.lh.dr sin(£My cosfO-O,) sin(i) sinfi,) 2 + sin/O-O,) cos(i) sin(i,) cosfi,) 
dt 4 h rj 


+ sin(26,) [ -cos(2(£l*Qj)) sin(i) cosfi,) + cos(fi-£2,) cos(i) sinfi,)] 
+ cos(20,)tsin(O-O,) cosCQ-n,) sin i + sinffl -ft,) cosffi -G,) 
x sin(i) cos(ij) 2 - sin(Q-£2,) cos(i) sm(i,) cos^)]] 


(3.2) 


where all of die T subscripted variables apply to the Moon's orbit. Equation 3.2 contains 
bod, secular and period rerms, die terms dependent on 0, am periodic with a penod of 14 
days, half of the lunar cycle. Equation 3.3 shows only the secular terms and is dependent 
on both the inclination and the the Moon's longitude of die ascending node. 

dL = lhd[ sin(Q-a) cosfQ-Q,) sin(i) sinfi,) 2 + sinffl-Q,) cos(i) sin(i,) cosfi,)] (3.3) 
dt 4 h 

The lunar orbital inclination takes 9 years ,0 complete a cycle of inclination change from a 
maximum value of 28.6° to a minimum value of 18.3°. during this nine years the longitude 
of the ascending node is varying ±13°. Fortunately, a. die maximum and minimum 
inclinations die longitude of die ascending node is zero. Figures 3.3 and 3.4 show the 
inclination change on die satellite as a function of the satellites longitude of die ascending 
node. The maximum value occurs when the Moon's inclination is a, 28.6° and contnbutes 
tt> the change in inclination .023 degrees per year. Because this is the maximum value 
encountered all cubital maintenance calculations wUl be based on die .023 degrees per year 

value. 
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didt(£l) 


0 . 03 


Maximum .023 degrees per year 



Figure 3.3 Moon's Influence on Di/dt (ii-28.6°) 
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3 Longitude of the Ascending Node 

‘ we arc interested in the longitude of the ascending node for the satellite for two 

reasons; firs,, all of the calculation for the inclination change arc based on dte spacecraft s 
longitude of dte ascending node and second, it can be ca.cuiated how much the ground 
track shifts on daily basis. The equation governing the change in the longttude of the 

ascending node can be approximated by 


dU = -2.062xl0 14 
dt 


cos( i ) 


3.5 


(3.4) 


For this orbit the longitude of the ascending node shifts -2.6 degrees per sidereal day. This 
orbit completes 15.14 orbits per sidereal day, computing the shift of the ascending node 
15 orbits (-2.57 degrees) and then computing the remaining degrees of rotation the Earth 
has at the 15 orbit point (+3.419 degrees) it can be shown tfia, the orbit ground mack shifts 
84 degrees east every 15* orbit. This orbit nearly passes over the same ground hack 
every day and a, tire equator the subsateUite point shite by 50 nm every 15* orbit. 

B. SWATH WIDTH 

Because tire satellite can only image within its visible horizon and nansmission of dam 
from tire satellite to a ground station must occur when they are visible to each other, tire 
swath width of tiie spacecraft is of interest. Defining three variables, p - angular radius of 

Earth, X - Earth central angle, and n - "ati* * e swath width Ca " * ca)culattd 
3.5). Equation 3.5 compares tire angular radius of the Earth based on the orbital radius (r) 

of the satellite, where R equals 6378.14 km. 

sin (p) = 7 3.5 

A, a 509.3 km altitude p (angular radius of Earth) is computed to be 67.8°, *is is the angle 
off nadir tire satellite sees the horizon. Because the grazing angle requirements for a ground 
station ,o communicate with tire satellite is different than tire graring angle of tire re observe 
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an object on the ground there will be two different swath widths to define. To compute the 
swath width the spacecraft will have for a ground station to observe the satellite at grazing 
angles greater than 10° a nadir angle of 65.7° will define the edge of the swath. While the 
statement of work restricts the grazing angle for observation of an object between 45° and 
90°. Using Equation 3.6 the remaining Earth central angle can be found, this angle 
represents half the angle the swath subtends on the Earth's surface. 

sin(p) sin(X) 3.6 

laJl<T ^ ) = 1 - sin(p) cos(X) 

The Earth’s central angle for the ground station to see the satellite is 14.2 and for 
observations it is 4.1°. This defines swath widths of 3165 km and 91 1 km for the ground 
station and observations respectively. Realizing that die subsatellite poin. is moving across 
the Earth's surface at 7.6 km/sec these distances become very important 



C RESULTS 

Drag on the satellite is a factor which must be accounted for in the fuel consumption 
for orbital maintenance, the predicted altitude loss for the 180 day period from the ASAP 
program was nearly equal to tile allowed altitude perturbation and tire sprreecraft must be 
reboosted to maintain the orbit. The maximum inclination change seen by the satellite did 
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not exceed the allowed tolerance, the total for both the Sun and Moon influence is .0314 
degrees per year and over the lifetime of the satellite is approximately 0.1 degree. No 
additional fuel is necessary for the shifting longitude of the ascending node, because it does 
not concern us that the ground track shifts and the calculations were done to inform on the 

nature of this orbit. 
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IV. RADIO 


FREQUENCY COMM SUBSYSTEM (RFCS) 


A. MISSION REQUIREMENTS „ the satelUl e and 

The RFCS communication subsystem is Teleme try Subsystem 

the ground station. The RFCS includes the Comman b and the 

«*™ «■«. 

gr „„„d station. As described^ the^ ^ ^ ^ (ollowing capabilities: 

Communicabon Subsystem ( 

. Omnidirectional command receiv.n g capability. 

. omnidirectional low da. rate (16 Kbps minimum) transmitting 

capability: 

. Directional high da. rate (150 Mbps maximum) transmitting 

capability. /vmRDSl is required to provide the 

The High Data Rate Down Link System (HDRDS) « requ 

following capabilities: 

. Optimal spectrum usage for 150 ps. 

. Bit Error Rate (BER)lOe-ll encrypted 

Bit Error Rate (BER) 10*9 unencrypted: 

. 99% availability due to rain attenuation: 

i» rtf 10 decrees at data collection 
. Minimum elevation angle of 10 degrees 

terminal: 

20 foot (6 meter) diameter receiver antenna is assum 
the ground station: 
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B. SYSTEM DESCRIPTION 
1. System Integration 

The communications subsystem is the interface between the satellite and 
the earth. The system's components are located as depicted in Figure 4.1. It 
generates a downlink RF signai that is phase coherent, transmits the downltnk error 



Figure 4.1 Communications Subsystem Components 
so its phase synchronizes with the received phase of the uplink carrier. This is 
referred to as the coherent turnaround or two-way-coherent mode. The turnaround 
ratio using SGLS is 256/205. The two-way-coherent mode allows the ground 
station to know more precisely the downlink signal's frequency and to measure the 
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Doppler range rare. This allows a sweep of fewer frequencies and thus acquire 
the spacecraft more quickly. The spacecraft will have large volumes of data and a 
short field-of-view time to the ground station. To transmit maximum data on a 
direct downlink, the ground station must acquire the signal in a minimum amount of 
time. By using a ranging method of navigation, range-rate information ,s obtained 

from the Doppler shift of the coherent signal . 

The communication subsystem s transponder is compatible with SGLS. 
tracks the uplink carrier, receives and detects commands, and transmits telemetry. 
The communication subsystem contains two transponders with parallel transmn 

and receive signal paths for redundancy. 

Command and Data Handling (C&DH) on board receives, decodes, 

processes and distribules spacecraft commands. It also gathers.formats, stores 

and transmits telemetry data from spacecraft measurement. As depicted in Figure 


4.2, the C&DH 
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AA 

DIU 

DT 

HD 

R 

R/D 

RIU 

T 

TCU 

WBT 


: Antenna Actuator 
: Data Interface Unit 

*: Downlink Transmitter 

: Hybrid Divider 
: Receiver 

: Receiver /Domod 

: Remote Interface Unit 
: Transmitter 

: Telemetry Command Unit 
: Wide Band Transmitter 


Figure 4.2 C&DH Subsystem 

subsystem consists of a Remote Interface Uni, (RIU), Tracking & Command Units 
(TCU), and Data Interface Units (DIU). Encryption and decryption devices are 
also included for which NSA distributes and regulates the keys. The TCU 
receives the demodulated uplink information and routes it to the DIU which also 
receives, formats and routes downlink telemetry, through encryptors, to the 
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transmitter. The TCU serves as decoder; authenticator for commands, a 

«, Remote interface units receive and process commands and revests * 

date. The TCU has a built in processor that performs calculations, 

and issues commands. Pulse Code Modulation (PCM) subcarrier provtded by D 

is sent to the downlink transmitter. (DT) Hybrid DividersfHP) 

or from the appropriate set of turnstile antennas The primary rece 

could be used as a back up transmit antenna, and the primary teansmu - 

could be used as a back up receive antenna. Dual redundancy ts provtded 

design. 

2. High Data Rate Down Link 

a. Description 

The Wide Band High Data Rate (WBHD) Transmitter ts destgned 
provide the capability of sending real time digital imagery data a. a rate of 150 

L„,„ «<*.>■ ~ **■ “-*• • > * 

a. - '""l 

This margin becomes important when you consider the uncertainties mvove ^ 
the receiver, a mobile platform (assumed to be a Van). There are tncrea 
problems with such platforms due to incalculable environmental interference » 
equipment calibration. This is especially true if you're us.ng a mob. 
diameter antenna (assumed to be a fold out umbrella design). The requ. 
energy per noise power spectra, density was calculated based on Quadrate 
Phase Shift Keying (QPSK) modulation. This consists of 
degree phase shifts for binary data transmission. Binary ase 
(BPSK) and Multiple Frequency Shift Keying (MFSK) were also considered 

“ a a Q«k - » « - " di '" 1 
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transmission. Using 20 GHz for a carrier frequency allows up to 1 MHz band 

width for data transmission. 

b. Antenna Design 

For the High Data Rate Wide Band Down-link antenna design, two 
parabolic antenna configurations were considered. Parabolic reflector antennas 
offer narrow beams over a wide range of frequencies. They are also simple to 
design, simple to construct and have a proven track record. For an offset parabolic 
antenna the reflected beam is not intercepted by the feed horn which reduces the 
side lobes The reduction in side lobes increases antenna efficiency. An offset 
parabolic antenna is depicted in Figure 4.3. Due to the fact that a feed horn can not 

be fixed 



Figure 4.3 Offset Parabolic Antenna 
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to the payload, the decision was made to use a 
Antenna design as depicted in Figure 4.4. 


Front Fed (symmetric) Parabolic 



Figure 4.4 Parabolic Antenna Design 
This simplified the antenna support structure design and allows a larger down-link 
FOV. The antenna will also have two three degree of freedom gimbal attachme 
to the supporting structure to provide a wider range of pointing angles. The Front 
Fed Parabolic Antenna will provide a 19 to 63 Km diameter footprint with a 2.1 
degree Half Power Beam Width (HPBW) as depicted in Figure 4.4. Antenna 
pointing control will consist of an open loop, on-board attitude and computer steer. 
This design provides +/- 0.2 degrees of pointing accuracy with about 1 db loss, 
assuming the attitude and control system maintains within spacecraft attitude 
design tolerances. Using a closed loop system would require the use of an auto- 
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licated three dimensional up-link beacon which 
.rack receiver and a more c ^ additional cost and complexity 

must be constantly tracked. It was 
of this system was not required for this mission. 

2. Tracking, Telemetry, And Contro 

o. Description telemetry are to 

equipment status and also alerts perso lhe events , and ways to 

failures or anomalies, data is used to det ® rm '" Te)emetry can also be 

counteract or alleviate the problems pr uc ^ and predict its 

used to analyze any degradation t at m of tw0 sets of 

effect on the spacecraft lifetime. Fo P Handling system (C&DH) 

omnidirectional turnstile antennas, C«-» 

and the Air Force Satellite Contro. Network (AFSC ). 

- — *— *• 

turnstile antenna consists o provides the necessary 180 

— —• « - “ TIC. “ L- ~ - - - - 

~ ~ 7 ;::;::; 

spacecraft a distance of X/4 (0.03 mi 
minimal gain as depicted in Fig 4.5. 
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coax carter corrector 
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detail of single dipole, folded BALUN 

Ref : (3) MBIgan, ’Modem Antenna Design’. 



Figure 4.5 Turnstile Antenna 
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The cone of coverage is +/- 125 degrees from center with a 4 db drop in gain at 
the 125 degree angle. The reason for using two pair of Turnstile Antenna is 
because of the different frequencies/wavelengths used for transmit and receive. 
The pattern for a three-axis-stabilized satellite is a cardroid of revolution and is an 
omnidirectional antenna except for a cone of silence where the spacecraft 

shadows the antenna 2. 

3. Ground Control and Frequency Selection 

Hie TT&C for the spacecraft was designed to be compatible with the Air 

Force Satellite Control Network (AFSCN)l. AFSCN consists of seven Remote 
Tracking Station (RTS) located throughout the world. The system of RF links with 
spacecraft is known as the S-Band Space-Ground Link Subsystem (SGLS). file 
downlinks provide telemetry and mission data rates up to 1.024 Mbps per channel. 
There are 40 downlink channels available with frequency range of 2197.5 2297.5 

MHz. There are 20 uplink channels available with frequency range 1763.721 - 
1839.795 MHz. TT&C was designed to operate at a primary SHF carrier 
frequency corresponding to channel 10 of the SGLS ground terminal: 

Command Uplink : 1799.756 MHz 

Telemetry Downlink: 2247.5 MHz 

A back-up channel could be channel 9 or 1 1 to remain close to dipole design 
antennas. The high data downlink was designed to be compatible with 20 foot 
maximum diameter antenna. The frequency was selected to 20 GHz. The Ku-band 
frequency was chosen to provide wideband transmit data capability as specified in 

mission requirements. 

D. MASS, POWER AND EQUIPMENT LIST 

The systems power and weight are summarized in the Table 4.1 and 4.2 

below. 
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Table 4.1 RFCS Power & Weight Summary 



Table 4.2 CATS Mass And Power Summary 
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V. 


electrical power subsystem (EPS) 


A. MISSION REQUIREMENTS 

ne electrical power subsystem must provide sufficient power for 30 minutes 
sensing time per day (15 minutes times two successive passes). Dunng non- 
sensing periods, the spacecraft will be sun oriented, and for active periods, 
spacecraft will slew to an earth seeking orientation, and will continue to track the 

earth until the active session is concluded. 

The spacecraft, being sun oriented for approximately 21 hours out of each day 

is an attribute that is fully exploited in this design. The solar panels are fixed, and 
during the intervals of earth orientation, the batteries play a big role in powering the 
spacecraft's active and, if necessary, it's housekeeping load. Fixing the solar 
arrays provides a higher degree of simplicity and reliability, which is probably 

desirable considering the payload involves rather high risk technology. 

The basic concept goes like this. Put fairly large batteries in the vehicle, 
capable of carrying the active load and housekeeping load for two consecutive 
orbits After this event, return the vehicle to sun orientation, and begin charging 
die batteries and be back to full charge in 21 hours. During the development of 
the design, worst case scenarios were always assumed, battery and solar cel. 

EOL, winter solstice, etc. 

B. EPS DESCRIPTION 
1. General 

The satellite electrical power subsystem consists of four solar arrays that 
nominally generate 34.2V a, 4.3A at end of life (EOL)., four lOAh NiH 2 batteries, 
and the regulation devices to fully regulate the bus voltage. 
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2. Solar Array 

There are four arrays mounted in a stationary, fixed manner. Each arr y 
capable of filling the role as either charge array, or housekeeping array. Two 
arrays are required normally for battery charging, and one array is requtred for the 
housekeeping electrical load. The manipulation of the solar array s power 
distribution is done from ground control, via the CATS link. 

3. Batteries 

The batteries are designed to carry the electrical demand during eclipse an 
the additional load that arises during sensing and data transfer. Four lOAh N1H2 


batteries are installed. 

After a maximum possible demand event, the batteries were designed to be 
at 60%depth of discharge (DOD). The batteries will have approximately 21 hours 

to fully recharge after completion of this event. 

4* Power Control Electronics 

The spacecraft bus is regulated to 28( + /-4)V. The regulation is achieved by 
use of two partial shunt regulators on each solar panel and one senes diss p 

regulator at the battery discharge terminal. 

a . Regulation of Housekeeping and Battery Charging Voltage 

The partial shunt regulators were designed to shunt the voltage spikes 

that are experienced at BOL eclipse emergence, and dissipate the power surplus to 

space in the form of thermal energy. 

One of the partial shunt regulators on the panel is to regulate the 

voltage supplied to the housekeeping bus when that array is aligned as such. The 

other partial shunt regulator is used to regulate the battery charging voltage, when 

the array is aligned for battery charging. 
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b. Regulation of Active Payload Voltage 

A series dissipative regulator is connected in series between the 

battery discharge terminal and the payload input terminal. It must regulate the 
voltage from the battery discharge voltage at full charge (33V) down to the payload 
operating voltage of 28V. This regulator also regulates the voltage to the 

housekeeping bus during eclipse. 

C. EPS DESIGN 

1. Solar Arrays 

Since each array must be capable of filling the role as either charge array, 
or housekeeping array, the voltage supplied from each array is fixed by the battery 
charging requirement (33.9V), and the current is fixed by the housekeeping load, 
since the housekeeping current draw (8.05A) is greater than the battery charging 

rate (4.0 A). 

2. Batteries 

The batteries are 22 cell NiH 2 in individually encapsulated pressure 
vessels (CPV). This allows for an open circuit failure of one cell per battery. 
Minimum battery discharge voltage in the event of a cell failure is 24.1V, a 
maximum battery output voltage at full charge, no cell failure is 33V. See 

Appendix F for detailed development of battery design. 

The storage device is divided into four 10 Ah batteries for purposes of 

mass distribution, battery margin, and redundancy. 
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Figure 5.1 Simplified EPS Dmgram 
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D. EPS PERFORMANCE 

1. General 

The worst possible case for EPS demand at end of life yields a 60% DOD 
and total battery recovery in a 24 hour cycle. 

2. Solar Array 

a. Radiation Degradation Of Solar Cells 

Table 5.1 depicts the anticipated outputs of each array at EOL and 
BOL. A detailed evaluation of radiation degradation is provided in Appendix F. 

Table 5.1. Subarray Output 



EOL 

“ BOL ~ 1 

Pmax(W) 

148 “ ~””~1 

158 

Voc (V) 

42.8 

44.4 I 

Isc (A) 

h 469 “ ““ "1 

4.74 

Imp (A) 

4.34 “ 1 

4.36 

Vmp(V) 

34.2 

36.6 


b. Temperature Effects 

Although the temperature range is severe (-60C to 50C), the most 
significant problem faced is the extreme voltage spike experienced when emerging 
from a 35. 5 .minute eclipse. The voltage output will nearly double it's nominal value 
for one or two minutes, and rapidly return to the nominal output value. The power 
control electronics are designed to accommodate this situation. 

c . Array Sizing 

As discussed, each subarray must provide one half of the current 
demand during spacecraft passive periods (8/2)A at the battery charging voltage 
(33.9V). The number of cells in series is determined by dividing the battery 
charging voltage by the voltage at max power point at EOL. 
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0.469 

The number of strings necessary is computed in the following manner. 


d. 


Np = fL - 

(104x2.5 

Solar Cell Arrangement 


= 6.45 * 7 strings 

x6.2 

And Panel Dimensions 


The cells are arranged on a 3/8 inch sheet of aluminum. Inter-cell 
spacing is 1 mm in both directions. There are 531 cells mounted to each panel, 
arranged with the 6.2 cm side pointing in a radial fashion from the yaw axis. 


3. Batteries 


The spacecraft's design life is three years, and the batteries will be 
nominally cycled one time per day, yielding approximately 1000 cycles. Nickel 
Hydrogen batteries will be used at 60% depth of discharge (DOD). With only 1000 
cycles, this battery could be cycled to a greater DOD, therein requiring a smaller 
battery, but selecting 60% DOD provides margin, and is an attempt to deal with the 
uncertainty that exists in the ability to maintain the battery's operating temperature 
because in the operating concept as outlined, no true "north and south faces are 
maintained, and the sun will impinge upon the battery's thermal radiator at times, 
the full consequence of this on battery life has not been completely analyzed.in this 


report. 

a. Requirements 

The energy storage must be sufficient to carry the entire spacecraft 
electrical load for 220 minutes, once per day. This event also, of course includes 
two 15 minute periods of active sensing and high speed data link. 

b. Battery Capacity And Sizing 

The capacity of the storage device may now be evaluated in the 


following manner, 



r - _ _ 1 x IQ — = 38 Ah 

V db DOD 24.1 x 0.6 


The storage device is divided into four 10 Ah batteries for purposes of 
mass distribution, battery margin, and redundancy. 
c. Battery Charging 

To determine a suitable charging voltage, the correct charging voltage 
(V C ) per cell must be obtained from the manufacturer of the battery, as well as the 
voltage drop across the bypass diodes (V DD ) in the event of a cell s open circuit 

failure. „ 

V BC = V C (N-1) + N D V DD = 1-5(21) + 3(0.8) - 33.9V 

Where Nd is the number of bypass diodes in series. 

A BASIC code was written that simulated a worst case 24 hr day, winter 


solstice, EOL scenario. (Appendix F) 

4. Power Control Electronics 
a. Shunt Regulator 


Each subarray is fitted with two partial shunt regulators (PSR). One to 
bring the supplied voltage down to the housekeeping load's bus voltage (nominally 
27V) when the subarray is aligned to feed the housekeeping load. This PSR is 


called PSRh- 

The housekeeping must be regulated to 27V so that diode 1 (Figure 1) 
will remain reverse biased during sunlight periods, ensuring that the housekeeping 

does not draw it’s power from the battenes. 

The other PSR on the panel is to protect the batteries from the voltage 

spike at emergence from eclipse when the subarray is aligned to charge battenes. 

This PSR is called PSRb. Their tap points are computed as follows: 
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Vbus ) 

TP (%) = 100 (1 - Voc (ECUPSE EMERGENCE) 


Table 5.2. Partial Shunt Regulator T ap Points 

Tp 


NUMBfcK ur ^ELLS 
STRING THAT ARE SHUNTED 


PSRh 

PSRb 


37 

24 


18 


b Battery Charge And Discharge Regulator 
' A series dissipative regulator is incorporated a, the battery dtscharg 
, m ,„ at full charge is well above the specified bus 

terminal, as the battery voltage output 8 33V (full charge voltage) to 

voltage. This regulator must reduce the voltage fro 

the 28V payload voltage. 

S System Integration and Failure Recovery 

Redundancy in the electrical power subsystem offers many °P u ° n ^° 
recover from component casualties. The output terminals on ali so ™ 
t0 three position switches .ha, will feed either the —mg lo * c r 
batteries or go to open circuit. These switches are ground c— vm 
allotted word in the TT & C housekeeping link. Two subarrays mu be 
parallel for carrying the housekeeping load. One subtly ,s su ficte ™ 

— 

the remaining three subarr y . , f electrical 

tr i Tahle 9 lists the effects of various levels of electric 
again through ground control. Table 9 lists 

power subsystem failure. 
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Table 5.3. EPS Degradation Effects Matrix 


Battery 

Loss 

Array 

Loss 

Effect 

0 

1 

None 

0 

2 

Mission Loss 

1 

0 

1) Limit active sensing to 15 min. per day, or 

2) Sacrifice lifetime due to higher battery DOD 
and higher charging rate 

1 

1 

Same as above 

1 

2 

Mission Loss 

2 

0 

Mission Loss 


D. MASS AND POWER SUMMARY 

Table 5.4 illustrates the mass budget of the electric power subsystem. 

Table 5.4 EPS Mass Budget 


BATTERIES 

43 

SOLAR ARRAYS 

17 

DISCHARGE CONTROL 

12 

CHARGE CONTROL 

10 

WIRING 

9 

HOUSEKEEPING REGULATION 

6 

TOTAL EPS MASS 

97 kg 
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VI. PROPULSION 


A. MISSION REQUIREMENTS 

A launch vehicle is required to place a satellite in a circular orbit that has a 
509.3 kilometer altitude and a 70 degree inclination. TOs satellite is antic.pa.ed to 
weigh between 800 and 1200 kilograms. The satellite should be des.gned for a 
.hree year lifetime. Additionally, the requirement exists for the satellite to rema.n 
within 1 nautical mile of its given altitude and within 0.5 degrees of its tnclmation. t 
is anticipated that these requirements will make a propulsion subsystem necessary. 

B. LAUNCH VEHICLE 
1. Description 

A Delta II (7320) launch vehicle was chosen for this mission. Appendtx F 
gives the demils of the selection process. This launch vehicle is a two stage liqutd 
propelled rocket with three solid strap-on boosters. Typically a Delta II has mne 
solid strap-ons but these are not required since this satellite has such a small mas 
If desired, the six other solids could remain attached and additional satellites could 
be included in the fairing space remaining. This multiple launch would reduce the 

launch costs. 

2. Satellite Integration 
a. Delta II Fairing 

The Delta II standard shroud has a 2.184 meter diameter for a he g 
Of , 448 meters. It then expands to 2.54 meters for 2.032 meters. The payload has 
dimensions of 1.524 meters by 1.524 meters and a heigh, of 0.508 meters. The 
necessary internal shroud diameter therefore is 2.155 meters. 
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Figure 6.1 The Satellite Inside the Delta II Standard Shroud 
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Figure 6.1 shows .ha. .he payload will to inside .he shroud. The 

internal dinrensions of .he shroud includes dre payload — 

— -c and dynamic ^flection ^ — 
spacecraft lateral vibration modes are above 15 Heftz), 

thickness. 

b. Spacecraft Attachment Assembly 

nelta 11 (7320) vehicle normally uses the 60 
The two stage Delta 11 ) 

. j u c o i S24 meter diameter (to tne 

attachment fitting. I. weighs 57 kilograms and has a 1.52 

tr mnrp detail see Structural Design, 

outside edge) at the attachment point. 

Section IX. 

3. Launch Profile Canaveral or Vandenberg 

The Delta II can be launched from eit er p 

Air Force Base. Since the desired launch inclination is 70 degrees a west coa, 

launch is the best choice. The De.ta .1 for this spacecraft will be launched 

Vandenberg pad SLC-2W at a 158 degree launch azimuth. 

The Delta ,1 (7320, has two ,it,uid stages as wel, as three solid strap-on. 

The tost stage RS-27 engine and the three solid rocket boosters are 
ground a. liftoff. Following burnout of the solids, the spent cases are then jet. 

8 The RS 27 engine continues to burn until main engi 

about one second later. The RS-2/ engine 

cutoff (MECO). This takes approximately 255 seconds. 

After a short coast period, the firs, to second stage separatum bolts . 

blown, followed by second stage ignition approximately five secon s ater.^ ^ 
next major event is the payload fairing separation, which occurs ear 

second sugeft^d ^ ^ 4I0 seconds, at which time 

rr repro n occurs The vehicle then follows a Hohmann 
stage two engine cutoff (SECO 1) occurs. o Fro 1 occurs 

transfer trajectory to the desired low Earth orbit altitude. After 
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approximately 670 seconds later, the second stage is re-igntted and completes 
itsburn to circularize to the desired orbit. Satellite separation then begins 
approximately 200 seconds after stage two engine cutoff command (SECO 2). 
Figure 6.2 graphically presents this sequence of events. 



Figure 6.2 Launch Sequence 
4. Spacecraft Separation 

The launch vehicle and the spacecraft are attached by three attachment 
bolts and bolt catcher assemblies. "Upon separation, the bolts and catcher 
assemblies are retained by the spacecraft.... Following release of the three 
explosive nuts, the spacecraft/launch vehicle is stabilized by the launch vehicle 
attitude control system. Subsequently, three retaining latches are released followed 
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by retrofire of the launch vehicle yielding a minimal separation tip-off of the 
spacecraft." 1 

Fifteen seconds after the explosive bolts are fired the latches are released. 
This delay allows the angular rates to dissipate. At this point the second stage 
retro rocket fires providing the required relative separation velocity from the 
spacecraft. Expected, angular velocities at separation are a little more than 0.2 
degrees per second. This can be reduced by employing additional steps in the 
separation process. The angular velocity can be increased to 30 degrees per 
second ( within a 5% accuracy) by using control jets. 

5. Launch Uncertainties 

Table 6.1 lists the three sigma (3a) injection error parameters. Worst case 
is the one guaranteed by McDonnell Douglas. The probable worst case is what 
they actually expect for this mission. The worst case was designed for. 


Table 6.1 3o Injection Parameters 



Worst Case 

Prnhahle Worst Case 

Altitude (km) 

±18.5 

±9.5 

Inclination (degrees) 

±0.05 

±0.02 


6. Launch Vehicle Performance 

The spacecraft load factors in table 6.2 represent the combined steady 

state and dynamic values with the three sigma maximum values applicable at the 
spacecraft center of gravity. These values should be multiplied by at least 1.25 to 

obtain ultimate loads. 


l rnmmerdal Delta II Pavloa H Planners Guide, MDC H3224B, December 
1989, p. 5-3. 
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Table 6.2 Spacecraft Limit Load Factors 


Lateral Load 


Stage (g'sl 
±2.0 
±2.5 1 
±2.21-02 2 


yfain Engin e Cut Off 


6.7 to 7.2 3 

Axial Load “^±1 — 

1 Lateral load factor to provide correct bending moment at separation plan 

2. Plus indicates compression load, minus indicates tension load. 

3. 6.7 is for a 2000 kg payload, 7.2 is for a 1000 kg payload. 

C. PROPULSION SUBSYSTEM 
1 Thruster Description 

' Eight MR- 1 11C one ib, thrusters built by Rocket Research Company 
were chosen. Figure 6.3 depicts the iayou, and physical dimensions of this thruster. 
Table 6.3 gives the design and performance charactensttcs. 
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Table 6.3 Characteristics of the MR-111C Thruster 


Propellant 

Hydrazine 

Catalyst 

Shell 405 

Steady State Thrust (N) 

5.338- 1.334 

Feed Pressure (MPa) 

2.7579-0.5516 

Chamber Pressure (MPa) 

1.2066 - 0.3447 

Expansion Ratio 

74:1 

Flow Rate (g/sec) 

2.404 - 0.635 

Mass (kg) 

0.33113 

Specific Impulse (sec) 

229 - 226 

Minimum Impulse Bit (N-s) 

0 0845 (5) 2.4132 MPa & 20 ms On 



Figure 6.3 MR-111C One Ibf Thruster 

2. Thruster Locations 

The analysis of the thruster placement was based on the mission requirements 
envisioned for the thrusters. The missions are broken down into primary and secondary 
missions with an emphasis on allowing for single point failure to complete all primary 

missions. 
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Primary Mission 

Detumble of Spacecraft on orbit 

Spin down 

Delta V maneuver required for oibit maintenance 
Secondary Mission 

Backup for desaturation of reaction wheels 

Slew maneuver 

Deorbit 

A single point failure of the taster system is defined as the inability of a single thruster to 
operate in the mode for which it was intended. Two separate designs are presented each 
having their own strengths and weaknesses. The firs, design meets the design cntena of 
redundancy within the thruster system, but is more compiicated to implement and requtres 
more tasters. A second design is also introduced which is much simpier to construct, but 
offers little redundancy. Both are proposed a. this time aliowing a made-off between the actual 
need for redundancy and the requirement for simplicity. 
a. Eight Thruster Design 

The eight thruster design for the spacecraft is described as two sets of four 
thrusters, each set oriented on either the roll or anti-roll face and bracketing the center of mass 
of the spacecraft. The thmsters are canted 35.8° off vertical and oriented to thrust in the x-z 
plane (Figure 6.4). The canting of the taster provides a moment arm with components m x, 
y and r. As shown in Table 6.4 any combination of two tasters on either face creates a 
moment only about a primary axis. To produce an incremental velocity 
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Table 6.4 Thruster Operations 



Table 6.5 Thruster Placement Summary 


Thruster Number 


Location 1 


(meters) 


Moment Arm 1 Moment Created 1 


(meters) 


(1) anti-roll face 

[-.762, -.427 ,-.241] 

(2) anti-roll face 

[-.762, -.427,-241] 

(3) anti-roll face 

[-.762, .427, -.241] 

(4) anti-roll face 1 

[-.762, .427, .241] 

(5) roll face 

[.762,- .427 ,-.241] 

(6) roll face 

[.762, -.427, .241] 



(7) roll face [.762, .427,-241] [.146,.427,.203] | [1.11.-UU.54] 


(8) roll face 1 [,762,. 427,.241] | [.146,.427,-.203] [ [-1.11,1.11,1-54]^ 

Note(l): All distances arc measured from the beginning of life center of masg for the spacecraft 
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change four thruster must be fired at once. In addition to the delta v created in the roll and anti- 
roll direction, a delta v can be created in the +Z and -Z direction by firing the top or bottom set 
of thruster, but this drops the effectiveness of the thruster down to below 60%. A summary of 
the thruster location, moment arm and moment created are shown in table 6.5. As stated earlier 
these eight thrusters give the satellite the ability to withstand a single point thruster failure. The 
thrusters are placed at positions and moments calculated which correspond to the beginning of 
life center of mass, the fuel requirements for this satellite are small compared to the mass of the 
satellite and will cause minimal changes to the center of mass. Appendix F shows a sample 


calculation of the moment created by the thruster #7. 



b. Six Thruster Design 

Four thrusters oriented on the negative pitch face perpendicular to the x-y plane 
and bracketing the center of mass, along with two additional thruster on the positive yaw axis 
also bracketing the center of mass describes the six thruster design. The four thrusters on the 
negative pitch face when fired in pairs would be able to create moments about the X and Y 
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axis, all four thruster would be able to cause and incremental velocity change to the spacecraft. 
The two thruster on the the positive yaw axis would create the moment required for the Z axis. 
While this simple design is effective it does not provide the redundancy for rolling maneuvers, 
redundancy is provided for the delta v maneuver in that two of the four thrusters on the pitch 

face could be fired to create a translational movement 

The disadvantages to the eight thruster design are the complexity of two 
additional thrusters, the movement of the center of mass (slight as it may be), and the accuracy 
of the cant angle required. These disadvantages for the eight thruster system must be weighed 
against the major disadvantage for using only six thruster, no redundancy within the 
propulsion system. While the secondary missions of the propulsion system are themselves 
redundant to some other system, the primary missions listed for the propulsion system have no 
other means to complete the assigned tasks. A trade off between full redundancy and necessity 
to complete the primary objectives is required to be evaluated at this time and beyond the scope 
of this preliminary design. 


61 


Thrusters 



Figure 6.5 Thruster Placement Summary 


3. Propellant Requirements 

The amount of fuel required for maneuvers was anticipated to e sma 
due to the mission, mission life, and orbit. One of the goals of the des.gn ts to eep 
everything as simple as possible. Combining these two thoughts ,t was ee, e 
use monopropellant hydraeine (N 2 Ha) as the fuel. For the antic.pated steady 
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low thrust maneuvers, hydrazine has a specific impulse, I s , of 225 seconds. 
Calculations are included in Appendix G. 

a. Spin Down and Detumble 

The first maneuver is to despin the satellite for three axis stabilization. 
Typically, the second stage is spinning at 0.033 rpm upon spacecraft separation. A 
worst case of 5 rpm was designed for to include the possibility of tumbling upon 
separation (or later during the mission). If reorientation is required, the satellite 
would be despun to less than 0.05 rpm, then reoriented. This uses only a negligible 
amount of fuel - even if a worst case maneuver of 180 degrees is performed. 

b. Orbital Insertion 

The second maneuver is to correct for possible errors in launch 
vehicle orbital insertion. The worst case of 18.5 kilometers from the expected orbit 
was designed for. The orbital section shows that the delta v required is 10.238 m/s 
(includes apogee and perigee firings). To do this maneuver four 1 lbf thrusters 
were used canted at 35.8 degrees from the thrust axis. 

c. Atmospheric Drag Corrections 

The next important consideration is correcting for atmospheric drag. The 
satellite must remain within 1 nmi (1.852 km) of its intended altitude. Drag 
calculations are included in the orbital section (Section III). The result is a delta v 
of 8. 172 m/s. 

d. Station Keeping 

The satellite must remain within ± 0.5 degrees of a 70 degree 
inclination. The Delta II inserts the satellite to within ± 0.05 degrees. The change 
in inclination due to effects from the earth, moon, and sun over the three year life 
are ± 0.0942 degrees (calculations in the orbital section). Therefore the worst case 
change in inclination with respect to time is ± 0.1492 degrees. This is well within 
the allowable tolerances so no extra fuel is necessary for station keeping. 
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e. Deorbit 

M the satellite's mission terminarion a decision is ,o be made regarding th 
disposition - dte spaeecta, Previous,, ihe common P^ce was to leave Ute spaeee^ 
dormant in their original orbits, considering them only as space ju . ow 
exploitation starting its fourth decade the heavens are begmnmg to ftllw.th ese e p 
satellites An investigation was conducted into the most practical method to d.sp 

satellite at the end of its lifetime. Three concepts governed the analysis; ^ 

back ta to the Earths atmosphete ,0 he burned up on teentry. second, minimize the ,m ac, the 

A controlled reentry into the Earths atmosphere is defined as: upon command to 

terminate the mission dte spacecraft would conduct a single deorbit bum bringing the * 

within 50 to 00 hm perigee altitude. This would enstne atmospheric capture amenable 
accurate placement of the geographic point in which dte destructton won occu ^ 

meet diese requirements die spacecraft would have to develop a !20.5 m/s de ta v a e 
To develop this much delta v the spacecraft would require approximately 75 kilograms 
fuel on board, this is five rimes the mission's fuel requirement. While thts tdea mee 
and third concepts of spacecraft disposal, it would require dte addmon of two addmona ue 

tanks and ultimately a larger spacecraft bus. . 

An uncontrolled reenny of the spacecraft would be defined as: a deorb, t of the 
spacecraft to a lower altitude and the satellite would complete mom than one orb,, beforo decay 
of the orbit was sufficient enough to cause die destruction of the satellite. Thts idea could be 
carried out with minimal impact to the mission requirement and while it is slower to achieve 
satellite destruction it is an otder of magnitude faster than leaving the satellite t„ te ongitra 
orbit. This idea would allow for a single deorbi. bum a. the end of mission . e, 
would expend all of the remaining fuel on board to achieve the highest delta v poss, «• 
to launch die fuel tank within the satellite would be filled to capacity, rather than i ingony 
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maximum mission requirement, this would give approximately 12.5 kg of fuel mass margin. 
Any fuel not used from the allotted maximum fuel requirement would only add to the fuel 
available at satellite end of life. Using the ASAP program the orbit was input for a satellite 
with a 509.3 km apogee and 430 km perigee, representing an end of life deorbit bum of 12.5 
kg of fuel. The ASAP program computed a decay of the orbit and destruction of the satellite m 
less than 9 years. This time frame appears to be very reasonable for a satellite of this size and 
weight. The LDEF satellite, which was much larger and heavier, was launched into an orbit 20 
nm less than this proposed orbit and its projected decay time to destruction was approximated 
to be 7 years. The idea of an uncontrolled reentry using the remaining on board fuel appears to 
be the better of the two concepts for destruction of the satellite. 

/. Total Fuel Required 

Since the worst case was designed for, the propulsion system is in 
reality a backup for the attitude control system. This provides a valuable 
redundancy. The total propellant required is calculated in table 6.6. 
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Table 6.6 Propellant Budget 


1. Propellant for delta v maneu vers and control 

2. Allowance for off-nominal performance 

3. Allowance for off-nominal op erations 

4. Mission margin (reserves) 

5. Contingency 

6. Total required propellant 

7. Residual propellant (trapped in motor case, tanks, lines, e tcQ 

8. Loading uncertainty 

9. Total propellant load 


11.03 
0.11 
0.11 
1.10 
1.10 
13.45 
0.27 
0.07 
13.79 kg 


4. Propellant/Pressurant Tank Selection 

To decide on a fuel tank from those available commercially, the required 

tank diameter was first determined to be 0.322 meters (see Appendix G). 

A survey was conducted of off-the-shelf positive expulsion tanks based on 

the the above minimum fuel tank diameter and the diameter of the satellite cylinder 
(s 0.5 meters) . The TRW 80225-1 sphere used by the OTS-Marex program was 
selected since it was the closes, match in propellant capacity. Several other tanks 
were in the same range. Table 6.7 describes the features of this tank. 
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Table 6.7 Propellant/Pressurant Tank Characteristics 


Expulsion Device 

AF-E-322 Diaphragm 

Volume (m 3 ) 

0.0306 

Propellant Capacity (kg) 

24.1 

Outside Diameter (m) 

0.39 

Minimum Wall Thickness (mm) 

0.48 

Weight (kg) 

3.70 

Pressure (P 0 ) (MPa) 

2.20 

Pressure (P D ) (MPa) 

3.52 

Pressure (Pb) (MPa) 

4.70 


5. Subsystem Operations 

Figure 6.5 is a schematic of the propulsion subsystem. This schematic 
shows a single tank connected to eight thrusters. The propulsion subsystem 
includes a fill/drain valve for the pressurant and one for the propellant. These will 
be manually operated during prelaunch operations. Two pressure transducers - 
one for each tank exit line - are included to measure absolute pressures. One 
pressure regulator for each line is included. Also a filter is included near each tank 
exit to remove impurities from the lines before they reach the thrusters. As a 
safety feature, each thruster has an isolation valve in its line in case of line breaks 
or failures. Also note that the propellant regulation valve is an integral part of each 

thruster. 
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D. MASS AND POWER SUMMARY 

Each thruster requires 9 watts during its operation. Since no more than four 
thrusters will be used at any one time, the power requirements will never exceed 
36 watts. Also the use of thrusters will only occur during the standby mode so only 
a minimal power of approximately 1 watt will be used during payload operation. 
Table 6.8 shows a breakout of the propulsion system masses. 


Table 6.8 Propulsion System Mass Summary (kg) 


1 lbf thrusters (8) 

2.65 

Propellant/Pressurant Tank 

3.70 

Tubim^^ 

1.0 

Valves, Filters, & Pressure Transducers 

1.0 

Propellant/Pressurant 

13.89 

Total 

22.24 
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VII. ATTITUDE CONTROL SUBSYSTEM 


A. MISSION REQUIREMENTS 

The Attitude Control System (ACS) is designed to control Spacecraft 
orientation, either automatically or on command, in the presence of external 
disturbances by using torque producing devices. The equations of motion are 
analyzed by the designer and the external disturbances modeled in order to select a 
proper control system to manage the spacecraft angular momentum. 

Listed below are the design specifications as provided by the statement of 

work: 

(a) Three axis stability with slew capability for sensing and standby modes, 

(b) Attitude accuracy: ± 0.5° / axis (3a), 

(c) Maximum time to slew: 90 degrees in 15 minutes, 

(d) Slew settling time: 1 minute, 

(e) Rate stability: 0.003° / second (3a) per axis, 

A payload supplied Inertial Measurement Unit (IMU) is available for short 
term attitude reference. It is aligned by a pair of star sensors contained within the 

payload. 

Further requirements are added due to the choice of satellite configuration and 
the power collection method. The design selection of fixed solar arrays requires the 
ACS to operate in multiple modes to collect power and to point the payload for 
sensing. It must also be able to slew between these modes accurately and with 
adequate speed. Intermediate modes include acquisition/de-tumble and boost 
maneuvers. Mass and power are kept to the minimum required within the 
constraints, specifications, and component availability. The majority of the satellite 


69 



lifetime is spent in a relatively quiet mode with respect to the ACS, but it must be 
extremely reliable for mission success. 

B. ACS CONCEPT 

1. Configuration and Operation 

The choice of configuration was driven mainly by the requirements to be 
able to slew about multiple axes in varying geometries and to maintain reasonably 
accurate three axis stability. Adequate redundancy is crucial because the ACS 
must be able to orient for both payload operations and power collection. Three 
orthogonal Reaction Wheel Assemblies (RWAs), one along each body axis, allows 
somewhat independent control of rotation about each axis. A fourth RWA is 
skewed 45 degrees out of plane with respect to the others to back up any single 
RWA failure. In the event of a RWA failure, the redundant rate will be commanded 
V3 times the normal rate to achieve the desired affect while the undesired 
components are temporarily taken up by the other wheels. Location of the major 
hardware is depicted in Figure 7.1 . 

Momentum dumping is accomplished by pairs of orthogonal magnetic torque 
rods. Normal operations call for the pairs to work simultaneously to rapidly 
desaturate the RWAs before slew maneuvers and to periodically activate to keep 
the RWA bias low. Specifically, automatic desaturation begins at 210 rpm. This 
minimizes dynamic coupling in the euler equations which are located in Appendix 
H. The rods can work singly in the event of a failure. Redundancy for momentum 
dumping is provided by the propulsion system, available in the event of complete 
torque rod failure. Both the magnetic dumping system and the propulsion system 
give limited three axis stable capability in the event of multiple RWA failure. The 
control system concept is depicted in Figure 7.2. Its algorithm contains all 
described failure modes. 
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IMU: 

Inertial Position 
and Rate for 
Slew Mode 


Earth Sensor: 

Roll, Pitch 
Errors for 
Pointing Mode 


Gyro 

Assembly: 

Body Rates for 
Pointing Mode 



(Command) 


(Tachometer) 


Torques to 



ACC #1: 

-Primary Sense 
Mode Processing 
-Assumes Both 
if #2 Fails 


ACC #2 

-Primary 
S untrack Mode 
Processing 
-Assumes Both 
if #1 Fails 



Pitch RWA 




Spacecraft 

Body 


Torques to 


Spacecraft 

Body 


Torques to 



Spacecraft 
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(Momentum 

Dumping)^ 


Sun Sensor: 


1 


Yaw Errors and 




Sun Alignment 





(back up Mom. 



Spacecraft 

Body 
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Figure 7.2 ACS Concept Block Diagram 
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Attitude errors are induced by solar, aerodynamic, gravity gradient, 
magnetic disturbance torques and perturbations during desaturation and Av 
maneuvers. The two Attitude Control Computers (ACC) receive data from the 
IMU, the earth sensor, one of two sun sensors in view, and the back up gyro 
assembly to compute and store two sets of these errors: (1) the euler angles with 
respect to the standard nadir pointing sensing coordinates, and (2) the euler angles 
with respect to the suntracking coordinates which track the incoming sun vector 
and can be considered "inertially" fixed. From these, each computer can calculate 
the direction cosine matrix used in its duty slew direction. In the event of a single 
computer failure, the other assumes the load for both. The first set of errors include 
the orbital rate while the second set are fixed with respect to the orbit normal 
coordinate system. Twelve independent transformations exist in each case ( sign 
ambiguities are removed ), and the computer defaults to the one which results in 
the minimum total correction path (or slew path) but any specific direction cosine 
matrix can be chosen. The duty ACC then commands the RWAs to perform the 
chosen sequence of single axis slews to reach the target axes (i.e. to zero the euler 
angles ). From here, smaller slews can be commanded to accomplish offset nadir 
pointing or to correct a thermal problem. Single axis slew sequences are not the 
fastest method but they simplify constraint checking, allow separate orthogonal 
error computation, and minimize dynamic coupling. A more detailed description of 
the coordinate systems used is in Appendix H. 

2. Operating Modes 

The spacecraft is deployed from the launch vehicle at or near desired orbit 
parameters with a spin of 0.2 RPM about the designated body Z- axis. Given the 
inertia properties of chapter II, this spin will remain stable in the absence of major 
dissipators. The propulsion system will slowly despin the satellite, if necessary, until 
a sun sensor acquires the sun. At this time the sun vector is used to control a slow 
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slew to roughly point the solar array face using the thrusters. At steady state, the 
solar arrays deploy, all components of the ACS power up and star data is linked to 
begin alignment of the IMU. The entire process is autonomous with the maneuvers 
time tagged shortly before launch. The detumble mode works in an identical 
manner. The RWAs shut down when the body rates indicate that for some reason 
three axis stability is lost and the last known sun position is used for search. If this 
doesn't work the ground controoler must link the data. The acquisition sequence 

then repeats to regain control. 

Seven to eight orbits will be available to power up and self test all 
subsystems due to the small magnitude of the disturbance torques. When ready, the 
ACS begins a slow slew to align with suntracking coordinates. The majority of the 
lifetime will be spent in this mode, tracking the sun to irradiate the solar arrays and 

maintain thermal balance. 

After successful power up, battery charge, and test of all systems the 
ACS is ready to perform slew maneuvers. When the command is received to 
sense, the ACS dumps any accumulated angular momentum to less than 2 RPM on 
each RWA and the commanded sequence of axis slews is performed to align with 
the sensing coordinates using the first set of errors. This maneuver may be time 
tagged if out of view of a controller. If offset nadir pointing is desired, a slow slew 
about the roll axis gives the commanded grazing angle. In this mode, noise effects 
suppression filtering is used to counteract sensor noise to provide the rate stability. 
Noise through motor vibrations and aeroelastic effects is damped by the control 
system. After a maximum of two orbits of sensing, the second set of errors is used 

to perform the slew back to suntracking. 

During momentum dumping, the ACS holds its current mode and the 

RWAs are given a command to a small pre-bias to ensure near zero RPM on 
completion of the process. All Av maneuvers are performed in the sun tracking 
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mode using the required offsets for thruster alignment and the ACS maintains an 
inertially fixed posture during bums. 

C. ACS DESIGN 

1. Component Sizing and Selection 

The equations of motion (EOM) for this attitude control configuration are 
well known and together with the specifications and the estimated disturbance 
torques, they are used to determine the required capabilities and sizes of the 
various components. Sensors must provide the required accuracy and no.se levels. 
RWAs and torque rods must supply adequate torque at rapid response, and 
electronics must have sufficient speed and memory. Off-the-shelf components that 
satisfy these findings are then chosen and used. The sizing and selection process, 
EOM, disturbance torque estimation ( all worst case scenarios ), and calculations 
are given in Appendix H. The results are summarized in Table 7.1 along with the 
mass and power consumption of each component. 


2. Control System Design 

The design of two primary controllers covers the operation of all the modes 
and will be discussed here: (1) the sensing mode controller, and (2) the slew 
controller. Only subtle changes separate the controllers between suntracking and 
sensing mode: suntracking and Av mode; slew and acquisition. Detailed analysis: 
assumptions, EOM reduction, controller type and gain selection is given in 
Appendix H. and the results are summarized below. 
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Table 7.1 (a) ACS Component Sizing and Selection Sensors and 

Electronics 


Component 

Mass 

(kg) 

Power (W) 

Manufacturer 

Back Up Spring Restraint 
frvrn Assembly 

13 

19.5 

Heritage: DMSP 

Precision Pointing Earth 
Sensor 

Attitude Control Computers 

3.8 

4.0 

Barnes 

2.5ea 

6.0ea 

MIL STD 1750 (UPS 
Version) 

(2) — 
Cnn Sensor ( Sense Mode) 

0.04 

1.0 

Adcole 

Sun Sensor (Sun-Track 

0.04 

1.0 

Adcole 


Mode) 


Table 7.1 (b) ACS Component Sizing and Selection Torque /Angular 

Momentum Devices 


Component 

i 


Roll RWA 


Pitch RWA 


Yaw RWA 


Redundant RWA 


Torque Rods (6) 


Storage Peak Mass 
Capacit Torque (kg) 
y (Nms) (Nm) 


19.9 I 0.3 


Power 

i(W) 



Manufacturer 

Honeywell 


N/A 


.003 
max @ 
lOAm^ 


1.76ea 


Honeywell 


Honeywell 


Honeywell 


1.6ea llthaco 
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In the sensing mode the attitude errors are limited to a small range leading 
to linearized and uncoupled EOM. Each axis can have its own proportional and 
rate feedback controller ( PD ). Using the disturbance torque estimates the gains 
are selected to provide the required pointing accuracy. Using noise models the rate 
stability characteristics are analyzed. The main source of noise is small and random 
structural mode vibration excited by the gimbaled transmit antenna which tracks 
during sensing. Inertia effects from aeroelasticity, and RWA motors plus noisy 
references from the sensors enter the control system. A circuit to detect jitter out 
of limits is included to signal a gain change adjusting controller bandwidth to 
suppress the noise effects. Control parameters for each axis are summarized m 
Table 7.2 

Table 7.2 Controller Parameters for Pointing Modes 


Loop 

Gain 

(NM/rad) 

Time 

Constant (s) 

Control Bandwidth by 
Cutoff Freq. 

roll 

46.37 

5.82 

0.0546 Hz 

pitch 

38.896 1 

6.395 

0.0498 Hz 

yaw 

26.45 1 

7.491 

0.0425 Hz 


For large angle slews the EOM are no longer linear and are doubly coupled: 
(1) significantly by the off-axis slew rates, and (2) weakly by small off-axis 
momenta stored in the RWAs ( not completely dumped prior to slew ). Given the 
appropriate set of errors for the slew direction, the computer transposes the 
direction cosine matrix to determine the target change in angle for each axis to 
align the body axes with the desired axes. During each single axis portion of the 
sequence, the controller exploits the independent axis capability and commands 
the off-axis rates to zero to minimize the time that coupling affects the single axis 
rotation . 
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The scheme for each axis is a nonlinear saturation feedback controller with 
variable gain on the rate depending on the initial errors. It is described by: 

^command = ■ N sgn ( position + gain *rate ) , 

where N is the selected torque level. Gain selection is also based on the axis 
( primary or off-axis ). Torque switching levels are also variable based on initial 
error magnitude and entry /exit torques are shaped to prevent structural excitation 
and actuator wear. Estimates from the structural team (chapter IX) indicate the 
fundamental frequency around 22-25 Hz, relatively rigid interms of the low 
operating torques and slew rates of the ACS. Enough so that rigid body 
assumptions are made for this preliminary design. 

A capture box is established at 2 degrees and 0.05 degrees per second 
signaling the slew controller transition to the pointing controller. This prevents 
chatter and actuator wear. The selection and mathematical models for this 
controller are given in Appendix H. 

D. PERFORMANCE PREDICTION 
1. Pointing Modes 

The well known EOM and designed controllers are simulated for various 
conditions using MATLAB and its associated SIMULAB for the MAC Ilsi. These 
programs are located in Appendix H. Additional results which may be of interest to 
the reader are also located in Appendix H, such as disturbance torques for various 
modes, off-axis responses during slew, phase plane plots, and wheel speed 
histories. 

The time responses to the estimated disturbance torques for both the 
sensing and suntracking modes are shown in figure 7.3. This is a worst case 
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scenario i.e. maximum disturbance torques which occur during RWA desaturation 
and a certain orbit configuration. The pointing accuracy specifications of less than 
one half a degree per axis are satisfied by these predictions. 



79 



Pointing Error, deg. 



The rate stability characteristics are shown in Figure 7.4 for each axis near 
its noise threshold i. e. that maximum level still within specs. A typical model for 
sensors is exponentially correlated noise. The threshold levels for this noise is 
approximately o = 67 rad/sec in the PSD which corresponds to 10* 4 degree rms 
magnitude in the time domain. Torque noise is modeled by white noise at low 
intensity. These level far exceed the normal levels encountered so the rate stability 
specs are met for this preliminary design. 
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Figure 7.4(b) Pitch Rate Stability Response Near Limit 


xlO-3 Yaw Axis Rate Stability Response 



Figure 7.4(c) Yaw Rate Stability Response Near Limit 
2. Sensing Modes 

A typical single axis slew is depicted in Figure 7.5. Slew rates of 
approximately one degree per second are attainable well within the required rate. 
Slower rates can be used and will consume less power. The time response shows 
nice damping characteristics. Once the spacecraft is slewed into the capture box, 
the individual axis PD controllers take over to settle the slew. The settling time is 
approximately twenty-five seconds, well within the requirements. 
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Time from Capture, sec. 


Figure 7.5(b) Capture Initial Conditions; 2 degrees and 0.5 dps 

As mentioned previously the ACS has redundancy exhibits excellent 
performance within the specifications in this preliminary design. The degraded 
modes all perform favorably but do not meet the normal specs. They can, however, 
maintain an operating satellite. 

The rate stability requirement is key and must be validated by a more 
advanced model. Further testing on subsequent design should focus here because 
from this report, it appears the other requirements are easily met. In addition, the 
slew rate and pointing accuracies far exceed the specs so that reduction in 
available control effort can be investigated. In its present condition, however, this 
bus could slew and control a much heavier payload. 


84 



VIII. 


THERMAL CONTROL 


A. MISSION REQUIREMENTS 

A spacecraft thermal control subsystem provides for adequate conduction and 
radiation heat paths between all spacecraft elements so that temperature sensitive 
components remain within specified temperature limits during all mission phases. 
This generic definition applies to this spacecraft bus as well as the cryogenic 
payload, however the bus thermal control subsystem is designed to function 
independently from that of the payload. The payload is equipped to provide its own 
cooling and heating for both sensing and stand-by modes. The only payload 
imposed thermal constraint on the bus was to maintain the bus/payload interface 
temperature between 0 and 30 degree Celsius. No bus components are required 
to maintain cryogenic temperatures, and in fact some bus components may go 
above the 30 degree Celsius but the interface between the bus and payload will 
always be in the required range. Specific temperature range requirements for 
selected bus component and subsystems are given as Table 8.1. The references 
for the table are given in Appendix I. 

B. THERMAL ENVIRONMENT 

Spacecraft orientation and configuration will vary throughout the lifetime of the 
system. The vehicle will go through various transfer, acquisition, and slewing 
modes. Solar arrays and the antenna will be deployed early in the mission. 
Thermal analysis, however, is done only for the two primary mission modes; on 
orbit stand-by and earth sensing. Theses two mode also represent the two 
extreme thermal cases for the system. 
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Table 8.1 Component Allowabel Temperature Ranges 


COMPONENT / REGION 

ALLOWABLE 
TEMPERATURES 
MIN / MAX ( *C ) 

PAYLOAD/BUS INTERFACE 

0/+30 

BATTERIES 

0/+10 

FUEL TANK 

-10/ + 50 

THRUSTERS/PIPING/ 

VALVES 

-10/ +50 

SOLAR ARRAYS 

-160/ +80 

ANTENNA 

-170/ +90 

TT&C ELECTRONICS 

-55/ +125 

EPS ELECTRONICS 

-55/ +125 

ADCS ELECTRONICS 

-55/ +125 

REACTION WHEELS 

-45 / + 55 

SENSORS 

-30 / +55 


The hot case is defined as the earth sensing mode during which the payload 
is acquiring and transmitting data. In this mode the high data rate transmitter is 
active and the batteries discharge providing their greatest thermal output. Other 
TT&C and ADCS components operate in this mode intermittently. The spacecraft 
is earth pointing with the payload side facing earth. With an orbit altitude of 275 
NM, inclination of 70 deg., and longitude of ascending node at 90 deg, we find that 
the spacecraft does not experience an eclipse period during this orbit. 

The cold case is defined as the standby mode during which the spacecraft is 
sun pointing (solar arrays are normal to sun vector). The batteries are being 
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charged in this mode and the high data rate transmitter is non-operational. The low 
data rate transmitter and a variety of other components will operate intermittently. 
The orbit altitude and inclination are the same as in the hot case, however with the 
longitude of the ascending node at 0 degrees the satellite experiences a 37 minute 

eclipse period per orbit. 

The significant thermal environment inputs of the system are given as Table 
8.2. A number of smaller electrical components are not identified 

Table 8.2 Spacecraft Bus Thermal Environment 


ELEMENT 

HOT CASE 
(EARTH SENSING) 

COLD CASE 
(STAND-BY) 

SOLAR INPUT 

1390 W/m 2 

1310 W/m 2 

EARTH REFLECTED IR 

76 W/m 2 

72 W/m 2 

ALBEDO 

325 W/m 2 

306 W/m 2 

BATTERIES 

27 W 


HIGH DATA RATE 
XMTR 

75 W 

5 W 

LOW DATA RATE 
XMTR 

5 W 

22 W 

GTU RFCS 

17 W 

17 W 

ASC 

25 W 

25 W 

RIU 

6 W 

P 25 W 

REACTION WHEELS/ 
TORQUE RODS 

20 W 

20 W 

GYRO ASSEMBLY 

19 W 

19 W 

UPLINK PROCESSOR 

9.5 W 

9.5 W 

DOWNLINK 

PROCESSOR 

9.5 W 

9.5 W 
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explicitly but are modeled later as an aggregate. Due to the 275 NM altitude the 
effects of aerodynamic heating have been ignored. Due to the location of nozzles 
plume impingement by the spacecraft propulsive elements is also ignored. 

C. THERMAL CONTROL SYSTEM DESIGN CONCEPT 

Thermal control of the spacecraft bus is achieved using a near passive thermal 
control system. Passive control systems are simpler, more reliable, and more cost 
effective than active systems. The most attractive of these features are the 
simplicity and reliability. It was decided that the majority of the mission risk should 
be with the experimental payload, therefore all of the bus subsystems are as simple 
and reliable as possible. Passive components of the system include Optical Solar 
Reflectors (OSR), Multi Layer Insulation (MLI), temperature sensors, fiberglass 
standoffs, and space qualified coatings/paints. Temperature sensors are used to 
control strip heaters located on critical elements. A pair of temperature sensor 
controlled louvers is used for thermal control of the batteries. All items are readily 
available and generally "off the shelf" hardware. A thermal control subsystem 
schematic is shown in Figure 8.1. 

The spacecraft bus radiates excess heat to cold space by way of OSR 
radiators. OSR's are bonded to the aluminum honeycomb outboard panels. They 
have a clear view to space and effectively dissipate heat since they have high (0.8) 
emissivity, while absorbing little incident energy due to low absorptivity (0.12). 
Heat dissipating hardware is mounted directly to the inside face of the honeycomb 
panels. A good conductive path is maintained by using thermal grease at the box 
to panel interface. The sizing of the radiators is based on the hot case. 

During the cold situations critical components are kept above their low 
temperature limits by applying heater power. The bus is equipped with 15 
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Figure 8.1 Thermal Control System 

heaters and range in size from 2.5 W to 20 W watts. Spacecraft heaters are 
electrical resistance elements which are mounted directly to the exterior surface of 
the critical component. The heaters are controlled by a closed loop system which 
includes local temperature sensors. Although the heaters function autonomously 
the provision has been made for ground controlled operation and adequate 
telemetry channels are available. 

The Nickel/Hydrogen batteries used in our spacecraft have a very narrow 
operating range (0-10 deg C). Based on preliminary hand analyses a simple 
OSR/heater system was deemed feasible. After further computer modeling and 
optimizing the batteries were placed on the ±X panels and the need for louvers for 
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was demonstrated. These active elements provide for a variable 
absorptivity/emissivity ratio which, in combination with battery mounted strip 
heaters, will ensure that the batteries will remain in the appropriate range (see 
analysis results in section D). Additionally the batteries are thermally decoupled 
from the rest of the spacecraft to eliminate spurious heat inputs from other 
components. The mounting of the battery, using thermal stand-offs, is shown as 
Figure 8.2. 



A similar isolation concept is applied at the payload/bus interface to isolate 
the payload. A detail of the interface is given as Figure 8.3. 


MLI BUS ATTACH PLANE 



THERMAL STANDOFF PAYLOAD ATTACH PLANE 

Figure 8.3 Bus Payload Isolation Schematic 
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To provide additional heat to cold sensitive components electrical strip 
heaters are applied directly to the surface of these components. Heaters are 
closed loop controlled by selected temperature sensors. Downlink telemetry is 
provided for critical temperature sensors while uplinks for the heaters and louvers 

provide for ground control of these components. 

A full summary of all thermal hardware is given as Table 8.3. All properties 
listed are BOL, the models include degraded values based on mission life and other 
factors. 

Table 8.3 Thermal Hardware Summary 


COMPONENT 

IMPORTANT CHARACTERISTICS 

HEATERS 

- Electrical strip heaters 

- Temperature sensor controlled 

- Ground control available 

- "Off the shelf', industry standard equipment 

- Total heater count of 120 W 

- Locations: 

Batteries / 40W (20 W per battery) 

Antenna drive / 20W 

Launch vehicle separation pyro / 5W 

Propellant tank / 20W 

Propellant lines and valves / 15W (total) 

Thrusters / 20W (2.5W each) 

LOUVERS 

- Only for battery OSR's 

- Temperature sensor controlled 

- Ground control available 

- Special made hardware for this mission 

- Developed and flight proven technology 

- 20 x 25 in. area covered 

- Surface properties when closed: 

absorptivity = 0.36 
emissivity = 0.45 

- Surface properties when open (underlying OSR 
properties): 

absorptivity =0.08 
emissivity = 0.8 
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OSR 


- Body mounted (exterior honeycomb panels) 

- "Off the shelf, industry standard equipment 

- Applied in 1 x 1 in. squares, allow for variable geometry 

- Roll panels (±X): 

20 x 25 in. area on each panel, used for 
battery heat rejection 

- Pitch panels (±Y): 

net 20 x 55 in. area on each panel, used for 
electronics boxes heat rejection 

- Surface properties: 

absorptivity = 0.08 
emissivity = 0.8 


MLI 


- Mounted on virtually all exposed areas: 

honeycomb panels, aft face of antenna, 
antenna support structure, bus/payload 
interface, used to closeout the anti-earth 
side of bus, propellant lines and thrusters 

- 10 to 15 layer construction 

- 1 to 2 mil aluminized kapton inner/outer layers 

- Dacron spacer layers 

- "Off the shelf", industry standard equipment 

- Surface properties: 

absorptivity = 0.35 

emissivity = 0.53 


TEMP. SENSORS 


- Used for control of heaters and louvers 

- "Off the shelf, industry standard equipment, available in 
a variety of temperature ranges. 


STAND-OFFS 


- Used to isolate batteries from bus, and payload from bus 

- G-10 of generic fiberglass construction 

- Not an "Off the shelf item but very easy manufacture 


COATINGS 


- A variety of anodyzes and paints used to tailor surface 
properties where required 


D. THERMAL ANALYSIS 

The thermal analysis done in support of this design is intended to be used as 
a preliminary sizing of hardware and an initial validation of the thermal control 
concept. 
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1. Preliminary Analysis 

Preliminary thermal hand analyses were done in the following areas: 1) 
Battery temperatures , and 2) Solar array temperatures. The results of these 
analyses, done early in the design process, indicated that a passive thermal control 
system was feasible. The important results were the following: 

Solar array temperature extremes: 

Hot 40 deg C 
Cold -44 deg C 

Solar cell performance is not degraded within these temperature limits. 

Battery temperature extremes: 

Hot 10 deg C 
Cold 0 deg C 

These limits are set as the required operating temperature range of the batteries by 
the manufacturer. To maintain these extremes each battery is required to have a 
dedicated 0.25 sq. meter OSR radiator for the hot case and 21 watts of heater 
power. The complete analyses and all simplifying assumptions are included in 
Appendix I. 

2. Thermal Analysis Modeling 

A detailed lumped-parameter analytical model was developed using the PC- 
ITAS (Version 7.0) thermal analysis software developed by ANALYTIX, Inc.. 
The software is a PC based, modular, menu driven finite difference program which 
is capable of interfacing with other industry standard codes such as SINDA and 
TRASYS. 

The model developed for the spacecraft consisted of a system of thermally 
coupled nodes which represent significant spacecraft components. The model is 
essentially divided into two submodels. The first is developed by the geometry 
definition module within PC-ITAS and is a geometric representation of the 
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spacecraft exterior. The interior of the spacecraft is modeled by the user as a set 
of nodes and associated conductances which represent the electronics boxes, 
batteries, fuel tanks, reaction wheels, etc., and the way in which they are mounted. 
The two submodels are then coupled using appropriate conductance terms. 
Thermal parameters such as optical properties and internal thermal dissipations are 
defined for all nodes in separate modules. Orbital parameters are set and then the 
entire model is submitted for processing to the internal analysis module. The model 
developed for this spacecraft is shown in Figure 8.4. 

The model used in our analyses consisted of 80 active nodes. Detailed 

model information is provided in Appendix I 

A full thermal analysis is carried out for only two orientations, the earth 
sensing and the on orbit stand-by modes. The spacecraft configuration is the same 
during both cases, antenna and solar arrays deployed. 

Based on current thruster placement plume impingement analysis is 
considered unnecessary at this time. Aerodynamic heating was not a concern at 
the 275 NM altitude. 
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Figure 8.4 Thermal Model 
3. Results of Detailed Analyses 

The results of the detailed model are given in Table 8.4. The temperature 
extremes experienced by specific components are given for the earth sensing and 
stand-by modes. 
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Table 8.4 Analysis Results 


COMPONENT 

EARTH SENSING 

STAND-BY MODE 


MIN /MAX CC) 

MIN / MAX ( *C ) 

PAYLOAD/BUS/ 



INTERFACE 

+3/+11 

0/+30 

BATTERIES 

0/+13 

0/+18 

FUEL TANK 

-11/+10 

-91+29 

SOLAR ARRAYS 

-57/+13 

-44/+66 

ANTENNA 

-14/- 13 

-A3 1-21 

ADCS ELECTRONICS 

-40/- 10 

-39/- 10 

TT&C ELECTRONICS 

-49/+91 

-46/- 14 

EPS ELECTRONICS 

-37/+33 

- 46/+4 

REACTION WHEELS 

-40/- 10 

-39/- 10 


4. Conclusions 


The results of the current models show that critical hardware is maintained 
within required temperature limits. 

The batteries are the most temperature sensitive components, and 
consequently required the most detailed work. After sizing the radiators for the 
sensing mode, which was the hot case for the batteries, the batteries were running 
too cold in stand-by mode. The two possible solutions were to include very large 
heaters (up to 60W per battery) or to have a louver system which would allow for 
different surface properties for the hot and cold cases. To keep the power 
requirements lower the louvers are chosen. 

The payload/bus interface is controlled by the payload. For early versions of 
the thermal model the payload was modelled using knowledge of it's exterior and 
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the electrical power drawn. Modeled in this manner the interface temperature was 
very cold, -40 deg C. After discussions with NRL we were told to treat the 
payload as a temperature controlled node, which meant that the ITAS program 
would not allow the node to make excursions beyond the 0 to 30 deg C limits. The 
mechanism by which this is accomplished is internal to the payload, and is 

therefore not a primary concern of the bus design. 

In general the bus electronics are running cold, but within limits. This 
problem would be alleviated by reducing the size of the radiators. Due to the 
coarseness of the model, and time constraints, the electronics boxes mounted to 
the honeycomb panels were lumped together and attached to the OSR s, modeled 
at 20 x 30 in areas. A finer detailed model could separate each box and its 
associated radiator patch. The overall radiator size would area would 
consequently go down, thereby the operating temperature ranges would drift 

upward. 

The solar arrays experience the greatest temperature fluctuations of the 
entire system, however they still remain within their required limits for optimal 
operation. 

Some temperature sensitive components such as thruster, sensors, and tank 
piping/valves are not explicitly defined in the model. These components are too 
small to have an effect on the over all system temperates but should be modeled at 
a later date due to their sensitive nature. 

E. MASS AND POWER SUMMARY 

A mass and power summary for the thermal control subsystem is given as 
Table 8.5. The total power shows all the component requirements, however not all 
the heaters will be active at the same time. The highest active power requirement 
is 80 W and occurs during the stand-by mode. 
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Table 8.5. Thermal Control System Mass And Power Summary 


COMPONENT 

MASS (KG) 

POWER (W) 

MULTI LAYER INSUL. 

8.2 

- 

OPTICAL SOLAR REFL. 

1.8 

- 

COATINGS / PAINTS 

0.45 

- 

HEATERS 

23 

120.0 

TEMP. SENSORS 

1.4 

2.0 

LOUVERS 

5.0 

5.0 

ELECTRONICS 

3.6 

5.0 

STAND-OFFS 

2.0 

- 

TOTALS 

24.8 

132.0 


98 

























IX. STRUCTURAL DESIGN 


A. MISSION REQUIREMENTS 

The structural subsystem must be designed to withstand the dynamic loads 
of the Delta II launch vehicle. During launch, the Delta II experiences axial 
loads of 7.2 g’s and lateral loads of 2.5 g’s. To ensure that the spacecraft bus 
can adequately withstand these loads, a factor of safety of 1.5 is used. The 
natural fundamental frequencies of the launch vehicle are 35 Hz axially and 15 
Hz laterally. To prevent dynamic coupling between the launch vehicle and the 
spacecraft, the spacecraft's modal frequencies must be above the launch 
vehicle's natural frequency in the thrust axis. If payloads are unable to meet 
this criteria, the structural design must be coordinated closely with the MDAC 
Delta Program Office so that appropriate analyses can be performed. 

B. DESIGN OVERVIEW 

The structural configuration of the spacecraft appears as Figures 9.1 and 
9.2. The primary load carrying members for axial compressive loads and 
bending moments are the central cylinder and conical adapter. These are 
monocoque shells of 6061-T6 Aluminum designed to withstand the ultimate 
quasi-static loads during launch. The loads are assumed to be a factor of safety 
of 1.5 times the axial and lateral limit loads of 7.2 and 2.5 g's of the Delta II 
launch vehicle. Aluminum 6061-T6 was chosen for its ease of machining and 
large strength-to-weight ratio. There are seven panels made of honeycomb 
sandwich material with aluminum skins and core. The top panel supports the 
load of the payload and transmits it to the central cylinder and cone. The two 
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equipment panels support the subsystem components. The remaining panels 
make up the body of the spacecraft bus. The panels are designed for stiffness 
to meet the design criteria for a minimum natural frequency of 50 Hz and stress 
due to dynamic loads. The adapter ring transfers the loads from the conical 
section to the 6019 adapter of the Delta II. The frame of the spacecraft bus is 
made of aluminum square tubing. It provides the mounting surfaces for the side 
panels and the parabolic antenna. Dimensions for the various components 
appear in Table 9.1. 

1. Central Support Assembly 

The central support assembly consists of the cylinder and conical 
sections. These are designed to support the dynamic loads on the spacecraft 
during launch. The procedure used to design these components appear in 
Appendix J. 

2. Honeycomb Panels 

The honeycomb panels are designed for stiffness to meet a natural 
frequency of 50 Hz. This frequency was chosen to prevent resonance between 
the launch vehicle and the spacecraft. The panel which must support the 
largest mass is the one between the payload and the bus. The procedure used 
to find the core depth assuming a skin thickness of 1 mm appears in Appendix 
J. All other honeycomb panels are made from the same material so that costs 
could be reduced by buying in larger quantities. Calculations in Appendix J 
show that the core depth chosen is adequate for all equipment and side panels. 

C. FINITE ELEMENT MODEL 

A finite element model of the spacecraft bus was created on GIFTS 
(Graphics-Oriented Interactive Finite-element Transportable System Version 
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6.5.2). GIFTS is a program used to analyze and design structures for static and 
dynamic loads. It is available on the VAX system of the Computational 
Instruction Laboratory of the Aero/Astronautics Department of NPS. GIFTS 
obtains displacements, stresses, stress contours, bending moment and shear 
forces, and associated normal and shearing stresses for each element. Free 
vibration response mode shapes and frequencies and forced responses may 
also be obtained.. 

All components listed in Table 9.1 were modeled on GIFTS. The central 
support assembly and all honeycomb panels were modeled as membranes, 
while all the struts were modeled as hollow beams. The payload and all other 
subsystem components such as attitude control and electric power were 
modeled as distributed masses on the equipment panels. The parabolic antenna 
was modeled as a point mass. The data files used to create the model appear 
in Appendix J. 

1. Modal Analysis 

The first eight modal frequencies of the spacecraft are given in Table 
9.2. In the first mode, the central support assembly moves opposite to the box 
framework which supports the side panels. The second through fifth modes 
show bending and twisting of the antenna support. These modes are similar to 
the first and second modes of a cantilever beam. Because the second and third 
modes have frequencies which are close to the launch vehicle natural 
frequency, the antenna support structure should be stiffened to raise these 
frequencies. This should be considered in refinement of this design but will not 
be considered at this time. The sixth through eighth modes show warping of 
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the equipment shelves. These modes are similar to the first and second modes 
of a membrane. 

2. Deflections and Stresses 

A static analysis of the spacecraft under loads equal to those 
experienced during launch was also performed. The maximum translational 
deflection occurs at the tip of the antenna support. This deflection is only 1.2 
mm. The maximum rotational deflection of 2.2 mm also occurs at this point. 
The maximum stress in the spacecraft is 87.9% of the maximum yield stress 
using Von Mises yield criterion. This results in a margin of safety of 13.8%. 

D. MASS SUMMARY 

The individual masses of each structural component appear in Table 9.1. 
These total to 164.97 kg. Assuming fittings, brackets, and other attachment 
devices add an additional 25% mass and assuming another 10% margin for 
growth, the total mass of the spacecraft bus structural subsystem is 226.8 kg. 
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Table 9.1 Structural Dimensions (length - m; mass • kg) 
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Table 9.2 Modal Frequencies (Hz) 


Mode 

Frequency 

1 

22.76 

2 

31.72 

3 

32.14 

4 

40.31 

5 

40.98 

6 

67.40 

7 

70.18 

8 

70.63 



Figure 9.1 Adaptor Cone and Cylinder 
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Figure 9.2 Honeycomb Panels and Box Frame 
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X. TESTING 


A. INTRODUCTION 

Due to the high risk nature of the proposed payload, the decision was made 
to minimize risk and cost for the bus design. To this extent, simplicity of design 
and testing was paramount. To minimize cost a modofied proto-flight test plan 
has been developed. The United States Air Force Military Standard Test 
Requirement for Space Vehicles (MIL-STD-1540B) was adopted as the basic 
guideline for developing testing requirements and acceptance levels. 

B. TESTING PHILOSOPHY 

All spacecraft development programs include intensive testing of all 
hardware to ensure proper operation during all mission phases. Testing is 
divided between system level tests and component/subsystem level tests. 

For system level tests a dedicated qualification unit is usually constructed 
and used to qualify the design before the flight unit is completed. The 
qualification unit is built to the same demanding standards as the flight article. 
With production satellites this cost is usually acceptable, however with a unique, 
one of a kind, spacecraft the cost is excessive. Under these circumstances a 
"proto-flight" test plan is usually developed. This plan calls for building only one 
vehicle, to flight hardware standards, which is used in all system level testing 
and subsequently flown as the mission flight unit. The testing is usually done to 
levels which are lower than qualification levels but higher than flight levels. This 
standard industry practice is outlined in great detail in MIL-STD-1540B. 
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The deficiency of the proto-flight testing philosophy is that certain tests 
cannot be undertaken until the proto-flight unit is at least partially completed. 
This generally requires a great deal of time due to the strict standards of 
procuring, manufacturing, assembling, inspecting, and accepting flight hardware. 
For this reason it is proposed that early in the development of this system a full 
scale "test" model is developed for early development tests. 

The test model would include primary structure such as the central cylinder 
and cone, external honeycomb panels, internal decks, and internal framework. 
Additionally mass simulators of the significant bus components would be 
included, as well as a mass model of the payload. This unit could be easily 
constructed and used in dynamic testing such as acoustic or random vibration 
tests which not only are an indicator of the structural system performance but 
can be used to derive component test levels for critical equipment. With slight 
modifications, such as the addition of some heaters, multi-layer insulation, and 
optical solar reflectors, the test model could be used as a development system 
for the thermal control system which may be critical. 

A flowchart for the proposed testing is shown as Figure 10.1. 
Appropriate functional tests are performed before and after all environmental 
tests. The modified proto-flight testing philosophy speaks primarily to the 
system level tests. Component level testing is done as required to accept 
critical hardware. All mission specific electronics packages are to be tested on 
a component, and subsystem level. All other components will be off the shelf 
and flight proven. These components will be tested on a subsystem and space 
vehicle level only. All component and subsystem testing of the payload is done 
by the 
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DYNAMIC TESTS 
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Figure 10.1 System Level Test Flow Chart 
payload supplier. Only those tests required to ensure proper integration of the 
payload and bus are considered in the bus test plan. 

C. SYSTEM LEVEL TESTING 

System level tests are done on both the test model as well as the proto- 
flight model. Not all tests are done on both. 
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1. Proto-Flight Unit 
a. Functional Test 


Purpose: This test verifies that mechanical and electrical 

performance of the space vehicle meets specification requirements and 
validates all test techniques and software algorithms used in computer assisted 
commanding and processing. 

Mechanical Functional Test: Mechanical devices, valves, 

deployables, and separable entities shall be functionally tested, with the 
spacecraft in the launch, or orbital configuration appropriate to the function. 
The maximum and minimum operating limits shall be determined with respect 
to mechanics, time, and other applicable requirements. Where operation at 1-g 
environment cannot be performed, a suitable ground test fixture shall be utilized 
to permit operation and evaluation of the devices. Fit checks shall be made of 
the space vehicle physical interfaces with the launch vehicle by means of 
launch vehicle master gages or interface assemblies. 

Electrical Functional Test: The space vehicle shall be in its 
flight configuration with all components and subsystems connected except 
pyrotechnic elements. The operation of all thermally controlled components, 
such as heaters and thermostats, shall be verified by test. The test shall verify 
the integrity of all electrical circuits in the initiating stimulus and confirmation of 
the successful completion of the event. Wherever possible, equipment 
performance parameters (such as power, voltage, gain, frequency, command 
and data rates) shall be varied over specification ranges to demonstrate the 
performance margins. A segment of this test shall operate the space vehicle 
through a mission profile with all events occurring in actual flight sequence to 
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the extent practical, this sequence shall include final countdown, launch, 
ascent, separation from the launch vehicle, orbital injection, and orbital 

operation as appropriate. 

b. Acoustic Test 

Piirnose: This test demonstrates the ability of the space vehicle 
to withstand the design level acoustic environment which is the maximum level 
imposed during launch plus a design margin. This test also verifies the 

adequacy of component vibration qualification criteria. 

Test Description: The mass model used for this test shall be of 

identical structural quality and design as the flight vehicle. It shall be mounted 
to a flight-type support structure. The model shall represent the space vehicle 
in the stowed launch configuration. Adequate dynamic instrumentation shall be 
installed to measure vibration responses at attachment points of critical and 
representative components. 

Test revel? and Duration; The acoustic test spectrum shall 
be the design environment which is the maximum predicted flight environment 
plus the design margin (3dB). The overall sound pressure level of the test shall 
not be less than 141 dB. Exposure time shall be at least three times the 
expected flight exposure time to the maximum flight environment, but not less 
than 3 minutes. Operating time shall be divided equally between redundant 

circuits. 

c. Pyro Shock Test 

Pnrnose: This test demonstrates the capability of the space 

vehicle to withstand the design level pyro shock environments which are the 
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levels predicted for flight plus a design margin. This test also verifies the 
adequacy of component pyro shock criteria. 

Test Description; In this test all pyrotechnically operated 
devices shall be operated. Separation subsystem shocks are often more severe 
than those from other pyrotechnic devices, and operation of the separation 
subsystem is therefore especially significant. For this test the vehicle shall be 
suspended or otherwise supported to prevent the possibility of recontact 
between the separated portions thereof. The launch vehicle separation adapter 
shall be attached and appropriate pyrotechnics used to impose the required 
shock. Adequate dynamic instrumentation shall be installed to measure pyro 
shock responses in 3 axes at attachment points of critical and representative 
components. 

Test Levels and Duration; All pyrotechnic devices shall be 
fired at once. Those pyrotechnic devices producing levels within 6 dB of the 
maximum shock response measured from any of the devices shall be fired an 
additional two times to provide the expected variability in the shock 
environment. 

d. Pressure Test 

Purpose: This test demonstrates the capability of fluid subsystems to 
meet the flow, pressure, and leakage rate requirements specified. 

Test Description: T ests shall be performed to verify compatibility 
with the test setup and to ensure that proper control of the equipment and test 
functions provided. The requirements of the subsystem including flow, leakage, 
and regulation shall be measured while operating applicable valves, pumps, and 
motors. The flow checks shall verify that the plumbing configurations are 
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adequate. Checks for subsystem cleanliness, moisture levels, and pH shall also 
be made. 

In addition to the high pressure test, propellant tanks and thruster 
valves shall be tested for leakage under propellant servicing conditions. The 
system shall be evacuated to the internal pressure normally used for propellant 
loading and the systems pressure monitored for any indication of leakage. 

Test Levels an d Duration; The subsystem shall be pressurized to 
1.5 times the maximum predicted operating pressure and held there for 5 
minutes, then the pressure shall be reduced to the maximum operating pressure. 
This sequence shall be repeated three times. Inspection for leakage after these 
cycles shall be made at the maximum operating pressure. The duration of the 
evacuated propulsion system leak test shall not exceed the time that this 
condition is normally experienced during propellant loading. 

e. Thermal Vacuum Test 

Purpose: This test demonstrates the ability of the space vehicle 
to meet design requirements under vacuum conditions and temperature 
extremes which simulate those predicted for flight plus a design margin. 

Test Description: The space vehicle shall be placed in a 

thermal vacuum chamber and a functional test conducted to assure readiness 
for chamber closure. Equipment that does not operate during launch shall have 
power applied after the test pressure has been reached. A temperature cycle 
begins with the space vehicle at ambient temperature (23 +/- 10 deg C). The 
temperature is reduced to the specified low level and stabilized. Temperature 
stabilization has been achieved when the rate of temperature change is no more 
than 3 deg C/hour. Following the cold soak, the temperature shall be raised to 
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the highest specified temperature and stabilized. Following the high 
temperature soak, the space vehicle shall be returned to ambient temperatures 
to complete one temperature cycle. Functional tests shall be performed during 
the first and last temperature cycle at both the low and high temperature limits 
with functional operation and monitoring of perceptive parameters during all 
other cycles. In addition to the temperature cycles, the chamber shall be 
programmed through various orbital operations. Strategically placed 
temperature monitors shall assure attainment of temperature limits. 
Strategically placed witness plates and quartz crystal microbalances shall be 
installed in the test chamber to assure that outgassing from the space vehicle 
and test equipment does not degrade system performance beyond specified 
limits. 

Test Levels and Duration; The temperature extremes shall be 
established by a survey of predicted temperatures in various equipment areas. 
Temperatures on the components shall not be allowed to exceed the design 
level of 0 - 30 deg C. The pressure shall be maintained at 0.0133 pascals 
(0.0001 Torr) or less. All orbital operational conditions and all equipment 
functional modes shall be tested. The test shall include six complete hot-cold 
cycles ate the maximum predicted orbital rate of temperature change and with 
an 8-hour soak at each temperature extreme. Operating time shall be divided 
equally between redundant circuits. 

1. Test Model 

a. Thermal Balance Test 

Purpose: This test verifies the analytical thermal model and 

demonstrates the ability of the space vehicle thermal control subsystem to 
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maintain components, subsystems, and the entire space vehicle within specified 
operating temperature limits. This test also verifies the adequacy of component 
thermal design criteria. 

Test Description: The mass model shall be tested to simulate 
the thermal environment seen by the space vehicle during the transfer orbit and 
orbital mission phases. Tests shall be conducted over the full range of seasons, 
equipments duty cycles, solar angles, and eclipse combinations so as to include 
the worst case high and low temperature extremes for all space vehicle 
components. Special emphasis shall be placed on defining the maximum and 
minimum battery temperatures. The power requirements for all 
thermostatically controlled heaters shall be verified during the test. The test 
chamber pressure, with the test item installed, shall provide a pressure of 
0.0133 pascals (0.0001 Torr). The space vehicle thermal environment will be 
supplied by the Incident Flux method, where the intensity, spectral content, and 
angular distribution of the incident solar, albedo, and planetary irradiation is 
simulated. 

Test Levels and Duration; For this space vehicle two modes 
of operation will exist. The first, Earth sensing, shall be used during sensing 
period which will be two consecutive orbits per day. The second, Sun sensing, 
shall be used during all other periods. The worst case temperature extremes 
for each case shall be predicted and used for the test. The worst case high 
temperature will be predicted with all heat producing components operating at 
full power and the space vehicle receiving maximum solar irradiation. The 
worst case low temperature will be predicted at maximum eclipse and 
minimum heat production by onboard equipment. The test duration will be 
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dependent on the thermal time constant of the subsystems and orbital 
maneuvering. 

All tests are to be conducted at the Naval Research Laboratory 

(NRL). 

D. COMPONENT AND SUBSYSTEM LEVEL TESTING 

A component is a functional unit that is viewed as an entity for purposes of 
analysis, manufacturing, maintenance, or record keeping. Examples are 
hydraulic actuators, valves, batteries, electrical harnesses, and individual 
electronics packages such as transmitters, receivers, or multiplexers. A 
subsystem is an assembly of two or more components including the supporting 
structure to which they are mounted, and any interconnecting cables or tubing. 
A subsystem is composed of functionally related components that perform one 
or more prescribed functions. Typical space vehicle subsystems are electrical 
power, attitude control, thermal control, telemetry, structure, and propulsion. 

Component level random vibration testing is done as required to accept 
critical hardware. All mission specific electronics packages and mechanical 
devices are to be tested to appropriate random vibration levels. Levels are to 
be determined from the test model acoustic and random vibe system level tests 
and any analytical models available. Electronics packages and mechanical 
devices should alsoAll other components will be off the shelf, flight tested and 
validated quality and will be tested on a subsystem and space vehicle level only. 
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XI. COST ANALYSIS 


A. INTRODUCTION 

Program costs can be estimated using one of three standard methods — 
engineering build-up, analogy, paramedics — or a combination of the three. 
Cost modeling for this project was performed using the Unmanned Space 
Vehicle Cost Model, Revision 6 (USCM6). 

B. PARAMETRICS 

The USCM6 is a parametric estimating tool. Parametric estimates, 
also known as "top-down" estimates, or Cost Estimating Relationship (CER) 
estimates, are based on mathematical expressions relating cost as the 
dependent variable to selected, independent cost driver variables. The 
mathematical expression is derived by statistically correlating the historical 
cost data of several related or similar systems to physical and/or 
performance characteristics of those same systems. The implicit assumption 
of this approach is that the forces that affect cost in the past will continue to 
affect cost in the future. 

1. Advantages 

Parametric estimating has four advantages. 

a. Provides unbiased, objective, and consistent estimates. 

b. Once the parametric relatinships have been defined, the 
estimating process can be performed numerous times varying the cost 
sensitive physical or performance characteristic. 
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c. Since paramedics are derived from actual program cost 
history, they inherently account for the impact of schedule changes, 
engineering changes, system growth, and redirection. 

d. Since most of the parametric relationships are based on 
the advance design characteristics, and most cost estimates are reguired 
during the concept formulation phase of the program, the parametric 
approach provides a feasible means for performing an estimate. 

2. Disadvantages 

The main disadvantage of parametric estimates stem from several 

factors: 

a. The underlying assumption that what affects cost in the 

past will affect cost in the future. 

b. The need for a sufficient data base that is representative 

of the particular products developement, production, and technology 
environment. 

c. The parametric relationships may become obsolete, and 
require continual re-evaluation of the CERs and expansion of the data base. 

d. The estimate lacks program peculiarity. 

e. There is the limitation that extrapolation beyond the 
range of historical data becomes risky. 

C. USCM6 SCOPE 

The USCM6 provides CERs for estimating hardware costs for earth 
orbiting, unmanned space vehicle programs. These CERs are limited within 
the following boundary conditions: 
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a. The model addresses only unmanned earth-orbiting space 
vehicles, with the following exceptions — the stucture and thermal CER data 
bases include data points for large structures such as the Apollo Adapter, 
Apollo Module, Lunar Module Descent, Shuttle , Spacelab, ERBS, STOTA, 
UARS, COBE, GRO, and Space Telescope. 

b. The model is an approximation of the real world based upon 
mathematical relationships derived from the analysis of historical cost data. 
Implicit in this approach is the assumption that historical costs will properly 
reflect current and future costs. 

c. All costs included in the model are end-of-program costs or 
estimates of mature programs (at least one launch). 

d. The model's emphasis is on space vehicle hardware costs. It 
does not address specific missin payload costs other than communications 
payloads. 

5) Launch vehicle, launch adapters, and their associated ground 
support are not addressed in the model. 

e. CERs are based on burdened cost (direct plus indirect) with 
General and Administrative (G&A) costs included. 

£ A 95% cumulative learning curve is used to derive data base 
first unit costs. 

g. All CERs with the exception of Program Level CERs represent 
first unit cost and require appropriate consideration for production learning. 
Program Level recurring CERs estimate total recurring costs. 
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h. The cost model provides for an estimate calculated in U.S. 
Government Fiscal Year 1986 dollars, and must have the appropriate inflation 
costs added if costs are to be expressed as future expenditures. 

i. Target and incentive fee costs are not included in the model. 

j. The model yields an "average" starting point estimate which 
represents and "average" program with "average" problems, "average" 
technology, "average" schedule, "average" engineering changes, etc. (Ref. 
USCM6) 

D. ANALYSIS 

The USCM6 has provisions for subsystem and component cost 
estimating. Due to the preliminary design nature of this project and the 
limited component information the subsystem cost estimating approach was 
determined sufficient. Since the USCM6 provides for cost estimates in FY 
1986 dollars, a constant 7 percent annual inflation rate was applied to convert 
to FY 1991 dollars. The complexity factor was also set to one (1) due to the 
simple design concept used for this project. Most cost drivers associated 
with the cost model were dependent on subsystem mass and/or power 
requirements. A summary of mass and power requirements can be found in 
Chapter II. 

1. Non-recurring Costs 

The non-recurring costs associated with the subsystems are listed 
in Table 11.1. Appendix L shows the actual non-recurring cost estimate with 
inflation and complexity taken into account. The High Estimate and Low 
Estimate are calculated by taking the Most Likely Estimate and either adding 
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or subtracting, as appropriate, a correction factor derived from the Standard 
Error (SE) as a percentage of the calculated subsystem cost. 


Table 11.1 Non-recurring Costs 


Cost Estimating 
Relationship 
(CER)Category 

Most Likely Cost 
($millions) 

Low Estimate 
($million) 

High Estimate 
($millions) 

Integration and 
Assembly (I&A) 

3.34 

2.81 

3.88 

Structure 

15.27 

12.60 

17.94 

Attitude Control 
System (ACS) 

10.41 

5.31 

15.50 

Thermal Control 

1.64 

0.86 

2.42 

Electrical Power 
Supply (EPS) 

4.20 

2.29 

6.11 

Communications/ 

TT&C 

20.97 

11.39 

30.55 

TOTAL 

55.83 

35.26 

76.40 


2. Recurring Cost 

The recurring costs associated with the subsystems are listed in 
Table 11.2. Appendix L shows the actual recurring cost estimate with 
inflation and complexity taken into account. The High and Low estimates 
are derived in the same manner as for non-recurring costs. 
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Table 11.2 Recurring Costs 


Cost Estimating 
Relationship 
(CER) Category 

Most Likely Cost 
($millions) 

Low Estimate 
($millions) 

High Estimate 
($millions) 

Integration and 
Assembly (I&A) 

1.20 

1.05 

1.36 

Structure 

2.60 

2.29 

2.92 

Attitude Control 
System (ACS) 

4.19 

2.05 

6.32 

Thermal Control 

0.49 

0.32 

0.66 

Electrical Power 
Supply (EPS) 

0.14 

0.08 

0.19 

Communications/ 

TT&C 

8.98 

6.92 

11.03 

TOTAL 

17.60 

12.71 

22.48 


3. Launch Vehicle Costs 

Table 11. 3. below is a cost comparison between the Titan 2 or 
Delta II launch vehicles. The Titan 2 figures are based on the Operational 
development Model (ODM) Freeflyer presentation given at the Naval 
Postgraduate School by Charlie Merk and P. Regeon of the Naval Research 
Laboratory on October 22, 1991. They were verified to be reasonable 
figures by Martin Marietta. The Delta II figures were provided by 
McDonnell Douglas. 
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Table 11.3 Launch Vehicle Cost Comparison ($millions) 


Titan nr.ni Pella 11(73201 

(2) 

Booster 

20.6 

Booster 

21.3 

Extended ACS 

1.0 

Failure Liability 

1.6 

Integration 

8.4 

Integration 

9.3 

Launch Services 

9.6 

Launch 

Operation 

3.9 



Range Support 

1.9 

Mgmt Reserve 
(5%) 

_LQ 

Mgmt Reserve 

Li 

TOTAL 

$40.6 


$39.1 


Notes: (1) Costs provided by Martin Marietta 
(2) Costs provided by McDonnell-Douglas 


E. CONCLUSIONS 

As with all cost estimating schemes the use of parametrics is 
not without errors. The payload, and on orbit operation support costs are 
not included with the overall cost figure. Launch cost, as previously noted, 
are approximaltely $40 million. Payload and on orbit operation support costs 
were unavailable for incorporation in this cost analysis. The spacecraft bus, 
design and production, and launch costs represent an overall most likely cost 
of approximatley $160 million. Program level, aerospace ground equipment 
and launch & orbital support are included in the total program costs along 
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with the non-recurring and recurring hardware costs. This cost is within 
realistic limits for a space vehicle of similar mass and complexity. A cost 
summary is provided in Table 1 1 .4. 

Table 11.4 Cost Summary 


Category 

Most LikelyCost 
($millions) 

Low Estimate 
($millions) 

High Estimate 
($millions) 

Satellite 

73.43 

47.97 

98.88 

Launch 

39.1 

39.1 

39.1 

Program Level 

29.59 

22.68 

36.50 

Aerospace 

Ground 

Equipment 

(AGE) 

9.64 

7.14 

12.15 

Launch & Orbital 
Ops Support 
(LOOS) 

3.53 

2.71 

4.36 

TOTAL 

155.29 

119.60 

190.99 
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APPENDIX A - STATEMENT OF WORK 


A. SCOPE 

This document establishes the preliminary requirements for an infrared (IR) 
surveillance spacecraft. 

B. REQUIREMENTS 

1. Mission 

This spacecraft will be used to sense areas of the earth subtended by a 4 
by 4 degree field of view (FOV) for footprint toe grazing angle of 90 to 45 degrees. 
The spacecraft orbit shall be circular, with a height of 275 nautical miles and an 
inclination of 70 degrees. 

2. Threat 
Classified 

3. System Functions 

a. Payload Subsystem 

The payload is an IR sensor which consist of the following: 

• Telescope 

• Dewar /IR focal plane array 

• Cryogenic cooler (65 degrees K) 

• Data Processor 

• Pointing controller and sensors 

• Structure/Optical Bench 

b. Physical Characteristics 

(1) Weight 800 pounds 

(2) Volume 

Width 60 Inches 
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Length 60 Inches 
Height 20 Inches 

Note: Use one 60 by 60 inch side for mounting the spacecraft buys 
and leave the other 60 by 60 inch side clear for payload thermal radiators and the 
telescope aperture. The boresight of the telescope FOV is at a 45 degree angle 
with respect to the plane of this side. 

(3) Power 

Standby Mode - 100 watts 
Sensing Mode - 150 watts 

c. Attitude Control Subsystem 

A three-axis stabilization capability shall be provide to provide nadir 
and offset-pointing for payload operation (sensing mode). In the nadir pointing 
mode, the spacecraft body frame tracks the rotating local vertical reference system 
(pitch, roll, and yaw = zero). The sensor boresight may then be offset, as require, 
about the roll axis to achieve grazing angles of between 90 degrees and 45 
degrees, for the toe of the footprint of the telescope FOV. 

Output form an inertial measurement unit (IMU), which is part of the 
payload subsystem, is available for short term attitude reference, The IMU is 
initialized by a pair of star sensor, which are also part of the payload. 

• Attitude Accuracy ±5 Degrees Per Axis (3a) 

• Time To Slew 90 Degrees in 15 minutes 

• Slew Settling Time 1 Minute 

• Rate Stability 0.003 Degrees/Sec (3a) Per Axis 

d. Command and Telemetry Subsystem (CATS) 

The CATS shall: 
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a. Decrypt, validate, decode, store and execute commands to the 
spacecraft subsystems in accordance with commands and data received from the 
ground sited via the RFCS (see below); and 

b. Process, format and encrypt spacecraft subsystem 
housekeeping data for transmission to the ground via the RFCS. 

e. Radio Frequency Communications Subsystem (RFCS) 

The RFCS shall provide redundant: 

a. Omnidirectional command receiving capability; 

b. Omnidirectional low data rate (16 Kbps minimum) transmitting 
capability; and 

c. Directional high data rate (150 Kbps maximum ) transmitting 
capability. 

For the high data rate downlink, design, for: 

a. Optimal spectrum usage; 

b. BER 10' ^ encrypted (10' 9 unencrypted); 

c. Minimum Elevation angle of 10 degrees at data collection 

terminal; 

d. 99 Percent link availability due to rain attenuation; 

e. Assume a 20 foot maximum diameter ground antenna. 

f. Electrical Power Subsystem (EPS) 

The electrical power system shall be designed to provide 30 minutes of 
sensing time per day. Initial sizing should assume 15 minutes of sensing time per 
orbit, for two contiguous orbits. This will allow some warm up time and as much 
data collection as the pass geometry allows. The solar array and battery shall be 
sized to proved energy balance for day and/or night sensing during the maximum 
orbital eclipse season. 
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APPENDIX C - ORBITAL ANALYSIS 


I. ESTIMATION OF ALTITUDE LOSS DUE TO DRAG USING 
MANUAL CALCULATIONS 

A back of the envelope prediction of the orbital decay encountered by drag forces is 

computed below. A simple atmospheric model based on an exponential atmosphere with a 

reference altitude and scale height is used and shown 

(3 = .0476 km 
r 0 = 140 km 

p 0 = 3.0 x 10' 6 kg/m 3 

P = Poe^<""> (C.l) 

Using the computed density of 7.0 x 10*^ kg/m A 3 for 509.3 km orbital altitude, the decay 
of the orbit can be estimated using Equation C.2. 


Ar = -2 


C d A 


m J 


f 


a (1 + e) I p dv per orbit 


(C.2) 


Ar = 1 .657 km/year 

Using a coefficient of drag equal to 2.4 and an effective area of 3.0 m A 2 the orbit will 
decay 1.657 km per year. Over a three year lifetime this will exceed the allowable altitude 
perturbation. 


II. ESTIMATION OF ALTITUDE LOSS DUE TO DRAG USING 
ASAP PROGRAM 

The ASAP program also computed the orbital decay caused by the atmospheric drag. 
The atmospheric model used within the program is the 1976 Standard Atmospheric model. 
The program was executed with all of the internal models, except drag, turned off to ensure 
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all the effects seen were due only to the atmosphere. The estimated altitude decay is 
approximated to be 4.5 km per year, indicating a need to reboost the satellite at six month 
intervals to maintain the orbit within defined specifications. Selected outputs from the 
ASAP program are include below and are summarized in Table C.2 

Table C.2. ASAP Orbital Perturbations from Drag 


Geocentric Radius (km) 

Change in Altitude (km) 

Time from epoch (days) 

6889.00 

0.00 

0 

6888.15 

0.85 

90 

6887.38 

1.62 

180 


A coefficient of drag on the satellite was estimated to be 2.4, slightly more than the 
recommended 2.0 to 2.2, this was to ensure that the perturbation was of a maximum value. 
The effective area used by the program was 3.0 m A 2, the average value of the minimum 
planform area of 2.25 m A 2 and the maximum value of 3.75 m A 2. An exponential 
extrapolation of the data yields a decay in the orbit of approximately 4.5 km per year. This 
decay clearly exceeds the allowable 1.852 km for the lifetime of the satellite. Because the 
value for altitude decay from the ASAP program is larger than the hand calculation value, 
the ASAP value will be used in on orbit fuel consumption calculations. 
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SELECTED ASAP OUTPUT 
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ORiGii^AL F,m {« 
OF POOR QUAUTY 



III. EFFECTS ON INCLINATION FROM EARTH’S 
OBLATENESS 

As mentioned earlier there is no cumulative effect on the inclination of the 
orbit. The Earth's equatorial bulge does exert an influence on the satellite, but this is small 
in nature and has a period equal to the orbital period of the satellite. This effect can be 
shown by relating two equations which evaluate the rotation of the longitude of the 
ascending node. Equation C.3 relates dQ/dt with the force the Earth's bulge exerts and 
Equation C.4 relates diT/dt with other orbital parameters. 

dQ _ r Fj 2 sin ( u ) 

n a 2 V(1 -e 2 ) sin ( i ) 


- 3 * ^[ffcosU) (C.4) 

dt 

These two equations are equated and an expression for the J2 force is found this expression 
is then substitute into Equation C.5 and an expression for di/dt can be found based only on 
orbital parameters (Equation C.6). 



di _ r Fj 2 cos ( u ) 
n a 2 V(1 - e 2 ) 


(C.5) 


di - ~ 3 tt J 2 cos( i ) sin( i ) 
dt tan( u)f \ I 

It now can easily be shown that di/dt varies with the parameter u (the argument of latitude). 
The effect is periodic with a period equal to the orbital period and singularities at u=0° and 
180°. 


R 


, 1.5 


(C.6) 
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IV. COMMENTS ON PROGRAMS USED 

1. ASAP 

This program was useful in the drag estimations. The input file was modified for 
our particular satellite and the requested data displayed on the screen in a scrolling fashion. 
The display of the data was very hard to cope with and modification should be made to 
have selected data dump into an output file for retrieval and manipulation. 

2. Orbital Workbench (Version 1.01 Cygnus Engineering) 

This program was useful in modeling the launch vehicle sequence (discussed 
later) and general orbit and ground track manipulation. The input sequence to execute a 
problem requires the entry and exit of about 8 submenus making it hard to keep track of 
what has been changed. The output data is displayed in a graphical format, which is very 
useful, and into a data file for further manipulation. The drag of the satellite was not 
modeled with this program because the time overhead associated with the graphical display 
slowed down the computation enough to make it impractical. 
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This MathCAD document calculates the effects of the sun and the moon 
on the inclination of a satellite. The first portion sets up necessary 
constants and then is followed by calculations for the sun and then the 
moon. 
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This MathCAD document calculates the precession of the longitude of the 
ascending node for a circular orbit. It also displays the degrees the ground 
track will move at the equator - based on the precession of the orbit and the 
rotation of the earth. 


For circular orbits the precession of the orbital plane can be expressed as: 
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APPENDIX D - LINK MARGIN CALCULATIONS 


Specification 

Value 

(Units) 

Omni-directional Low Data Rate 

16 

kbps 

Directional High Data Rate 

150 

Mbps 

Bit Error Rate (BER): 



Encrypted 

IE-11 


Unencrypted 

IE-09 


Data Collection Terminal (DCT): 



minimum elevation angle 

10 

deg 

Satellite Orbit Altitude 

275 

nmi 


509643.39 

m 

link distance (R) 

1718156.3 

m 

Assume receiver antenna diameter 

20 

ft 


6.096 

m 

Assume antenna effeciency of 

0.45 


Effective receiver antenna area 

13.133858 

m A 2 

99% Avail, due to rain attenuation 

-10 

dB 

Omni-directional Link Calculations: 



Down-Link Telemetry Frequency 

2247.5 

MHz 

Corresponding wavelength 

0.1334816 

m 

Signal: 



Power of Transmitter (Pt) 

43.537235 

dBm 


22.579976 

Watts 

Turnstile Antenna Gain (Gt) 

1 

dB 

Gain of Receiver (Gr) 

39.667591 

dB 

Directional Link Calculations: 



Down-Link transmitter freq 

20 

Ghz 

corresponding wavelength 

0.015 

m 

Signal: 



Power of Transmitter (Pt) 

48.757616 

dBm 


75.121038 

Watts 

Gain of Transmitter (Gt) 

37.804199 

dB 

antenna diameter (m): 0.5 



Gain of Receiver (Gr) 

37.362723 

dB 
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APPENDIX D 

Frequency Communication Subsystem (RFCS) 


Losses: 

Loss due to medium (rain) 

-10 

dB 

Signal/Path Loss (Ls) 

-183.16363 

dB 

Incidental Loss (Li) 

-1 

dB 

Transmission Line (Lt) 

-3 

dB 

Actual Transmitted Power (Ca) 

Low Data Rate 

-112.9588 

dBm 

High Data Rate 

-73.239087 

dBm 

Receiver Noise Density: 

k 

-198.6 

dBm/Hz-K 

250W Traveling Wave Tube Amp. (Te) 

1000 

K 

10*log(Te) 

30 

dBK 

Noise Power Spectral Density (No) 

-168.6 

dBm/Hz 

Data Rate (B): PCM 

Omni-directional 1 E-09 

42.0412 

dB-Hz 

Directional 1 E-1 1 

81.760913 

dB-Hz 

Eb/No Available 

16.6 

dB 

Eb/No = (Ca - B) - No 

13.6 

dB 

Margin 

3 

dB 
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APPENDIX E - ELECTRICAL POWER SUBSYSTEM (EPS) 


A. RADIATION DEGRADATION 

From the “Solar Array Radiation Degradation Handbook”, for a 275 NM, 70 
inclination, circular orbit, the annual 1 MeV equivalent electron fluences are as 
shown in Table 4. 

Table E.l. Annual 1 MEV Equivalent Electron Fluence 



isc 

Voc & Pmax 

Electrons 

3.95 x 10 11 

3.95 x 10 11 

Protons 

2.11 x 10 12 

4.08 x 10* 2 

TOTAL 

2.51 x 10 12 

4.48 x 10^ 2 


After three years (spacecraft design life) the anticipated fluences are as shown 
in table E.2. 

Table E.2 1 MEV Equivalent Electron Fluence At EOL 


Isc 

Voc & Pmax 

7.53 x 10 12 

1.34 x 10 13 


The cells selected are 2.5 x 6.2 cm Si, back surface reflector, 0.02 cm thick 
(8 mil), 10 O | x, with a 3 mil shield and Ti Ox AI2O3 anti-reflective coating. 

The effects of radiation on the solar cells are summarized in table E.3 below. 
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Table E.3 Radiation Effect On Solar Cells 



EOL 

BOL 

Pmax (mW/cm 2 ) 

18.7 

20 

Voc (mV) 

587 

608 

ISC (mA/cm 2 ) 

43.2 

43.7 

lMp(mA/cm 2 ) 

40 

40.2 

Vmp (mV) 

469 

502 


B. BATTERY DESIGN 

The payload voltage specified is 28 (+ 4) VDC. Since 24V is the minimum 
allowable bus voltage, Vdb is 24V. Assuming an open circuit failure of one battery 
cell * 

v DB = (N-i)v D -v DD 

Vd = 1.2V (Ref: AN IMPROVED, HIGH SPECIFIC ENERGY 
( LIGHTWEIGHT j SPACE BATTERY SYSTEM , Mr. Lee Miller, Eagle-Picher 
Industries) 

Vdd — 1.1V (Bypass diode voltage drop in event of open circuit failure of one 
cell) 

Therefore, N = 21.9 or 22 

Now, recomputing for Vdb yields Vdb = 24.1V. 

To determine the battery capacity necessary to accomplish this objective, 
first the initial "energy deficit" (E 0 ut) must be determined by computing the energy 
that is drawn by the payload and bus in the earth oriented mode with the two 
active events mentioned earlier and subtracting the incidental solar input during this 
period. 
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Eou-SPbth+SPA 

Table E.4 Housekeeping Power Budget (Watts) 


( 1 ) 


LOAD 

Ph 

Payload (Sensor Idle) 

100 

Attitude Determination and Control 

50 

Thermal Control 

80 

Propulsion 

1 

TT & C 

5 

TOTAL 

256 


Table E.5 Active Power Budget (Watts) 


LOAD 

Pa 

Payload (Sensor Active) 

50 

Datalink 

80 

Margin 

20 

TOTAL 

150 


Table E.6 Time Intervals During Earth Seeking Event 


EVENT 

TIME (MINUTES) 

Passive (th) 

194.6 

Active (t a ) 

30 


Evaluating Equation 1 and using data contained in Tables E.1-E.3 yields an 
energy deficit of approximately 2 X 10^ Joules. The capacity of the storage device 

may now be evaluated in the following manner, 

p E oirr _ 2 x 10 _ oo u 

V db EX>D 24.1x0.6 
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C. EPS SIMULATION 


A simulator was developed to evaluate the demands placed on the energy 
storage system and to visualize the nature of the DOD cycle for a nominal active 
event under worst case conditions (EOL, winter solstice, maximum possible 
eclipse, etc.) 


10 REN M«tMM*M»t**MmM*MM*m**«**l«»«**«M»* 

20 REM SIMULATION OF SATTERY OOO 

SO REM aaaa*aa*«*aa*aaaaaaa»»»»**»*a»**aaa*aaaaaaaaaaaaaa 

40 REM SIMULATION RUNS FOR A 24 HOUR DAY (1440 MINS), MAX ECLIPSE SEASON, AT EO 
SATTERY DOO STARTS AT 2ERO ANO PROCEEDS ACCORDING TO PARAMETERS THAT HAVE BEEN 
INPUT 

60 REM MtttMttMttMVMtattMMtamaMMtMtltaMttMMMMVfniMtMMVVVM 
SO REM aaa»»aaaaaa*#*Ma**a*a*«*-*aaa*aa 4 *aaa*aaa*aaaaaa 4 *aaaaaaa»aaaaaaaaaaaaa 

SO REM VDB = MINIMUM BATTERY DISCHARGE VOLTAGE 
70 REM V6C = SATTERY CHARGING VOLTAGE 

60 REM ERATE IS ANGULAR VELOCITY OF SPACECRAFT AROUND EARTH 
•0 REN C s SATTERY CAPACITY (AH) 

100 REM PH s HOUSEKEEPING POWER 
110 REM PA = ACTIVE POWER 

120 REM UNITPWR ■ MAX POWER AT EOL («*/ CH A 2) 

130 REM aaaaaaataaaaaataaaaaaaataaa*naaaaaaaaaaaaaaaaaaaa 
140 REM OUTPUT GOES TO A FILE CALLED DCALC.PRN IN ASCII 
ISO REM aa*aaaaaaaaaaaaaaaaaaaaa******aaaaaaaaa*aaa*aaaasa 

ISO OIM 000(1440) 

170 VDB*24 . 1 
ISO VBC-34.2 
190 PI«3. 141662664* 

200 ERATE*. 066466733* 

210 C=40 
220 TC*0 
230 PC*C/10*VBC 
240 PH-160 
260 PA* 230 
200 PC=PC/.85 

270 UNITPWR* IS. 7/1000*100*2 

2S0 AREA* PC /UNITPWR 

260 KaVDfltCMO 

300 000( 0 )*0 

310 FOR T*1 TO 226 

320 IF T<6 THEN SFR* 1 : PO*0 : GOTO 400 

330 IF T <47 THEN SFR-0 : PO-PH : GOTO 400 

340 IF T< 62 THEN SFR*COS( ERATE*( T-47 )-PI/2 ) : PO»PA*PH*SOR( 1 -SFR*2 > : GOTO 400 
360 IF T< 65 THEN SFR-COS C ERATE*( T-47 )-PI /2 ) : PO*PH»8GflM -BFR~2 ) :GOTO 400 
360 IF T< 142 THEN 6FR*0 : PO*PH : GOTO 400 

370 IF T< 1 67 THEN 6FR-COS( ERATE* ( T-1 42 )-PI/2 ): P0*PA+PH»80fl( 1 -SFR*2 ) :GOTO 400 
360 IF T< 166 THEN SFR-COS (ERATE* (T-1 42) -PI /Z ) :PO*PHaSQR( 1-SFR*2) : GOTO 400 
360 8FR-0: PO*PH 

400 DOO( T )*DOO( T- 1 )- < 8FR*UNITPWR*AREA-PO )/K 
410 NEXT T 

420 FOR T =230 TO 1440 

430 IF T=230+TC*85 THEM SFR* 1 : PO=0 : GOTO 460 
440 IF T=288*TC«66 THEN SFR*0 : PO»PH : TC-TC+ 1 
460 DOO(T )*DOOf T-1 )-( SFR*UNITPWR*AREA-PO)/K 
460 NEXT T 

470 OPEN "C:\DCALC\OCALC.PRN" FOR OUTPUT AS il 
460 FOR I«0 TO 1440 
480 PRINT * 1 ,DOD(I) 

600 NEXT I 
610 END 


Figure E.l Battery Simulator 
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APPENDIX F - LAUNCH VEHICLE SELECTION 

A survey of available launch vehicles was conducted (see table F.l). The 
preliminary launch vehicle selection was based on a spacecraft mass of 800 to 
1400 kilograms, a circular orbit of 509.3 kilometers, and a 70 degree inclination. 
Based on the sensitive nature of the payload, choices were limited to American 
vehicles. Launch vehicles such as Scout and Pegasus were found to be too small 
while Atlas Centaur, Titan 3, and Titan 4 were too big. The closest matches were 
Adas E, Titan 2, and Delta II. The Atlas E was immediately eliminated since the 
program is soon to be cancelled and has very few launch vehicles available. 
Additionally, the Delta II (6920) will not be available when the spacecraft is 
launched. The best choices were the Titan 2 and the Delta II (7320). 

Additionally, the option of using the space shuttle was considered. The shuttle 
normally releases satellites between 200 to 400 kilometers at about a 28.5 degree 
inclination. A kick motor such as the PAM series (developed by McDonnell 
Douglas for the Delta launch vehicles) would be required to adjust the altitude and 
inclination. This option was decided against due to the extra costs associated with 
the human safety features required by space shuttle payloads. 

Since both Titan 2 and Delta II can place the spacecraft at the correct 
altitude, they were compared in four areas: orbital insertion, cost, limit loads, and 
performance records. A two stage Delta II rocket would place the satellite in a 
circular orbit very close to that desired. The Titan 2 requires an apogee kick motor 
to circularize its orbit, adding a measure of complexity. Both launch vehicles cost 
virtually the same (see cost section Chapter 11). The limit loads of both launch 
vehicles are similar on the lateral axis but vary greatly on the vertical axis - Delta 
II is 7.2 g’s while Titan II is 12.5 g's (based on a 1000 kg payload). The Delta 


161 



program has had over two hundred launches since 1960 and has a 98% success 
rate over the last twelve years. Titan 2 has had only three launches since the 
Gemini program. These factors lead to the choice of the Delta II (7320). 


Table F.l Launch Vehicle Survey* 


Launch 

Vehicle 

Launch 

Sites 

Satellite 
Mass (kg) 

Insertion 
Altitude (km) 

Notes 

Atlas Centaur 



185 

1 

Atlas E 

E/WSMC 

2350 

Atlas 1 

ESMC 

5900 

185 

2 

Atlas 2 

ESMC 

6780 

185 


Scout 



550 


Scout 1 

WSMC 

1 65-220 

3 

Scout 2 

WSMC 

353-451 

555 


Titan 



185 polar 


Titan 2 

WSMC 

2177 

4 

Titan 2 

WSMC 

3400 

648 km & 70 deg incl 

5 

Titan 3 

E/WSMC 

14742 

185 


Titan 4 

ESMC 

17690 

185 


Delta II 



185 

A 

Delta 6920 

E/WSMC 

3983 

6 

Delta 7320 

E/WSMC 

2100 

510 km & 70 deg incl 


Delta 7920 

E/WSMC 

6700 

648 km & 70 deg incl 


Pegasus 

NA 

408 

463 



Note 1 : Includes a Star 20 motor. Only two vehicles remain as of June 1991. 


Program scheduled to be terminated in 1993. 

Note 2: ESMC = Cape Canaveral 

Note 3: WSMC = Vandenberg Air Force Base 

Note 4: Only three Titan 2 launches since 1966. 

Note 5: Includes a Star 37 kick stage. 

Note 6: Replaced by the 7000 series in late 1992. 

* Only US launch vehicles were considered due to the sensitive nature of the payload. 
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APPENDIX G - PROPULSION SUBSYSTEM 


A. CALCULATION OF REQUIRED PROPELLANT 

Monopropellant hydrazine (N 2 H 4 ) with a specific impulse, I s , of 225 seconds 
was used as the fuel. An estimated mass of 1050 kg for the satellite was used 

throughout the following equations. 

1. Spin Down and Detumble 

A worst case of 5 rpm was designed for to include the possibility of 
tumbling upon separation (or later on). The following inputs and calculations 
determined the propellant expended by this maneuver: 


Approximate Spacecraft Spin Moment of Inertia — Im— 475 kg-m 
Thruster Moment Arm = r = 0.5 m 


= Im* 

Is * go * r 


475 * (5 * (^)) 
oU 

225 * 9.806 * 0.5 


0.225 kg 


2. Orbital Insertion 

The worst case of an 18.5 km insertion error was designed for. The 
orbital section (Section III) shows that the delta v required is 10.238 m/s (includes 
apogee and perigee firings). To do this maneuver four 1 lb f thrusters were used 
canted at 35.8 degrees from the thrust axis. The following equations determine the 
fuel required and the time to station. 

F = m*a = m*^p- 
at 

T * cos 0 =4rL*v 2 *cos0 = I s *g o *(^ L )*cos0 
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Where F is the axial force, T is the thrust produced, and 0 is the cant angle 
of the thrusters. When the cant angle is zero, F equals T. Otherwise F equals T * 
cos 0. Plugging into the equations: 

T = 4 * (1 lb f ) * (4.448 N/lb f ) = 17.792 N 

T = m * (Av.) = 1050 * ( 10,238) _ 17 792 

cos 0 At cos 35.8 At 

At = 745 sec 

T = I s * go * (^J 31 ) = 17.792 = 225 * 9.806 * (^-) 

Am = 6.01 kg 

3. Atmospheric Drag Corrections 

Drag calculations are included in the orbital section (Section III). The 
result is a delta v of 8.172 m/s. The equations above apply. Using compact forms 
of these equations: 

Am = — mJLte =4.80 kg ^ = Am *I S * g 0 = 595 sec 

Is * go * cos 0 T 


B. PROPELLANT/PRESSURANT TANK REQUIRED DIAMETER 

To decide on a fuel tank from those available commercially, the required tank 
diameter was determined below. The tank volume was assumed to be 90% fuel 
and 10% pressurant. 

Density = p = 876.2 kg/m 3 

Tank Volume = V = 21 = JUL22. = 0.0157 m 3 
p 876.2 

~~V~ 


Tank Diameter = D = 2 * 3 


0.9 * (|) * 7t 


= 0.322 m 
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C. DEORBIT CALCULATIONS USING THE ASAP PROGRAM 


RUN STARTS ON JULIAN DATE -» .2..44992430000000QOE+06 

RUN ENDS ON JULIAN DATE 3.4572295O00OOQ0O0E+06 

REFERENCE JULIAN OATE OF PN AND. EPHEM * ' 2.4499245000000000B+06 

O DAYS O. OOOOOO HOURS PROM EPOCH " 

' CARTES IAN COORD X, Y, Z, XD, ID * ZD 
6.81293B265O4205E+03 O. OOOOOOOOOOOOOOE+OO '0;00000000600000E+00 0. OOOOOOf 

Pi MEAN LONG 

6.8564737Q178488E+03 0. OOOOOOOOOOOOOOE+OO ^.349537478936888-03 O. OOOOOOC 

CLASSICAL ELEMENTS E r I, NODE, W* MA, TA, BA, TRUE LONG, ECC LONG, ARG C 
6.349S3747893688E-03 7 .001307171693088+01 O.OOOOOOOOOOOOOOB+OO 0.0000001 

-O.OOOOOOOOOOOOOOE'fOOO.OOOOOOOOOOOOOOS+OO-OiOOOOOOGOOOOOOOB+OO O.OOOOOOt 

OTHER PARAMETERS R>- V , ENODE , ' LONG , LAT, HOUR ANGLE/ ALT, ELLIPSOIDAL A3 
«. 812938265042058+03 7.67319245670377E+00 -.■•3- 15000000000000E+02 3 ; 15Q000< 

4 .34798265042048E+02 4 . 34?98265042048B+02. . 4 .34798265042048E+02 1 . 569494: 

365 DAYS 0.000000 HOURS PROM EPOCH 
CARTESIAN COORD X, Y, Z, XD, YD, ZD .. 

-2.67245331340190E+03 6 .113073826233708+03 1.60102362248374E+03 -1.77138H 

EQUINOCTIAL ELEMENTS A, H, K, P, Q, -MEAN LONG 

6 .848 182 764829 83 E+O 3 -4 . 38277941839643B-03 -3 * 23545774199522E-03 6.637399C 

CLASSICAL ELH4ENTS E, t, NODE, K, MA, TA, BA, TRUE LONG, ECC LONG, ARG C 
5.44785474590269E-03 7 . 00103118090051E+01 I.08605572273694E+02 1.249589; 

2 . 494 18414128877E+02 2 .497109010506548+02 1, 22982909259365E+02 “ 1 . 232753S 

- OTHER PARAMETERS R, V, ENODE, LONG, LAT, HOUR ANGLE/ ALT, ELLIPSOIDAL A1 
6 . 86111907468798E+03 7 . 61484086803961B+00 6.38570964976474E+O1 6.886503C 

4.82979074687983E+02 4.8414907S656875E+02 4 .32736229490193E+02 1.566648: 

730 DAYS '• 0 . OOOOOO HOURS PROM EPOCH -V. 

CARTESIAN COORD X, Y, Z, XD, YD, ZD 

' ■ 2 74 6 38 112869 590E+03 -9.055811845367868+02 6.22696343905423E+03 5.741067: 

EQUINOCTIAL ELEMENTS A, H, R, P, Q, MEAN LONG 
6.82S56149289003E+03 5.46896376789193B-03 -3.02071378898686E-03 -3.8S4251C 

CLASSICAL ELEMENTS K, I, NODE, tf, MA, TA, EA, TRUE LONG, ECC LONG, ARG C 
6.24774171116973E-03 6.998«704?062779B+01 2 . 13406815154228B+02 2.655Q67C 

X.99640768542655E+02 X. 99761445986203K+02 3 . 18554285373984E+02 3.186749c 

OTHER PARAMETERS R, V, ENODE, LONG, LAT, HOUR ANGLE/ ALT, ELLIPSOIDAL A1 
6.8656944482814 0E+O3 7.597069599905288+00 X.68909863602135E+Q2 2.9725375 

4.87554448281402E+02 5. 05161338085Q73E+02 4 .04777147648743E+02 1.558892C 

. 1095 DAYS 0.000000 HOURS FROM EPOCH 

;CAKTESIAN COORD X, Y, Z, XD, YD, ZD 1 I.'.’.. 

4.67261772554053E+03 -4 . 93634260Q84055E+03 -1 . 61596798173121E+02 2.013502C 

EQUINOCTIAL ELEMENTS A, H, . K, P , Q, MEAN LONG 
6.83231582868964E+03 -2.24178026696187E-03 4 .75798163233606B-O3 -5.044273' 

CLASSICAL ELEMENTS E, I, NODE, M, MA, TA, BA, TRUE LONG, ECC LONG, ARG C 
5.25965473952302E-03 7.00D96508283337E+01 3 . 13923405235067E+02 2.084856: 

3.37702193789004E+02 3.37816257723682B+02 3 . 12474160777036E+02 3.125882; 

OTHER PARAMETERS R, V, ENODE, LONG, LAT, HOUR ANGLE/ ALT, ELLIPSOIDAL M 

6.79904024161415E+03 7.67538372182930E+00 2.69677977906928E+O2 2.691824' 

4. 2090024 1614150E+O2 4 .20912383026455E+02 4 ,182402063593 54E+02 1.561206' 

1460 DAYS 0. OOOOOO HOURS FROM EPOCH 

CARTESIAN COORD X, Y, Z, XD, YD, ZD 

3.3178543939844SE+03 5.38518767676673E+02 -5.95655875473560E+03 3.85B722; 

EQUINOCTIAL ELEMENTS A, H, K, P, Q* MEAN LONG . 

6.80702412514308E+03 -1. 064133938654911-03 -5. 396330756714Q9E-03 5.316834c 

CLASSICAL ELQtENTS E, I, NODE, V, MA, TA, EA, TRUE LONG, ECC LONG, ARG C 
5. 50025151018 166E-03 6.99856599846952E+01 ... 4.94223376122076E+01 1.417330: 



APPENDIX H - ATTITUDE CONTROL TOPICS 


A. COORDINATE SYSTEMS 

Four coordinate systems were used to perform the analyses in the different 
modes: (1) the spacecraft body, (2) the nadir pointing/sensing, (3) the sun tracking, 
and (4) the orbit normal. 

The body Z-axis is the outward normal from the payload with the origin at the 
center of mass. The body X-axis is the outward normal from the freestream face 
and the Y-axis is orthogonal to X and Z. 

The sensing coordinates are the standard satellite nadir pointing coordinates 
with the Z-axis nadir pointing, X in the direction of flight, and Y orthogonal to both. 
The origin is coincidental with the body coordinates. 

The Z-axis of the sun tracking coordinates tracks parallel with the incoming 
sun vector, X is orthogonal in the equatorial plane , and Y completes the right hand 
set. For the purposes of analysis, these coordinates are an "inertially" fixed frame 
because of the slow change in the earth-orbit-sun orientation. 

As a base for coordinate transformations between modes, the orbit normal 
system was used. The Z- axis in this system points out the normal for the orbital 
plane, X is orthogonal out the ascending node, and Y completes the right hand set. 
This also can be considered "inertially" fixed because of the slow rotation of the 
orbital plane. 

B. THE EQUATIONS OF MOTION 

The general equations for a three axis stabilized, three reaction wheel system 
in terms of euler angles are quite lengthy and involved. For this particular design, 
assumptions reduced these equations to two forms: (1) the pointing mode form 
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where the attitude errors remain in a small linear range, and (2) the slew mode 
form where the disturbance torques are small compared to the slewing torques. 

In the pointing mode equations the small products of inertia are neglected as 
well as the small stored angular momentum in the RWAs. The control torque 
produced by the RWA is governed by the standard position plus rate feedback law 
leading to an uncoupled second order system for each axis. These are described by 
the following: 


(H.l) I X x <t» + K$ T* <J) + Kft = I M x 
(H.2) i yy e + K e t 0 e + K e e = z M y 
(H.3) IzzV "*■ Ky'tyV KyY = ^ M z , 

where Mi are the external torques acting on the spacecraft, In are the principal 
moments of inertia, and Kjand Tj are the respective system gain and time 
constant.The derivation of these equations can be found in ref. (1). 

In the slew mode equations, the gravity gradient, aerodynamic, and solar 
disturbance torques are 3 - 4 orders of magnitude lower than the slewing torques 
and are neglected for the slew analysis. Assuming the products of inertia are 
negligable, the three axis / three reaction wheel equations reduce to: 

(H.4) I xx^x + a y Oz( Izz - Iyy ) + Gyh z - = *hx 

(H.5) IyyQy + Ox^zC IxX - I Z Z ) + ^X ’ ^xh Z = ‘hy 
(H.6) I ZZ Q Z + Q x O y( Iyy - Ixx ) + ^xhy - Qyh x = -h z , 
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where Qi are the angular velocities about each body axis described in terms of one 
of 12 euler transformations in terms of <l>, 0, and y, the roll, pitch and yaw errors 
respectively. The hj are the angular momenta produced about each axis by the 
RWAs. The various euler forms of the Q, , along with the general forms of 
equations G4 - G6 can be found in ref. (2). It is easy to see that even these 
simplified forms are highly coupled. 

C. DISTURBANCE TORQUE MODELS 

The following MATLAB code contains the disturbance torque models used in 
this design for both the sensing and sun tracking modes. The sense mode gravity 
gradient equations are derived in reference (1). The gravity gradient torque model 
for the sun track mode was developed by this designer. The aerodynamic moments 
model was developed by this designer assuming the cubic geometric center to be 
the center of pressure. Solar torque equations are derived in Ref (1) for 
geosynchronous spacecraft. Using suitable coordinate transformations these 
equations are applied to this design. Finally , the magnetic torques available for 
momentum dumping are included. This model was also developed by the designer. 
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% PROG DISTURBANCE TORQUES 

% THIS PROGRAM CALCULATES THE DISTURBANCE TORQUES FOR 
HTSI IN BOTH SENSING AND % SUNTRACKING MODES. ( NOTE IF 
USED AS A SUBROUTINE, THE SUBROUTINE FLAG MUST % BE SET 
AND A TIME VECTOR DEFINED ) 


sflag=input(' subroutine use = 1 , not = O'); 


% A. SENSING MODE : 

% THE GRAVITY GRADIENT CONTRIBUTION IS SECULAR IN THIS 
MODE FOR STEADY STATE 

% TRACKING ERRORS AS IN OFFSET NADIR POINTING. WERE 
NEGLECTING THE CONTRIBUTION 

% FROM SMALL TRANSIENT ERRORS. USE WORST CASE SCENARIO. 


err=input('steady state offset angles as a row vector; roll pitch yaw in deg.'); 
err=(pi/180)*err; 

h=input('enter orbital altitude in meters'); 

I=input(’enter principal moments as a row; x y z in kg-m A 2’); 

w0=sqrt(3.987464el4/((6378e3 + h) A 3)); 

Mg(l) = (3 * wO A 2)*(I(3)-I(2))*sin(err( 1 ))*cos(err( 1 ))*cos(err(2)); 

Mg(2) = (3*wO A 2)*(I(3)-I(l))*sin(err(2))*cos(err(2))*cos(err(l)); 

Mg(3) = (3*wO A 2)*(I(l)-I(2))*sin(err(2))*cos(err(2))*sin(err(l)); 
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% THE AERODYNAMIC DRAG IS ALSO SECULAR IN THIS MODE. BEGIN 
BY ESTIMATING THE 
% DRAG FORCE: D=Cd*rho/2*v A 2*A 


cd=input('enter drag coefficient'); 

rho=input('enter estimated density at altitude, kg/m A 3’); 

A=input('enter crosssectional area of spacecraft, m A 2'); 


v=(6370e3+h)*w0; 

D=.5*cd*rho*A*v A 2; 

% NEXT DEFINE THE FORCE DIRECTION AND MOMENT ARM. IN THIS 
MODE BOTH ARE FIXED IN % THE BODY NEGLECTING SMALL 
CHANGES DUE TO TRACKING ERRORS. 


Du=input('enter drag direction cosines as a row, x y z'); 
cp=input('enter the cp-cm offset as a row vector; x y z’); 
Md=D.*[0 -cp(3) cp(2);cp(3) 0 -cp(l);-cp(2) cp(l) 0]*Du’; 


% THE SOLAR TORQUES ARE BOTH CYCLIC AND SECULAR IN THIS 
MODE. THEY ARE 

% CALCULATED USING AGRAWAL P ADJUSTED FOR COORDINATE 
SYSTEMS AND SELECTING THE % WORST CASE ORBIT AND SUN 
ANGLE 

% SECULAR TERMS AND CYCLIC AMPLITUDES 
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del=input(’enter sun declination in degrees'); 
del=(pi/180)*del; 

ps=input('enter coefficient of specular reflection'); 
pd=input('diffuse reflection'); 
inc=input('enter orbit inclination in deg.'); 
inc=(pi/180)*inc; 
dele= -(inc-del); 

kl =(( 1 -ps)*cos(dele) + 2*(ps +.33*pd))*cos(dele); 
k2=(l-ps)*cos(dele)*sin(dele); 

PA=4.644e-6*A A 2; 

Mxcos= PA*cp(2)*kl ;Mbx=-PA*cp(3)*k2; 
Mysin=PA*cp<3)*kl;Mycos=-PA*cp(l)*kl; 

Mbz=PA*cp( 1 )*k2;Mzsin=-PA*cp(2)*kl ; 

if sflag==l % COMPUTE TIME HISTORY OF TORQUES AND SUM 

t=input('enter time as a row vector in seconds'); 

Ms(:,l )=Mxcos*cos(wO.*t')+Mbx; 

Ms(:,2)=Mysin*sin(wO.*t')+Mycos.*cos(wO.*t'); 

Ms(:,3)=Mzsin*sin(wO.*t')+Mbz; 

Msum=Mg+Md'; 

Madd(:,l)=Msum(l)*ones(t)';Madd(:,2)=Msum(2)*ones(t)';Madd(:,3)=Msum(3)*o 

nes(t)’;Mtot=Ms+Madd; 

end 
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% NOISE DISTURBANCES ARE CALCULATED AND USED IN 
% ANOTHER ROUTINE DUE TO THEIR VARIABLE AND SPECIAL CASE 
NATURE. BE AWARE 

% THAT Mtot CONTAINS AERO, SOLAR AND GRAVITY GRADIENT 
ONLY. 

% NEXT WE PERFORM SIMILAR ANALYSIS IN THE SUNTRACKING 
MODE, WHERE THE DIS. 

% TORQUES WILL CHANGE DUE TO DIFFERENT ORIENTATION AND 
THE "INERTIAL" NATURE OF % THE COORDINATES 

% SUNTRACKING MODE 

% THE SOLAR TORQUES ARE STRICTLY SECULAR IN THIS MODE 


F=PA*cos(del). * [0;( 1 -ps-pd)*sin(del);(l -ps-pd)*cos(del)+2*ps] ; 
Msst=([0 -cp(3) cp(2);cp(3) 0 -cp(l);-cp(2) cp(l) 0]*F)'; 


% GRAVITY GRADIENT TERMS ARE CYCLIC AND QUITE INVOLVED IN 
THIS MODE AND THE 

% DERIVATION IS GIVEN IN THIS APPENDIX. THERE ARE NO PURELY 
SECULAR TERMS. 


if sflag ==1 


nu=input('enter angle between sunvector projection and ascend, node in deg.'); 
nu=(pi/180)*nu - 90; 
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arg=180-inc; 

ste=input(' small steady state tracking errors in deg. as a row’); 
ste=(pi/180)*ste; 

% CALCULATE AVAILABLE TORQUES FOR MAGNETIC 
DESATURATION ALSO 

myu=input('earth rotation angle wrt ascending node, in deg.'); 
myu=(pi/l 80)*myu; 
bmag=8e 1 5/((6370e3+h) A 3); 

lintran=[l ste(3) -ste(2);-ste(3) 1 ste(l);ste(2) -ste(l) 1]; 
xrot=[l 0 0;0 cos(del) sin(del);0 -sin(del) cos(del)]; 
yrot=[cos(nu) 0 sin(nu);0 1 0;-sin(nu) 0 cos(nu)]; 
xrot2=[l 0 0;0 cos(arg) sin(arg);0 -sin(arg) cos(arg)]; 
ttran=xrot*yrot*xrot2; 

sighz=size(t); 

for k=sighz(l):sighz(2) 

c=cos(wO*t(k)); 

s=sin(wO*t(k)); 

r0=(-(6370e3+h).*ttran'*[c;s;0])’; 

rOt=(lintran*rO')’; 

magr0=(r0t(l) A 2 + r0t(2) A 2 + rOt(3) A 2) A 2.5; 
mul=3*3.987464el4/magrO; 
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Mgst(k,l)=mul*(iOt(2)*rOt(3)*(I(2)-I(3)»; 

Mgst(k,2)=mul*(r0t( 1 )*rOt(3)*(I(3)-I( 1 ))); 

Mgst(k,3)=mul*(rOt( l)*rOt(2)*(I( 1 )-I(2))); 

% NEXT CALCULATE AERO DRAG TERMS WHICH ARE CYCLIC IN THIS 
FRAME. THE DERIVATION % IS GIVEN IN THIS APPENDIX 


d vec=(D . * ttran' * [ s ; -c ; 0] )’ ; 

Mdst(k,:)=([0 -cp(3) cp(2);cp(3) 0 -cp(l);-cp(2) cp(l) 0]*dvec')’; 
Mtotst(k,:)=Mgst(k,:)+Mdst(k,:)+Msst; 

%CONTINUE WITH AVAILABLE DESATURATION TORQUE (PLUS OR 
MINUS) 

Bl=bmag%9816*sin(inc)*c-.1908*sin(myu)*c-.1908*cos(myu)*cos(inc)*c); 

B2=-bmag*(.9816*cos(inc)+.1908*sin(inc)*cos(myu)); 

B3=2*bmag*(.9816*sin(inc)*s+.1908*sin(myu)*c-.1908*cos(myu)*cos(inc)*s); 

Bst=( ttran' *[B1 B2 B3]')'; 

Mmag(k,:)=20*[B3-B2 B1-B3 B2-B1]; 

Mmagst(k,:)=20*[Bst(3)-Bst(2) Bst(l)-Bst(3) Bst(2)-Bst(l)]; 

end 

end 

save dtq t Mg Md Mbx Mbz Ms Mtot Msst Mgst Mdst Mtotst Mmag Mmagst 
clear 
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torque, Nm 


Figure H-l shows the behavior of these disturbance torques for each mode 
over one orbit for a typical orbit - sun orientation. 



time, sec. 


Figure H-l(a) Sense Mode Sum Of Disturbance Torques For One Orbit 

In ATypical Orbit-Sun Orientation 
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Torque, Nm 



Figure H-l(c) Available Fixed Polarity Magnetic Torque Available for 
Desaturation in a Typical Orbit-Sun Orientation 

D. SIZING OF RWAS AND TORQUE RODS 

The required torque producing and angular momentum storage capabilities for 
the RWAs are governed by the slew rate and pointing accuracy requirements. A 
simple model for the slew dynamics is used for the purposes of sizing: 

(H.7) Is A(0 S / At = SM W , 


where the subscript s indicates the properties about the slew axis, the subscript w 
indicates the properties of the wheel, and to is the respective angular velocity . 
When driven by constant torque levels, the angular position and velocity of this 
simplified model will follow parabolic trajectories in the phase plane. This allows 
the designer to obtain average slew rates as a function of constant torque for 
fixed inertia properties. 
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Once the required torque level is determined a 100% margin for design 
uncertainty is added and a search for an off-the-shelf RWA is made that satisfies 
the design torque. In addition to the torque level, these factors must be considered: 
(1) mass, (2) power consumption, and (3) angular momentum storage capacity 
(peak RPM). Designers of the other subsystems need to know about the first two 
and the RWA must be able to absorb the increase in angular momentum over the 
time the slewing torque is applied. 

The final check is in the area of pointing accuracy. Given the selected RWA 
and the disturbance torque models a check is made to see if the RWA can produce 
the required response to maintain pointing accuracies within the specified goal. 
This procedure is explained in the next section. The magnetic torque rods must be 
able to dump the accumulated angular momentum due to disturbance torques. The 
torque produced by a magnetic torque rod is described simply by. 

(H.8) T = m x B, 

where T= torque produced, m= dipole moment of the torque rod, and B = the earth 
magnetic field strength. Using the canted dipole model for the earth magnetic field, 
and the disturbance torque models, the dipole moment size of the torque rods can 
be found as follows: 

(H.9) m x = -(B z M y - B y M z ) / I B 1 2 

(H.10) m y = (B Z M X - B X M Z ) / IbI 2 

(H.ll) m z = -(B y M x - B x M y ) / I B I 2 
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The canted dipole model and detailed derivation of these equations can be 
found in ref (3). Once the dipole moment required is determineed, as in the RWA 
an off-the -shelf component is sought after a design margin of 100% is added. 

E. . CONTROL SYSTEM DESIGN 

The gain Ki, and time constant Ti are determined by the following equations 

found in ref (1): 


(H.12) 

x = 2ale / M , 

(H.13) 

K = I / (x / 2) 2 - 


where 1= the moment of inertia of the axis in question, a = the maximum allowable 
pointing error, e = the natural exponent, and M = the maximum torque impulse 
acting about the axis. These gains are then fed into the following simulation code 
written in MATLAB. It takes the previously calculated disturbance torque models, 
adds in the noise models discussed in chapter VII, and simulates the system for the 
control parameters selected. 
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% THIS PROGRAM DOES THE STEADY STATE POINTING ANALYSIS 
% TO BE CALLED FROM DTORQUE IF NO DTQ FILE RESIDENT 


rand(’normal') 

mflag=input('sense mode=l, strack mode=0'); 
z=input('enter axis no. 1,2 or 3'); 

I=input('enter I,K,tau as a row vector in consistent units’); 
n=input('enter white noise tau and intensity as arow'); 

An=[0 1 0;-(I(2)/I(l)) -(I(2)*I(3)/I(1)) I(2)/I(l);0 0 -l/n(l)]; 
Bn=[0 0;1/I(1)0;0 1]; 

Ab=An(l:2,l:2); 

Bb=Bn( 1:2,1); 

Bnw=Bn(l:3,2); 

cb=eye(Ab);cn=[l 0 0;0 1 0];d=[0 0]';dn=[0 0;0 0]; 


load dtq 

clear Ms Msst Mg Mgst Md Mdst 

% MODE RESPONSES DEPENDING ON FLAG 
% (A) BASIC DTORQUES, NO NOISE 

xO=input('enter initial error in deg. as a row'); 

x0=(pi/180)*x0'; 

if mflag==l 

nflag=input('add a noisy torque? 1-y 0-n’) 
if nflag==l 
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Mtot=Mtot+ 1 e- 8 *rand(Mtot); 
end 

[eab,stb]=lsim(Ab,Bb,cb,d,Mtot(:,z),t,xO); 

else 

[eab,stb]=lsim(Ab,Bb,cb,d>Mtotst(:,z),t,xO); 

end 

% (B) BASIC W/ DESATURATION 


if mflag==l 

M=Mtot+ [Mmag(:,l) -Mmag(:,2) -Mmag(:,3)l; 
Msim=M(:,z); 

[ead,std]=lsim(Ab,Bb,cb,d,Msim,t,xO); 

else 

M=Mtotst-Mmag; 
clear Mmag 
Msim=M(:,z); 
clear M Mtotst 

[ead,std]=lsim(Ab,Bb,cb,d,Msim,t,xO); 

end 

ead=(180/pi)*ead; 
eab=( 1 80/pi) 1|, eab; 

%plot(t,eab); 

%pause 

%plot(t,ead); 
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% WE ARE ONLY INTERESTED IN JITTER EFFECTS FOR THE SENSING 
MODE 


w=n(2)*rand(t); 
if mflag==l 


% NOISE EFFECTS ONLY 

[ean,stn]=Isim(An,Bnw,cn,d,w,t); 

% NOISE PLUS BASIC DTORQUES 

unb=[Mtot(:,z) w']; 
[eanb,stnb]=lsim(An,Bn,cn,dn,unb,t); 

% NOISE PLUS DESATURATION ADDED 

unb=[Msim w']; 

[eand,stnd]=lsim(An,Bn,cn,dn,unb,t); 

ean=(180/pi)*ean; 
eanb=(180/pi)*eanb; 
eand=(l 80/pi)*eand; 


%plot(t,ean(:,2)) 

end 
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This code generated figures 7.3 and 7.4 in chapter VII. Figure H-2 shows 
the control system response to a thruster impulse in the event of torque rod failure 
and in the presence of system noise. 



time, sec. 

Figure H-2 Response to Thruster Impulse as Back Up Desaturation 

The design of a saturating non-linear controller is found in ref (4). The gain 
selection algorithm is given in the following MATLAB code. It is adaptively 
selected from the initial errors and effort desired to steer phase trajectories with 
minimum chatter. 
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% STORES DIFFERENTIAL EQUATIONS TO BE SOLVED 


function xdot=slew(t,x) 


xdot(l)=x(2); 

xdot(2)=kl.*x(4).*x(6) + (Iw/I(l))*(x(6).*x(8)-x(4).*x(9))- ... 

(N(l)/I(l))*sign(x(l) + gain.*x(2)); 
xdot(3)=x(4); 

xdot(4)=k2.*x(2).*x(6) + (Iw/I(2))*(x(2).*x(9)-x(6).*x(7))-... 

(N(2)/I(2))*sign(x(3) + (gain/2).*x(4)); 
xdot(5)=x(6); 

xdot(6)=k3.*x(2).*x(4) + (Iw/I(3))*(x(7).*x(4) - x(2).*x(8))-... 

(N(3)/I(3))*sign(x(5) + (gain/2). *x(6)); 
xdot(7)= -(N( 1)/I( 1 ))*sign(x( 1 ) + gain.*x(2)); 
xdot(8)= -(N(2)/I(2))*sign(x(3) + (gain/2).*x(4 )); 
xdot(9)= -(N(3)/I(3))*sign(x(5) + (gain/2).*x(6)); 

% THIS PROGRAM EXECUTES THE R-K INTEGRATOR 

I=[392.7 397.7 371.1] %input('enter princ. inertias, 1 2 3'); 
Iw=.0884 %input('enter wheel inertia’); 
kl =(I(2)-I(3))/I( 1 );k2=(I(3)-I( 1 ))/I(2);k3=(I( 1 )-I(2))/I(3); 

%hnom=Iw*wnom; 

N=input('enter torque limits, Nm'); 
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global I Iw kl k2 k3 N gain 


t0=0;tf=225; 

x0=[(pi/2) .001 .26 .00005 .26 .00005 0 0 0];gain=22.5; 

[t,x]=ode23(’slew’,tO,tf,xO,le-4); 
x=(l 80/pi). *x; 
x(:,7)=(l/6)*x(:,7); 

plot(t,x(:,l )),title('roll slew angle response'), xlabel('time’) 
pause 

plot(x(:,l),x(:,2)),title(' phase plane response’),xlabel('xl') 
pause 

plot(t,x(:,3)),title('2-axis intermediate response’) 
pause 

plot(t,x(:,5)),title('3-axis intermediate response’) 
pause 

plot(t,x(:,7)),title('wheel speed') 
pause 


Figure H-3 shows the phase plane response of the slew axis. Figure H-4 shows 
the corresponding wheel speed response, and Figure H-5 shows a typical off-axis slew 
response. 
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Figure H-4 Wheel Speed for a Typical Slew 
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APPENDIX I - THERMAL ANALYSIS 

This Mathcad calculation sheet is used to solve thermal balance equations set 
for the spacecraft bus solar arrays. These sheets were done prior to the 
development of the PC-ITAS stand alone thermal model. The sheets included s 
several used in a trade study conducted to see what the temperature of the s< 
arrays would be given certain conditions. Design variations included array 
locations (body mounted .vs. deployed), area of array, density of cell arranger 
coatings to be used on aft side of deployed arrays, and others. The associated 
geometries are given in Figure A8.1 


Geometry A 

This calculation is for a 60 in diam cell layout ignoring emmitance of corners, 
is similar to mounting the cells on the central bay since no radiation of heat 
be allowed into the S/C bus. All values EOL and sun line is normal to surface 


Area (60 in Diam) 

A := 1.824 

Emittance 

E := 0.8 

Absorptance 

a := 0.75 

Cell Packing Factor 

F := 0.9 

Cell Efficiency 

N := 0.15 

Boltzmans Constant 

Sig := 5.67-10 

Solar Flux 

S := 1397 

Effective Solar 
Absorptance 

a := a - N-F 

e 




"a *A S " 


Operating 


L e 


T : = 

(E -A) - Sig 


Temperature. 

op 

- 


T.qp = 370.979 _ 
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Geometry B 

This calculation is for a 60x60 in panel which is covered with a 60 in diam 


Area A := 2.323 


Emittance 

E 

cell 

:= 0.8 

E 

OSR 

:= 0.8 

Absorptance 

a 

cell 

:= 0.75 

a 

OSR 

:= 0.145 

Cell Packing Factor 

F := 0. 

,9 



Cell Efficiency 
Boltzmans Constant 

N := 0. 
Sig : = 

,15 

5.67-10 

8 


Solar Flux 

S := 1397 



Averaged Solar 
Absorptance 


a 

cell 

-1.824 

+ a - . 

OSR 

& 

avg 


2.323 

Effective Averaged 
Solar Absorptance 

a 

eavg 

:= a 

avg 

- N-F 


Operating 

Temperature. 

m • _ 

a 

eavg 

-A-S 

0.25 

T : = 

op 

fE 

_|_ cell 

Aj - Sig 



T 

1 ?? 


= 349.603 
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This calculation is for a 60x60 in panel which is covered with a 60 in diam 
patch of solar cells and the rest covered with OSR, with an additional 10 in 
strip of white paint added to the sides of the S/C looking into space and 
assumed to be totally effective as a radiator (conductively coupled to the 
solar cell area). Values are EOL 


Areas (front and white 
strip) 

Emittances (front and 
white strip) 

Absorptance 


:= 2.323 A 


:= 1.548 


E := 0.8 E := 0.9 

f b 

a := 0.75 a := 0.145 

cell OSR 


Cell Packing Factor 


F := 0.9 


Cell Efficiency 


N := 0.15 


Boltzmans Constant 


Sig := 5.67-10 


Solar Flux 


S := 1397 


Averaged Solar 
Absorptance 


avg 


i -1.824 + a 
cell OSR 

2.323 


.499 


Effective Averaged 
Solar Absorptance 


:= a 


eavg 


- N-F 


Operating 

Temperature. 


i -A S 

eavg f 


E -A + E -A 1 -Sig 

f f b b J 
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= 303.974 
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Geometry D 

This would be the solar array steady state temp if we use the deployed peta 
array configuration ( 4 arrays each 60x12 in ). Front covered with solar cell, I 
painted white. The main advantage to this design is that the arrays are only 
conductively coupled to the S/C . Values again EOL. 


Areas (cells front and 

A 

:= 1.858 A 

:= 1.858 

white back) 

f 

b 


Emittances (front and 

E 

:= 0.8 E 

:= 0.9 

back) 

f 

b 


Absorptance 

a : = 

0.75 


Cell Packing Factor 

F : = 

0.9 


Cell Efficiency 

N : = 

0.15 


Boltzmans Constant 

Sig 

-8 

:= 5.67-10 


Solar Flux 

S : = 

1397 


Effective Averaged 

a 

:= a - N-F 


Solar Absorptance 
Operating 

e 

a -A 

-s 

Temperature. 

T : = 

op 

e f 

[E -A + E 

-A 1 * Sig 


T 

_|_ f f b 

= 307.263 

b J J 


pp 
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Geometry E 


This would be the solar array steady state temp if we use the deployed petal 
solar array configuration ( 4 arrays each 60x15 in ). Front covered with solar 
and 3 inch strip of white paint, back painted white. The main advantage to th 
design is that the arrays are only partially conductively coupled to the S/O . 
Values again EOL. 


Areas (cells/paint 
front and white back) 


A.f := 2.32 A.b := 2.32 


Emittances (front and 
white strip) 

Average Emmitance 
of the Front 


E . cell := 0.82 E. paint := 0.9 


E.eff 


E. cell-1. 858 + E. paint- 0.4 
2.32 


Absorptance 


a. cell := 0.75 a. paint := 0.4 


"a. cell-1. 858 + a. paint-0. 4< 

Average Absorptance a.avg := — “ 

for Front Side 2,32 


Cell Packing Factor 
Cell Efficiency 
Boltzmans Constant 
Solar Flux 


F := 0.9 
N := 0.15 

-8 

Sig := 5.67-10 
S := 1397 


Effective Averaged 
Solar Absorptance 


a 

eavg 


:= a 

avg 


- N-F 
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a -A -S 

Operating Temperature 


eavg f 


T : = 

OP 

[E -A + E -A 1-Sig 


_|_ eff f paint b J 


T.op = 296.626 c 







Figure A8.1 
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This Mathcad calculation sheet is used for estimating the operating 
temperature of the batteries. These sheets were done prior to the developmen 
of the PC-ITAS stand alone thermal model. The sheets were used for initial 
sizing estimates of the battery radiator area and battery heaters. 


Battery Temperature File 


Battery parameters: 


Radiator mounted surface 
0.23x0.32 m 

o 

• 

o 

II 
• • 

H 

4 

Required watts dissipated 
27 watts at 5 deg C 

W1 := 27 
T1 := 283 

OSR radiator parameters: 


Emissivity (EOL) 
Absortivity (EOL) 
Radiator efficiency 

e := 0.8 
a := 0.145 
n := .95 

Other parameters: 


Boltzmans constant 
Sun incident angle 
Solar flux 
Albedo flux 

sig := 5.67 
E := 90 
S := 1397 
A Lb := 79.4 


For a first cut at sizing assume that the radiator is the same size as 
battery face and that it looks into space without any solar inpingemeni 
irradiation by other components . The heat balance equation is then: 


1 


T 


e 


4 


W1 


- sig-n-Al_ 
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T = 353.775 Too Hot! 


APPENDIX I - THERMAL ANALYSIS 


Redefining the problem, what radiator area is required for 54 watts a 
5 deg C? 


A : = 


W1 


[e-sig-n-Tl J 


A = 0.098 


This means that the radiator is 3 times bigger than the conduction 
path area. This can be done but radiator efficiency will be lower so 
try problem with lower assumed efficiency . 


n .85 


A : = 


W1 


4 1 

e-sig-n*Tl 


A = 0.109 „ 

Side := Ja Side = 0.33 t 


The battery heats up during discharge which means that it is in earth 
sensing mode not sun pointing. In earth pointing it is likely that th< 
radiator will be in sun view at least for some part of the orbit 
therfore we need a heat balance equation that includes sun impingement 
(assume a 90 deg incidence) and earth reflection (assume all reflected 
energy hits the radiator) . 


1 


T 


_ L 4 

(W1 + a-A-S-sin(El) + a-A-Alb)] 


(e-sig-n-A) 


El : = 


E 

2 • 

360 


T = 330.74! 

Too Hot! 


If we are at an average theta of only 30 degrees, which is reasonable : 
when in ear th pointing mode we are not inertially fixed, then 


E := 30 
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1 


T 


-L«. 

f(Wl + a- A- S* sin (£l) + a-A-Alb) 1 


(e- sig-n*A) 


T = 310. 87< 

Too Hot! 


Again redefining the problem, what radiator area is required for 54 
watts at 5 deg C with sun impingement at theta = 30 deg ? (even 
smaller efficiency) 

n := .75 


A : = 


wi 


e- sig-n-Tl - a-S-sin(El) - a-AlbJ 


A = 0. 
Side := 
Side = 


Finally we need to consider any other heat sources which may affect th< 
battery such as internal radiation from other components or . conduction 
heat leaks from caiponents mounted on the same panel. Combining these 
make a guess of an additional 50 watts which would need to be dissipat< 
(still maintain the radiator at 5 deg C) . 


Wl + 50 


e-sig-n-Tl - a-Ssin(El) - a-AlbJ A = 0.73 

Side : = \ 

Side = 0. 
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This results in a radiator that is approximately 0.92 square meter (pe: 
battery) . 

Since only the battery is required to be at less than 5 deg C we may si 
able to go with a 0.30 sq. meter radiator if we were to mount the batt< 
that it was adiabatic with respect to the S/C bus. 

Additionally what power is required for the heater to maintain the batteries at 
C in full shadow , ie now other thermal inputs except charging heat plus, done 1 
radiator of 0.256 sq m area 

P := 21.5 

A := 0.256 


T 


1 

4 

' W1 + P 
e- sig- n- A_ 


T = 273.176 
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The PC-ITAS thermal analysis program was used to develop an analytical model of the 
spacecraft. The model consists of user defined geometric nodes and associated internal 
nodes which together represent the spacecraft system. Table A8.1 is a listing of major 
components and the nodes which represent them in the model. 


NODE NUMBER 

COMPONENT 

63 

SPACE SINK 

101 - 106 

PAYLOAD BOX 

199 

PAYLOAD BOX CENTRAL NODE (FOR MASS) 

301 - 304 

+X HONEYCOMB PANEL 

501 - 504 

-X HONEYCOMB PANEL 

701 - 704 

+Y HONEYCOMB PANEL 

901 - 904 

-Y HONEYCOMB PANEL 

1101 - 1104 

PAYLOAD INTERFACE PANEL 

1301 - 1304 

CLOSEOUT BLANKETS ON ANTI EARTH SIDE 

1401 - 1408 

SOLAR ARRAYS 

1199,1299,1399 

CENTER CYLINDER AND INTERNAL HARDWARE 

2201 - 2299 

ANTENNA 

20304 and 20504 

BATTERIES 

20702 and 20902 

EPS ELECTRONICS 

20704 and 20904 

TT&C ELECTRONICS 

All other 20xxx 

VARIOUS LOW DISSIPATION HARDWARE 


TABLE A8.1 
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Additionally values need to be set for a number of parameters. These include: 


1) Orbital parameters 

2) Incident flux settings 

3) Component thermal masses 

4) Optical surface properties 

5) Conductive and radiative conductances 

The parameters are different for the hot (earth sensing) and cold(stand-by) cases. The 
following pages are a summary of parameters used in the development of the PC-ITAS 

thermal model. 

1) Orbital parameters : 

Hot case: 

Earth pointing vehicle 
Non-spinning 
Altitude = 275 NM 

Longitude of ascending node = 90 deg 
Inclination = 70 deg 
Sun inclination = 23 deg 
S/C rotations X=90/Y=00/Z=00 
Cold case: 

Sun pointing vehicle 
Non-spinning 
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Altitude = 275 NM 

Longitude of ascending node = 0 deg 

Inclination =70 deg 

Sun inclination = 23 deg 

S/C rotations X=0/Y =-90/Z=90 

2) Incident flux settings: 

Solar flux hot case = 1340 W/m^ 

Solar flux cold case = 1310 W/m^ 

Albedo hot case: 



Albedo cold case: 



Thermal radiation of earth hot case: 


= 325 


= 306 
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Thermal radiation of earth hot case: 





6 378 2 

6886 2 , 



3) Component thermal masses: 

The inner workings of most electrical boxes and mechanical devices vary in composition. 

It was not possible to determine exact materials used in all components so for developing 
their thermal masses it is assumed that their total mass is made from aluminum. Most 
components use aluminum as part of the housing at the minimum. For a more realistic 
analysis each component thermal mass would need to be calculated individual based on the 
inner workings and construction of the piece. The following formula was used to generate 
our thermal masses: 

, W- min ._ , . , kJ N 1000 

Thermal mass ( — ^ — ) = Mass (kg) • A1 Specific heat ( kg~r°c ' 50 


4) Optical surface properties: 

The PC-ITAS software come with a 350+ catalog of often used surface materials. 
The following materials were used through out the development of the model: 


MATERIAL 

ABSORPTIVITY 

EMISSrVITY 

OPTICAL SOLAR REFLECTORS 

0.08 to 0.12 

0.8 
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MULTI LAYER INSULATION 

0.35 

0.5 

WHITE PAINT 

0.5 

0.8 

BLACK PAINT 

0.94 

0.9 

ANODYZE 

0.36 

0.68 

ALUMINUM 6061-T6 

0.41 

0.04 





TABLE A8.2 


5 ) Conductances: 

In addition to the radiative coupling that occurs between components aboard the 
spacecraft a number of conductive coupling terms have to be defined as well. Among these 
are the panel conductances between adjacent panel nodes, electronics box to mounting 
surface, payload to bus, etc.. Examples of these types of calculations are given below. 

Panel conductances terms: The conductive heat path on a 4 node pane is a function of the 
cross sectional area of the path, the distance between the nodes, and the material of which 
the panel is constructed. For the side honeycomb panels of the bus we have the following 

situation: 
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For conductance path 1 we have, 


Conductance =■ 


K • A 
L 


204 5S~r^c- »5-2cm 2 

50.8 cm 


0.61 


W 

<c 


similarly for path 2 we have, 


Conductance = 


2.04 


W 


cm.°C- 101cm2 n w 
76.2 cm — 0.27 ^ 


assuming a 1 mm skin thickness and ignoring the cell material contribution. 
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Hardware mounting terms: It is assumed that all hardware mounting is done with #10 

screws or 0.25 in bolts. The number of bolts attaching the hardware will depend on the 

mass of the box and was estimated for the model. The conductance is again a function of 

the area of the heat path, the distance traversed, and the bolt material. Some boxes are 

mounted so that one of their sides is flush against the mounting surface. In these cases the 

temperature gradient between the box and the mounting surface would be negligible, this 
means that the conductance term is very high. We used 100 ^ for the conductance which 

is several orders of magnitude higher than if just the bolts were used as the heat paths. 

In situations where the hardware is mounted using stand-offs, or if only a small 
area of the equipment comes in contact with the mounting surface, then the conductance is 
determined by looking at the attachment bolts. Looking at Figure A8.3 
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The following output is generated by the PC-ITAS program the first file is the stand 
alone thermal model output which lists all the information regarding the node network and 
important parameters. The second file is the results file which is a tabular record of the 
results from a steady stste analyses. Tabular hardcopy of transient analyses was not 
available in ascii format. 

These files are representative of the output form PC-ITAS but are not the latest 
versions of the model, and are included for completeness only. 


FILE 1: 


i #OF THERMAL NODES, #OF SURFACES, ORBIT PERIOD(MIN) 
84 84 94.7700 


ACTUAL NODE NUMBERS 






101 

102 

103 

104 

105 



106 

301 

302 

303 

304 

501 


502 

503 

504 

701 

702 

703 


704 

901 

902 

903 

904 

1101 



1102 

1103 

1104 

1301 

1302 

1303 


1304 

1401 

1402 

1403 

1404 

1405 


1406 

1407 

1408 

2201 

2202 

2203 


2204 

2205 

2206 

2207 

2208 

2209 


2210 

2211 

2212 

2213 

2214 

2215 


2216 

2217 

2218 

2219 

2220 

2221 


2222 

2223 

2224 

63 

199 

2299 


20301 

20302 

20303 

20304 

20501 

20502 

20503 

20504 

20701 

20702 

20703 

20704 

20901 

20902 

20903 

20904 

1199 


1299 

1399 

ACTUAL NODE NUMBER POINTERS 





1 

2 

3 4 

5 


6 



7 

8 

9 10 11 

12 


13 

14 

15 

16 

17 


18 


19 

20 

21 

22 

23 


24 


25 

26 

27 

28 

29 


30 


31 

32 

33 

34 

35 


36 
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37 38 39 

43 44 45 

49 50 51 

55 56 57 

61 62 63 

67 68 69 

73 74 75 

79 80 81 


POWR: orbital_absorbed_POWeR 
MASS: thermal_MASS 
COND: linear.CONDuctance 
RAD : RADiative_conductance 
NODE: boundary _NODE 


40 

41 


46 

47 

48 

52 

53 

54 

58 

59 

60 

64 

65 

66 

70 

71 

72 

76 

77 

78 

82 

83 

84 


POWR 


POWR 


101 94.77 110 

0.00 

48.83 

7.90 

48.71 

15.79 

49.45 

23.69 

50.09 

31.59 

142.31 

39.49 

209.12 

47.38 

234.76 

55.28 

216.37 

63.18 

153.65 

71.08 

62.44 

78.97 

60.10 

86.87 

54.32 

94.77 

48.83 

102 94.77 51 . 

0.00 

48.15 

7.90 

48.20 

15.79 

48.18 

23.69 

48.10 

31.59 

48.18 

39.49 

48.20 

47.38 

48.16 

55.28 

51.81 
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63.18 

56.91 

71.08 

58.82 

78.97 

56.91 

86.87 

51.81 

94.77 

48.15 

103 

94.77 459.76 

0.00 

412.25 

7.90 

471.44 

15.79 

514.73 

23.69 

536.07 

31.59 

522.74 

39.49 

476.31 

47.38 

412.66 

55.28 

426.44 

63.18 

438.28 

71.08 

442.52 

78.97 

437.86 

86.87 

425.87 

94.77 

412.25 

104 

94.77 109.51 

0.00 

234.02 

7.90 

208.31 

15.79 

135.72 

23.69 

49.28 

31.59 

48.61 

39.49 

47.84 

47.38 

47.92 

55.28 

53.31 

63.18 

59.09 

71.08 

61.59 

78.97 

152.84 

86.87 

215.62 

94.77 

234.02 

105 

94.77 562.42 

0.00 

558.84 

7.90 

556.14 

15.79 

556.61 

23.69 

556.87 

31.59 

556.69 

39.49 

556.32 
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47.38 

559.11 

55.28 

566.31 

63.18 

571.51 

71.08 

573.31 

78.97 

571.33 

86.87 

566.03 

94.77 

558.84 

106 ( 

94.77 74. 

0.00 

7.46 

7.90 

94.51 

15.79 

51.16 

23.69 

47.66 

31.59 

51.16 

39.49 

94.53 

47.38 

7.49 

55.28 

148.96 

63.18 

81.73 

71.08 

75.30 

78.97 

81.73 

86.87 

148.97 

94.77 

7.46 

301 94.77 37. 

0.00 

14.92 

7.90 

15.31 

15.79 

16.21 

23.69 

17.18 

31.59 

49.42 

39.49 

72.31 

47.38 

80.46 

55.28 

74.02 

63.18 

52.02 

71.08 

20.01 

78.97 

19.12 

86.87 

16.97 

94.77 

14.92 

302 94.77 40.: 

0.00 

15.83 

7.90 

18.48 

15.79 

22.62 

23.69 

27.32 
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31.59 

59.47 

39.49 

79.80 

47.38 

83.77 

55.28 

69.43 

63.18 

47.21 

71.08 

22.53 

78.97 

21.37 

86.87 

18.57 

94.77 

15.83 

303 1 

94.77 

0.00 

14.92 

7.90 

15.31 

15.79 

16.21 

23.69 

17.18 

31.59 

49.42 

39.49 

72.31 

47.38 

80.46 

55.28 

74.02 

63.18 

50.39 

71.08 

20.01 

78.97 

19.12 

86.87 

16.97 

94.77 

14.92 

304 

94.77 

0.00 

22.33 

7.90 

23.55 

15.79 

25.44 

23.69 

27.51 

31.59 

42.28 

39.49 

51.58 

47.38 

53.39 

55.28 

40.36 

63.18 

30.34 

71.08 

25.32 

78.97 

24.86 

86.87 

23.59 

94.77 

22.33 

501 

94.77 

0.00 

80.40 

7.90 

72.22 


37.19 


32.54 


37.25 
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15.79 

49.32 

23.69 

17.10 

31.59 

16.12 

39.49 

15.22 

47.38 

14.82 

55.28 

16.86 

63.18 

19.02 

71.08 

19.95 

78.97 

51.96 

86.87 

73.96 

94.77 

80.40 

502 

94.77 

0.00 

83.77 

7.90 

79.77 

15.79 

59.42 

23.69 

27.29 

31.59 

22.59 

39.49 

18.45 

47.38 

15.80 

55.28 

18.54 

63.18 

21.33 

71.08 

22.52 

78.97 

47.21 

86.87 

69.43 

94.77 

83.77 

503 ! 

94.77 

0.00 

80.45 

7.90 

72.28 

15.79 

49.27 

23.69 

17.16 

31.59 

16.19 

39.49 

15.29 

47.38 

14.90 

55.28 

16.94 

63.18 

19.10 

71.08 

20.00 

78.97 

50.39 

86.87 

74.01 

94.77 

80.45 


504 94.77 32.52 


210 


POWR 


POWR 


APPENDIX I - THERMAL ANALYSIS 


0.00 

53.37 

7.90 

51.55 

15.79 

42.24 

23.69 

27.49 

31.59 

25.40 

39.49 

23.51 

47.38 

22.30 

55.28 

23.55 

63.18 

24.83 

71.08 

25.30 

78.97 

30.32 

86.87 

40.35 

94.77 

53.37 

701 

94.77 196 

0.00 

194.93 

7.90 

194.88 

15.79 

195.72 

23.69 

196.06 

31.59 

195.73 

39.49 

194.89 

47.38 

194.95 

55.28 

197.43 

63.18 

199.22 

71.08 

199.85 

78.97 

199.20 

86.87 

197.41 

94.77 

194.93 

702 

94.77 101 

0.00 

107.87 

7.90 

109.88 

15.79 

111.83 

23.69 

112.56 

31.59 

111.84 

39.49 

109.89 

47.38 

107.87 

55.28 

95.42 

63.18 

86.28 

71.08 

82.82 

78.97 

87.76 

86.87 

96.86 
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94.77 107.87 

703 94.77 196.71 


0.00 

194.94 

7.90 

194.91 

15.79 

195.76 

23.69 

196.10 

31.59 

195.77 

39.49 

194.92 

47.38 

194.96 

55.28 

197.43 

63.18 

199.23 

71.08 

199.86 

78.97 

199.22 

86.87 

197.42 

94.77 

194.94 

704 

94.77 101 

0.00 

107.89 

7.90 

109.92 

15.79 

111.87 

23.69 

112.59 

31.59 

111.87 

39.49 

109.92 

47.38 

107.90 

55.28 

96.90 

63.18 

87.80 

71.08 

82.84 

78.97 

86.31 

86.87 

95.45 

94.77 

107.89 

901 

94.77 15 . 

0.00 

14.66 

7.90 

14.72 

15.79 

14.71 

23.69 

14.64 

31.59 

14.71 

39.49 

14.72 

47.38 

14.66 

55.28 

16.07 

63.18 

17.98 

71.08 

18.69 
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78.97 

17.97 

86.87 

16.07 

94.77 

14.66 

902 

94.77 

0.00 

22.16 

7.90 

22.27 

15.79 

22.26 

23.69 

22.12 

31.59 

22.26 

39.49 

22.27 

47.38 

22.16 

55.28 

23.10 

63.18 

24.25 

71.08 

24.67 

78.97 

24.25 

86.87 

23.10 

94.77 

22.16 

903 

94.77 

0.00 

14.66 

7.90 

14.72 

15.79 

14.71 

23.69 

14.64 

31.59 

14.71 

39.49 

14.72 

47.38 

14.66 

55.28 

16.07 

63.18 

17.98 

71.08 

18.69 

78.97 

17.97 

86.87 

16.07 

94.77 

14.66 

904 

94.77 

0.00 

22.16 

7.90 

22.27 

15.79 

22.26 

23.69 

22.12 

31.59 

22.26 

39.49 

22.27 

47.38 

22.16 

55.28 

23.10 


22.91 


15.80 


22.91 
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63.18 24.25 

71.08 24.67 

78.97 24.25 

86.87 23.10 

94.77 22.16 

1101 94.77 

0.00 0.00 

7.90 0.00 

15.79 0.00 

23.69 0.00 

31.59 0.00 

39.49 0.00 

47.38 0.02 

55.28 48.31 

63.18 83.66 

71.08 96.59 

78.97 83.64 

86.87 48.28 

94.77 0.00 

1102 94.77 

0.00 0.00 

7.90 0.00 

15.79 0.00 

23.69 0.00 

31.59 0.00 

39.49 0.00 

47.38 0.02 

55.28 48.31 

63.18 83.66 

71.08 96.59 

78.97 83.64 

86.87 48.28 

94.77 0.00 

1103 94.77 

0.00 0.00 

7.90 0.00 

15.79 0.00 

23.69 0.00 

31.59 0.00 

39.49 0.00 


30.04 


30.04 


30.04 
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47.38 0.02 

55.28 48.31 

63.18 83.66 

71.08 96.59 

78.97 83.64 

86.87 48.28 

94.77 0.00 

1104 94.77 

0.00 0.00 

7.90 0.00 

15.79 0.00 

23.69 0.00 

31.59 0.00 

39.49 0.00 

47.38 0.02 

55.28 48.31 

63.18 83.66 

71.08 96.59 

78.97 83.64 

86.87 48.28 

94.77 0.00 

1301 94.77 

0.00 3.90 

7.90 42.92 

15.79 31.25 

23.69 29.31 

31.59 28.75 

39.49 41.18 

47.38 3.85 

55.28 26.36 

63.18 16.23 

71.08 15.37 

78.97 16.28 

86.87 26.45 

94.77 3.90 

1302 94.77 

0.00 3.87 

7.90 41.20 

15.79 28.78 

23.69 29.33 


30.04 


23.49 


23.51 
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31.59 31.28 

39.49 42.94 

47.38 3.91 

55.28 26.46 

63.18 16.30 

71.08 15.38 

78.97 16.25 

86.87 26.38 

94.77 3.87 

1303 94.77 

0.00 3.78 

7.90 36.20 

15.79 13.85 

23.69 9.96 

31.59 10.52 

39.49 33.89 

47.38 3.73 

55.28 26.12 

63.18 15.98 

71.08 15.12 

78.97 16.03 

86.87 26.20 

94.77 3.78 

1304 94.77 

0.00 3.74 

7.90 33.88 

15.79 10.52 

23.69 9.96 

31.59 13.85 

39.49 36.20 

47.38 3.78 

55.28 26.20 

63.18 16.03 

71.08 15.11 

78.97 15.98 

86.87 26.13 

94.77 3.74 

1401 94.77 

0.00 96.41 

7.90 118.76 


17.62 


17.62 


132.22 


216 



APPENDIX I - THERMAL ANALYSIS 


POWR 


POWR 


POWR 


15.79 151.03 

23.69 188.52 

31.59 187.22 

39.49 161.80 

47.38 120.29 

55.28 119.86 

63.18 117.26 

71.08 112.69 

78.97 109.56 

86.87 103.20 

94.77 96.41 

1402 94.77 82.37 

0.00 0.01 

7.90 0.01 

15.79 0.01 

23.69 0.01 

31.59 0.01 

39.49 0.01 

47.38 0.05 

55.28 132.44 

63.18 229.35 

71.08 264.81 

78.97 229.31 

86.87 132.38 

94.77 0.01 

1403 94.77 131.16 

0.00 119.34 

7.90 160.82 

15.79 186.16 

23.69 187.39 

31.59 149.94 

39.49 117.66 

47.38 95.32 

55.28 102.07 

63.18 108.47 

71.08 111.60 

78.97 116.29 

86.87 118.91 

94.77 119.34 

1404 94.77 82.37 
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0.00 0.01 

7.90 0.01 

15.79 0.01 

23.69 0.01 

31.59 0.01 

39.49 0.01 

47.38 0.05 

55.28 132.44 

63.18 229.35 

71.08 264.81 

78.97 229.31 

86.87 132.38 

94.77 0.01 

1405 94.77 193.30 

0.00 165.99 

7.90 210.25 

15.79 244.27 

23.69 256.56 

31.59 244.27 

39.49 210.24 

47.38 166.04 

55.28 165.11 

63.18 164.34 

71.08 163.20 

78.97 164.28 

86.87 165.06 

94.77 165.99 

1406 94.77 82.37 

0.00 0.01 

7.90 0.01 

15.79 0.01 

23.69 0.01 

31.59 0.01 

39.49 0.01 

47.38 0.05 

55.28 132.44 

63.18 229.35 

71.08 264.81 

78.97 229.31 

86.87 132.38 
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94.77 

0.01 

1407 

94.77 

0.00 

94.16 

7.90 

96.01 

15.79 

96.05 

23.69 

93.66 

31.59 

96.05 

39.49 

96.01 

47.38 

94.16 

55.28 

100.48 

63.18 

106.76 

71.08 

108.80 

78.97 

106.75 

86.87 

100.48 

94.77 

94.16 

1408 

94.77 

0.00 

0.01 

7.90 

0.01 

15.79 

0.01 

23.69 

0.01 

31.59 

0.01 

39.49 

0.01 

47.38 

0.05 

55.28 

132.44 

63.18 

229.35 

71.08 

264.81 

78.97 

229.31 

86.87 

132.38 

94.77 

0.01 

2201 

94.77 

0.00 

0.34 

7.90 

0.75 

15.79 

1.04 

23.69 

1.19 

31.59 

1.09 

39.49 

0.77 

47.38 

0.33 

55.28 

0.43 

63.18 

0.51 

71.08 

0.54 


99.11 


82.37 


0.66 
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78.97 0.51 

86.87 0.43 

94.77 0.34 

2202 94.77 

0.00 2.57 

7.90 3.45 

15.79 4.38 

23.69 4.69 

31.59 3.75 

39.49 2.88 

47.38 2.56 

55.28 2.77 

63.18 2.95 

71.08 3.02 

78.97 2.96 

86.87 2.78 

94.77 2.57 

2203 94.77 

0.00 0.37 

7.90 0.70 

15.79 0.95 

23.69 1.06 

31.59 0.99 

39.49 0.72 

47.38 1.22 

55.28 2.04 

63.18 2.43 

71.08 2.26 

78.97 1.57 

86.87 0.52 

94.77 0.37 

2204 94.77 

0.00 2.55 

7.90 2.69 

15.79 3.02 

23.69 3.19 

31.59 2.87 

39.49 2.67 

47.38 2.54 

55.28 2.74 


3.23 


1.24 


2.82 
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63.18 2.92 

71.08 2.99 

78.97 2.93 

86.87 2.76 

94.77 2.55 

2205 94.77 

0.00 2.01 

7.90 1.24 

15.79 0.90 

23.69 0.87 

31.59 1.26 

39.49 2.06 

47.38 2.85 

55.28 3.68 

63.18 4.18 

71.08 4.21 

78.97 3.76 

86.87 2.95 

94.77 2.01 

2206 94.77 

0.00 2.54 

7.90 2.67 

15.79 2.80 

23.69 2.85 

31.59 2.80 

39.49 2.67 

47.38 2.53 

55.28 2.72 

63.18 2.89 

71.08 2.96 

78.97 2.90 

86.87 2.73 

94.77 2.54 

2207 94.77 

0.00 2.71 

7.90 1.86 

15.79 1.01 

23.69 0.58 

31.59 0.62 

39.49 1.00 


2.50 


2.75 


2.29 
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47.38 

1.83 

55.28 

2.76 

63.18 

3.56 

71.08 

4.01 

78.97 

4.00 

86.87 

3.52 

94.77 

2.71 

2208 94.77 

0.00 

2.53 

7.90 

2.67 

15.79 

2.80 

23.69 

2.85 

31.59 

2.80 

39.49 

2.67 

47.38 

2.54 

55.28 

2.73 

63.18 

2.90 

71.08 

2.96 

78.97 

2.89 

86.87 

2.72 

94.77 

2.53 

2209 

94.77 

0.00 

1.05 

7.90 

0.29 

15.79 

0.32 

23.69 

0.31 

31.59 

0.26 

39.49 

0.19 

47.38 

0.11 

55.28 

0.22 

63.18 

1.22 

71.08 

1.88 

78.97 

2.11 

86.87 

1.81 

94.77 

1.05 

2210 

94.77 

0.00 

2.54 

7.90 

2.67 

15.79 

2.87 

23.69 

3.19 


2.75 


0.82 


2.82 
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31.59 3.01 

39.49 2.69 

47.38 2.55 

55.28 2.76 

63.18 2.93 

71.08 2.99 

78.97 2.92 

86.87 2.74 

94.77 2.54 

221 1 94.77 

0.00 0.26 

7.90 0.47 

15.79 0.60 

23.69 0.64 

31.59 0.56 

39.49 0.39 

47.38 0.18 

55.28 0.23 

63.18 0.27 

71.08 0.29 

78.97 0.31 

86.87 0.30 

94.77 0.26 

2212 94.77 

0.00 2.56 

7.90 2.88 

15.79 3.75 

23.69 4.69 

31.59 4.38 

39.49 3.45 

47.38 2.57 

55.28 2.78 

63.18 2.96 

71.08 3.02 

78.97 2.95 

86.87 2.77 

94.77 2.56 

2213 94.77 

0.00 0.68 

7.90 1.26 


0.37 


3.23 


1.15 
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15.79 1.69 

23.69 1.90 

31.59 1.76 

39.49 1.30 

47.38 0.68 

55.28 0.82 

63.18 0.95 

71.08 0.99 

78.97 0.94 

86.87 0.81 

94.77 0.68 

2214 94.77 

0.00 2.99 

7.90 7.41 

15.79 10.08 

23.69 9.84 

31.59 8.91 

39.49 6.04 

47.38 2.97 

55.28 3.23 

63.18 3.44 

71.08 3.52 

78.97 3.45 

86.87 3.25 

94.77 2.99 

2215 94.77 

0.00 0.68 

7.90 1.09 

15.79 1.40 

23.69 1.54 

31.59 1.45 

39.49 2.03 

47.38 2.74 

55.28 3.43 

63.18 3.40 

71.08 2.64 

78.97 1.07 

86.87 0.78 

94.77 0.68 

2216 94.77 


5.43 


1.85 


3.56 


224 



APPENDIX I - THERMAL ANALYSIS 


POWR 


POWR 


0.00 

2.96 

7.90 

4.26 

15.79 

4.92 

23.69 

4.47 

31.59 

3.44 

39.49 

3.13 

47.38 

2.93 

55.28 

3.16 

63.18 

3.37 

71.08 

3.46 

78.97 

3.40 

86.87 

3.20 

94.77 

2.96 

2217 

94.77 

0.00 

4.49 

7.90 

3.83 

15.79 

3.63 

23.69 

3.93 

31.59 

4.65 

39.49 

5.60 

47.38 

6.54 

55.28 

7.30 

63.18 

7.57 

71.08 

7.30 

78.97 

6.54 

86.87 

5.51 

94.77 

4.49 

2218 

94.77 

0.00 

2.92 

7.90 

3.12 

15.79 

3.30 

23.69 

3.36 

31.59 

3.29 

39.49 

3.11 

47.38 

2.91 

55.28 

3.11 

63.18 

3.32 

71.08 

3.40 

78.97 

3.33 

86.87 

3.13 


5.57 


3.19 
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94.77 2.92 

2219 94.77 5.26 


0.00 

6.27 

7.90 

5.29 

15.79 

4.30 

23.69 

3.56 

31.59 

3.26 

39.49 

3.49 

47.38 

4.20 

55.28 

5.22 

63.18 

6.24 

71.08 

6.99 

78.97 

7.28 

86.87 

7.01 

94.77 

6.27 

2220 

94.77 

0.00 

2.91 

7.90 

3.11 

15.79 

3.29 

23.69 

3.36 

31.59 

3.30 

39.49 

3.12 

47.38 

2.92 

55.28 

3.13 

63.18 

3.33 

71.08 

3.40 

78.97 

3.32 

86.87 

3.11 

94.77 

2.91 

2221 

94.77 

0.00 

2.39 

7.90 

1.27 

15.79 

0.37 

23.69 

0.35 

31.59 

0.32 

39.49 

0.28 

47.38 

0.26 

55.28 

0.27 

63.18 

0.49 

71.08 

2.04 
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78.97 

2.86 

86.87 

2.99 

94.77 

2.39 

2222 

94.77 

0.00 

2.93 

7.90 

3.13 

15.79 

3.44 

23.69 

4.47 

31.59 

4.92 

39.49 

4.26 

47.38 

2.96 

55.28 

3.20 

63.18 

3.40 

71.08 

3.46 

78.97 

3.37 

86.87 

3.16 

94.77 

2.93 

2223 

94.77 

0.00 

0.49 

7.90 

0.67 

15.79 

0.78 

23.69 

0.79 

31.59 

0.72 

39.49 

0.56 

47.38 

0.37 

55.28 

0.42 

63.18 

0.47 

71.08 

0.49 

78.97 

0.53 

86.87 

0.52 

94.77 

0.49 

2224 

94.77 

0.00 

2.97 

7.90 

6.04 

15.79 

8.91 

23.69 

9.84 

31.59 

10.08 

39.49 

7.41 

47.38 

2.99 

55.28 

3.25 


3.56 


0.57 


5.43 
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63.18 3.45 

71.08 3.52 

78.97 3.44 

86.87 3.23 

94.77 2.97 


MASS 

0 

31.64371000 
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MASS 

0 
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0 
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0 
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RAD 00 0. 1021 829000E+03 2206 

RAD 00 0. 1 1 1 5072000E+03 2207 

RAD 00 0.2086344000E+01 2208 

RAD 00 0.25 0422 1000E+01 2208 

RAD 00 0.3077459000E+01 2208 

RAD 00 0.3025947000E+01 2208 

RAD 00 0.2839348000E+01 2208 

RAD 00 0.2053 159000E+01 2208 

RAD 00 0.2839348000E+01 2208 

RAD 00 0.3025947000E+01 2208 

RAD 00 0.1021829000E+03 2208 

RAD 00 0. 1083426000E+03 2209 


2216 

2218 

2220 

2222 

2224 

63 

63 

2206 

2208 

2210 

2212 

2214 

2216 

2218 

2220 

2222 

2224 

63 

63 

2208 

2210 

2212 

2214 

2216 

2218 

2220 

2222 

2224 

63 

63 

2210 

2212 

2214 

2216 

2218 

2220 

2222 

2224 

63 

63 
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RAD 

00 

0.2086344000E+01 

2210 

2212 

RAD 

00 

0.3025 947000E+01 

2210 

2214 

RAD 

00 

0.3077459000E+01 

2210 

2216 

RAD 

00 

0.3025947000E+01 

2210 

2218 

RAD 

00 

0.2839348000E+01 

2210 

2220 

RAD 

00 

0.2053 159000E+01 

2210 

2222 

RAD 

00 

0.2839348000E+01 

2210 

2224 

RAD 

00 

0.1021829000E+03 

2210 

63 

RAD 

00 

0.9220393000E+02 

2211 

63 

RAD 

00 

0.2839348000E+01 

2212 

2214 

RAD 

00 

0.3025947000E+01 

2212 

2216 

RAD 

00 

0.307745 9000E+01 

2212 

2218 

RAD 

00 

0.3025947000E+01 

2212 

2220 

RAD 

00 

0.2839348000E+01 

2212 

2222 

RAD 

00 

0.2053 159000E+01 

2212 

2224 

RAD 

00 

0. 1 02 1 829000E+03 

2212 

63 

RAD 

00 

0. 1006598000E+03 

2213 

63 

RAD 

00 

0.54845 10000E+01 

2214 

2216 

RAD 

00 

0.4255825000E+01 

2214 

2218 

RAD 

00 

0.4144866000E+01 

2214 

2220 

RAD 

00 

0.4255 825000E+01 

2214 

2222 

RAD 

00 

0.5484494000E+01 

2214 

2224 

RAD 

00 

0. 1 4 1 5008000E+03 

2214 

63 

RAD 

00 

0. 1 1 98694000E+03 

2215 

63 

RAD 

00 

0.54845 10000E+01 

2216 

2218 

RAD 

00 

0.4255825000E+01 

2216 

2220 

RAD 

00 

0.4144866000E+01 

2216 

2222 

RAD 

00 

0.4255825000E+01 

2216 

2224 

RAD 

00 

0.1415008000E+03 

2216 

63 

RAD 

00 

0. 1 590222000E+03 

2217 

63 

RAD 

00 

0.5484494000E+01 

2218 

2220 

RAD 

00 

0.4255825000E+01 

2218 

2222 

RAD 

00 

0.4144866000E+01 

2218 

2224 

RAD 

00 

0.141 5008000E+03 

2218 

63 

RAD 

00 

0.1725270000E+03 

2219 

63 

RAD 

00 

0.54845 10000E+01 

2220 

2222 

RAD 

00 

0.4255825000E+01 

2220 

2224 

RAD 

00 

0.1415008000E+03 

2220 

63 

RAD 

00 

0. 1 690626000E+03 

2221 

63 

RAD 

00 

0.54845 10000E+01 

2222 

2224 


238 



APPENDIX I - THERMAL ANALYSIS 


RAD 

00 

0. 1 4 1 5008000E+03 

2222 

63 

RAD 

00 

0.1423734000E+03 

2223 

63 

RAD 

00 

0.141 5008000E+03 

2224 

63 

NODE 

1 1 

63 1.000 0.001 



$$$$ 

! END OF ITAS-FORMAT THERMAL MATH MODEL 
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FILE 2: 


TIME= 0.00 NODE IR_POWER UV_POWER TOT_POWER 


101 

48.12 

0.71 

48.83 

102 

48.12 

0.03 

48.15 

103 

409.42 

2.83 

412.25 

104 

47.18 

186.84 

234.02 

105 

47.17 

511.68 

558.84 

106 

6.47 

0.99 

7.46 

301 

14.65 

0.27 

14.92 

302 

15.45 

0.38 

15.83 

303 

14.65 

0.27 

14.92 

304 

22.15 

0.17 

22.33 

501 

14.54 

65.86 

80.40 

502 

15.39 

68.38 

83.77 

503 

14.62 

65.83 

80.45 

504 

22.12 

31.25 

53.37 

701 

14.54 

180.39 

194.93 

702 

22.05 

85.81 

107.87 

703 

14.62 

180.33 

194.94 

704 

22.09 

85.80 

107.89 

901 

14.65 

0.02 

14.66 

902 

22.13 

0.03 

22.16 

903 

14.65 

0.02 

14.66 

904 

22.13 

0.03 

22.16 

1101 

0.00 

0.00 

0.00 

1102 

0.00 

0.00 

0.00 

1103 

0.00 

0.00 

0.00 

1104 

0.00 

0.00 

0.00 

1301 

3.51 

0.40 

3.90 

1302 

3.54 

0.33 

3.87 

1303 

3.54 

0.24 

3.78 

1304 

3.54 

0.19 

3.74 

1401 

95.30 

1.12 

96.41 

1402 

0.00 

0.01 

0.01 

1403 

94.07 

25.27 

119.34 

1404 

0.00 

0.01 

0.01 
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1405 

92.59 

73.41 

165.99 

1406 

0.00 

0.01 

0.01 

1407 

93.78 

0.38 

94.16 

1408 

0.00 

0.01 

0.01 

2201 

0.31 

0.03 

0.34 

2202 

2.51 

0.06 

2.57 

2203 

0.25 

0.12 

0.37 

2204 

2.51 

0.04 

2.55 

2205 

0.10 

1.91 

2.01 

2206 

2.51 

0.02 

2.54 

2207 

0.05 

2.66 

2.71 

2208 

2.51 

0.02 

2.53 

2209 

0.06 

0.99 

1.05 

2210 

2.51 

0.03 

2.54 

2211 

0.16 

0.10 

0.26 

2212 

2.51 

0.05 

2.56 

2213 

0.64 

0.03 

0.68 

2214 

2.90 

0.09 

2.99 

2215 

0.51 

0.17 

0.68 

2216 

2.90 

0.06 

2.96 

2217 

0.26 

4.23 

4.49 

2218 

2.90 

0.02 

2.92 

2219 

0.22 

6.05 

6.27 

2220 

2.90 

0.02 

2.91 

2221 

0.23 

2.16 

2.39 

2222 

2.90 

0.03 

2.93 

2223 

0.35 

0.13 

0.49 

2224 

2.90 

0.07 

2.97 


TIME= 7.90 NODE IR_POWER UV_POWER TOT_POWER 


101 

48.14 

0.57 

48.71 

102 

48.14 

0.06 

48.20 

103 

409.41 

62.03 

471.44 

104 

47.20 

161.11 

208.31 

105 

47.18 

508.95 

556.14 

106 

6.47 

88.04 

94.51 

301 

14.65 

0.66 

15.31 

302 

15.45 

3.03 

18.48 

303 

14.65 

0.66 

15.31 

304 

22.16 

1.38 

23.55 

501 

14.54 

57.68 

72.22 
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502 

15.40 

503 

14.63 

504 

22.13 

701 

14.55 

702 

22.06 

703 

14.63 

704 

22.10 

901 

14.65 

902 

22.14 

903 

14.65 

904 

22.14 

1101 

0.00 

1102 

0.00 

1103 

0.00 

1104 

0.00 

1301 

3.51 

1302 

3.55 

1303 

3.55 

1304 

3.55 

1401 

95.22 

1402 

0.00 

1403 

94.01 

1404 

0.00 

1405 

92.60 

1406 

0.00 

1407 

93.78 

1408 

0.00 

2201 

0.31 

2202 

2.51 

2203 

0.25 

2204 

2.51 

2205 

0.10 

2206 

2.51 

2207 

0.05 

2208 

2.51 

2209 

0.06 

2210 

2.51 

2211 

0.16 

2212 

2.51 

2213 

0.64 


64.38 

79.77 

57.65 

72.28 

29.43 

51.55 

180.33 

194.88 

87.82 

109.88 

180.28 

194.91 

87.82 

109.92 

0.06 

14.72 

0.13 

22.27 

0.06 

14.72 

0.13 

22.27 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

39.41 

42.92 

37.65 

41.20 

32.65 

36.20 

30.33 

33.88 

23.55 

118.76 

0.01 

0.01 

66.80 

160.82 

0.01 

0.01 

117.65 

210.25 

0.01 

0.01 

2.23 

96.01 

0.01 

0.01 

0.44 

0.75 

0.94 

3.45 

0.46 

0.70 

0.17 

2.69 

1.15 

1.24 

0.16 

2.67 

1.81 

1.86 

0.16 

2.67 

0.22 

0.29 

0.15 

2.67 

0.31 

0.47 

0.37 

2.88 

0.62 

1.26 
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2214 

2.90 

4.52 

7.41 

2215 

0.51 

0.58 

1.09 

2216 

2.90 

1.37 

4.26 

2217 

0.26 

3.57 

3.83 

2218 

2.90 

0.22 

3.12 

2219 

0.22 

5.07 

5.29 

2220 

2.90 

0.22 

3.11 

2221 

0.23 

1.05 

1.27 

2222 

2.90 

0.23 

3.13 

2223 

0.35 

0.32 

0.67 

2224 

2.90 

3.14 

6.04 


TIME= 15.79 NODE IR_POWER UV_POWER TOT.POWER 


101 

48.13 

1.32 

49.45 

102 

48.13 

0.06 

48.18 

103 

409.41 

105.32 

514.73 

104 

47.19 

88.54 

135.72 

105 

47.17 

509.43 

556.61 

106 

6.47 

44.69 

51.16 

301 

14.65 

1.56 

16.21 

302 

15.45 

7.17 

22.62 

303 

14.65 

1.56 

16.21 

304 

22.16 

3.28 

25.44 

501 

14.55 

34.78 

49.32 

502 

15.39 

44.03 

59.42 

503 

14.62 

34.64 

49.27 

504 

22.11 

20.13 

42.24 

701 

14.54 

181.18 

195.72 

702 

22.05 

89.78 

111.83 

703 

14.62 

181.14 

195.76 

704 

22.09 

89.77 

111.87 

901 

14.65 

0.06 

14.71 

902 

22.13 

0.13 

22.26 

903 

14.65 

0.06 

14.71 

904 

22.13 

0.13 

22.26 

1101 

0.00 

0.00 

0.00 

1102 

0.00 

0.00 

0.00 

1103 

0.00 

0.00 

0.00 

1104 

0.00 

0.00 

0.00 

1301 

3.51 

27.74 

31.25 

1302 

3.54 

25.23 

28.78 
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1303 

3.54 

10.31 

13.85 

1304 

3.54 

6.98 

10.52 

1401 

95.21 

55.83 

151.03 

1402 

0.00 

0.01 

0.01 

1403 

94.01 

92.15 

186.16 

1404 

0.00 

0.01 

0.01 

1405 

92.64 

151.63 

244.27 

1406 

0.00 

0.01 

0.01 

1407 

93.86 

2.19 

96.05 

1408 

0.00 

0.01 

0.01 

2201 

0.31 

0.74 

1.04 

2202 

2.51 

1.87 

4.38 

2203 

0.25 

0.70 

0.95 

2204 

2.51 

0.50 

3.02 

2205 

0.10 

0.80 

0.90 

2206 

2.51 

0.29 

2.80 

2207 

0.05 

0.96 

1.01 

2208 

2.51 

0.29 

2.80 

2209 

0.06 

0.26 

0.32 

2210 

2.51 

0.36 

2.87 

2211 

0.16 

0.44 

0.60 

2212 

2.51 

1.24 

3.75 

2213 

0.64 

1.05 

1.69 

2214 

2.90 

7.18 

10.08 

2215 

0.51 

0.89 

1.40 

2216 

2.90 

2.03 

4.92 

2217 

0.26 

3.36 

3.63 

2218 

2.90 

0.40 

3.30 

2219 

0.22 

4.08 

4.30 

2220 

2.90 

0.39 

3.29 

2221 

0.23 

0.14 

0.37 

2222 

2.90 

0.55 

3.44 

2223 

0.35 

0.43 

0.78 

2224 

2.90 

6.02 

8.91 


TIME= 23.69 NODE IR_POWER UV_POWER TOT_POWER 


101 

47.87 

2.22 

50.09 

102 

48.09 

0.01 

48.10 

103 

409.43 

126.64 

536.07 

104 

46.97 

2.31 

49.28 

105 

47.15 

509.72 

556.87 
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106 

6.44 

301 

14.57 

302 

15.38 

303 

14.57 

304 

22.05 

501 

14.48 

502 

15.33 

503 

14.55 

504 

22.02 

701 

14.53 

702 

22.05 

703 

14.61 

704 

22.08 

901 

14.64 

902 

22.12 

903 

14.64 

904 

22.12 

1101 

0.00 

1102 

0.00 

1103 

0.00 

1104 

0.00 

1301 

3.50 

1302 

3.53 

1303 

3.53 

1304 

3.53 

1401 

95.64 

1402 

0.00 

1403 

94.37 

1404 

0.00 

1405 

92.45 

1406 

0.00 

1407 

93.64 

1408 

0.00 

2201 

0.31 

2202 

2.51 

2203 

0.25 

2204 

2.51 

2205 

0.10 

2206 

2.51 

2207 

0.05 


41.22 

47.66 

2.60 

17.18 

11.94 

27.32 

2.60 

17.18 

5.46 

27.51 

2.61 

17.10 

11.95 

27.29 

2.60 

17.16 

5.46 

27.49 

181.53 

196.06 

90.51 

112.56 

181.49 

196.10 

90.51 

112.59 

0.00 

14.64 

0.00 

22.12 

0.00 

14.64 

0.00 

22.12 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

25.81 

29.31 

25.80 

29.33 

6.44 

9.96 

6.44 

9.96 

92.88 

188.52 

0.01 

0.01 

93.02 

187.39 

0.01 

0.01 

164.11 

256.56 

0.01 

0.01 

0.02 

93.66 

0.01 

0.01 

0.88 

1.19 

2.18 

4.69 

0.82 

1.06 

0.67 

3.19 

0.77 

0.87 

0.33 

2.85 

0.53 

0.58 
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2208 

2.51 

0.33 

2.85 

2209 

0.06 

0.24 

0.31 

2210 

2.51 

0.67 

3.19 

2211 

0.16 

0.48 

0.64 

2212 

2.51 

2.18 

4.69 

2213 

0.64 

1.25 

1.90 

2214 

2.90 

6.94 

9.84 

2215 

0.51 

1.03 

1.54 

2216 

2.90 

1.57 

4.47 

2217 

0.26 

3.67 

3.93 

2218 

2.90 

0.46 

3.36 

2219 

0.22 

3.34 

3.56 

2220 

2.90 

0.46 

3.36 

2221 

0.23 

0.12 

0.35 

2222 

2.90 

1.57 

4.47 

2223 

0.35 

0.44 

0.79 

2224 

2.90 

6.94 

9.84 


TIME= 31.59 NODE IR_POWER UV_POWER TOT_POWER 


101 

48.13 

94.18 

142.31 

102 

48.13 

0.06 

48.18 

103 

409.41 

113.33 

522.74 

104 

47.18 

1.42 

48.61 

105 

47.17 

509.51 

556.69 

106 

6.47 

44.69 

51.16 

301 

14.65 

34.77 

49.42 

302 

15.45 

44.02 

59.47 

303 

14.65 

34.77 

49.42 

304 

22.16 

20.12 

42.28 

501 

14.55 

1.58 

16.12 

502 

15.40 

7.19 

22.59 

503 

14.63 

1.57 

16.19 

504 

22.12 

3.29 

25.40 

701 

14.54 

181.19 

195.73 

702 

22.06 

89.78 

111.84 

703 

14.62 

181.15 

195.77 

704 

22.09 

89.77 

111.87 

901 

14.65 

0.06 

14.71 

902 

22.14 

0.13 

22.26 

903 

14.65 

0.06 

14.71 

904 

22.14 

0.13 

22.26 
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1101 

0.00 

0.00 

0.00 

1102 

0.00 

0.00 

0.00 

1103 

0.00 

0.00 

0.00 

1104 

0.00 

0.00 

0.00 

1301 

3.51 

25.25 

28.75 

1302 

3.54 

27.73 

31.28 

1303 

3.54 

6.98 

10.52 

1304 

3.54 

10.31 

13.85 

1401 

95.21 

92.01 

187.22 

1402 

0.00 

0.01 

0.01 

1403 

93.99 

55.95 

149.94 

1404 

0.00 

0.01 

0.01 

1405 

92.64 

151.63 

244.27 

1406 

0.00 

0.01 

0.01 

1407 

93.86 

2.19 

96.05 

1408 

0.00 

0.01 

0.01 

2201 

0.31 

0.78 

1.09 

2202 

2.51 

1.24 

3.75 

2203 

0.25 

0.74 

0.99 

2204 

2.51 

0.36 

2.87 

2205 

0.10 

1.17 

1.26 

2206 

2.51 

0.29 

2.80 

2207 

0.05 

0.57 

0.62 

2208 

2.51 

0.29 

2.80 

2209 

0.06 

0.20 

0.26 

2210 

2.51 

0.50 

3.01 

2211 

0.16 

0.40 

0.56 

2212 

2.51 

1.87 

4.38 

2213 

0.64 

1.12 

1.76 

2214 

2.90 

6.02 

8.91 

2215 

0.51 

0.94 

1.45 

2216 

2.90 

0.55 

3.44 

2217 

0.26 

4.39 

4.65 

2218 

2.90 

0.39 

3.29 

2219 

0.22 

3.04 

3.26 

2220 

2.90 

0.40 

3.30 

2221 

0.23 

0.09 

0.32 

2222 

2.90 

2.03 

4.92 

2223 

0.35 

0.37 

0.72 

2224 

2.90 

7.18 

10.08 
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TIME= 39.49 NODE IR_POWER UV.POWER TOT.POWER 


101 

48.14 

160.98 

209.12 

102 

48.14 

0.06 

48.20 

103 

409.40 

66.91 

476.31 

104 

47.20 

0.65 

47.84 

105 

47.19 

509.13 

556.32 

106 

6.48 

88.05 

94.53 

301 

14.65 

57.65 

72.31 

302 

15.45 

64.35 

79.80 

303 

14.65 

57.65 

72.31 

304 

22.16 

29.41 

51.58 

501 

14.54 

0.68 

15.22 

502 

15.40 

3.05 

18.45 

503 

14.63 

0.67 

15.29 

504 

22.12 

1.39 

23.51 

701 

14.55 

180.34 

194.89 

702 

22.07 

87.83 

109.89 

703 

14.63 

180.29 

194.92 

704 

22.10 

87.82 

109.92 

901 

14.65 

0.06 

14.72 

902 

22.14 

0.13 

22.27 

903 

14.65 

0.06 

14.72 

904 

22.14 

0.13 

22.27 

1101 

0.00 

0.00 

0.00 

1102 

0.00 

0.00 

0.00 

1103 

0.00 

0.00 

0.00 

1104 

0.00 

0.00 

0.00 

1301 

3.51 

37.67 

41.18 

1302 

3.55 

39.39 

42.94 

1303 

3.55 

30.34 

33.89 

1304 

3.55 

32.66 

36.20 

1401 

95.22 

66.58 

161.80 

1402 

0.00 

0.01 

0.01 

1403 

93.98 

23.68 

117.66 

1404 

0.00 

0.01 

0.01 

1405 

92.57 

117.67 

210.24 

1406 

0.00 

0.01 

0.01 

1407 

93.78 

2.23 

96.01 

1408 

0.00 

0.01 

0.01 

2201 

0.31 

0.46 

0.77 
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2202 

2.51 

0.37 

2.88 

2203 

0.25 

0.48 

0.72 

2204 

2.51 

0.15 

2.67 

2205 

0.10 

1.96 

2.06 

2206 

2.51 

0.16 

2.67 

2207 

0.05 

0.95 

1.00 

2208 

2.51 

0.16 

2.67 

2209 

0.06 

0.13 

0.19 

2210 

2.51 

0.17 

2.69 

2211 

0.16 

0.23 

0.39 

2212 

2.51 

0.94 

3.45 

2213 

0.64 

0.66 

1.30 

2214 

2.90 

3.14 

6.04 

2215 

0.51 

1.52 

2.03 

2216 

2.90 

0.23 

3.13 

2217 

0.26 

5.34 

5.60 

2218 

2.90 

0.22 

3.11 

2219 

0.22 

3.27 

3.49 

2220 

2.90 

0.22 

3.12 

2221 

0.23 

0.05 

0.28 

2222 

2.90 

1.36 

4.26 

2223 

0.35 

0.21 

0.56 

2224 

2.90 

4.51 

7.41 


TIME= 47.38 NODE IR_POWER UV_POWER TOT_POWER 


101 

48.12 

186.64 

234.76 

102 

48.12 

0.03 

48.16 

103 

409.42 

3.24 

412.66 

104 

47.18 

0.74 

47.92 

105 

47.17 

511.94 

559.11 

106 

6.47 

1.03 

7.49 

301 

14.65 

65.82 

80.46 

302 

15.45 

68.32 

83.77 

303 

14.65 

65.82 

80.46 

304 

22.15 

31.23 

53.39 

501 

14.54 

0.29 

14.82 

502 

15.40 

0.41 

15.80 

503 

14.62 

0.27 

14.90 

504 

22.12 

0.18 

22.30 

701 

14.54 

180.41 

194.95 

702 

22.06 

85.81 

107.87 
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703 

14.62 

180.34 

194.96 

704 

22.09 

85.81 

107.90 

901 

14.65 

0.02 

14.66 

902 

22.13 

0.03 

22.16 

903 

14.65 

0.02 

14.66 

904 

22.13 

0.03 

22.16 

1101 

0.00 

0.02 

0.02 

1102 

0.00 

0.02 

0.02 

1103 

0.00 

0.02 

0.02 

1104 

0.00 

0.02 

0.02 

1301 

3.51 

0.34 

3.85 

1302 

3.54 

0.37 

3.91 

1303 

3.54 

0.18 

3.73 

1304 

3.54 

0.23 

3.78 

1401 

95.30 

24.99 

120.29 

1402 

0.00 

0.05 

0.05 

1403 

94.06 

1.26 

95.32 

1404 

0.00 

0.05 

0.05 

1405 

92.58 

73.46 

166.04 

1406 

0.00 

0.05 

0.05 

1407 

93.78 

0.38 

94.16 

1408 

0.00 

0.05 

0.05 

2201 

0.31 

0.02 

0.33 

2202 

2.51 

0.05 

2.56 

2203 

0.25 

0.97 

1.22 

2204 

2.51 

0.03 

2.54 

2205 

0.10 

2.75 

2.85 

2206 

2.51 

0.02 

2.53 

2207 

0.05 

1.78 

1.83 

2208 

2.51 

0.02 

2.54 

2209 

0.06 

0.05 

0.11 

2210 

2.51 

0.04 

2.55 

2211 

0.16 

0.02 

0.18 

2212 

2.51 

0.06 

2.57 

2213 

0.64 

0.03 

0.68 

2214 

2.90 

0.07 

2.97 

2215 

0.51 

2.23 

2.74 

2216 

2.90 

0.03 

2.93 

2217 

0.26 

6.28 

6.54 

2218 

2.90 

0.02 

2.91 
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2219 

0.22 

3.97 

4.20 

2220 

2.90 

0.02 

2.92 

2221 

0.23 

0.03 

0.26 

2222 

2.90 

0.06 

2.96 

2223 

0.35 

0.01 

0.37 

2224 

2.90 

0.09 

2.99 


TIME= 55.28 NODE IR_POWER UV.POWER TOT_POWER 


101 

48.14 

168.23 

216.37 

102 

48.14 

3.66 

51.81 

103 

409.40 

17.03 

426.44 

104 

47.20 

6.11 

53.31 

105 

47.19 

519.12 

566.31 

106 

6.48 

142.48 

148.96 

301 

14.65 

59.37 

74.02 

302 

15.45 

53.98 

69.43 

303 

14.65 

59.37 

74.02 

304 

22.16 

18.20 

40.36 

501 

14.54 

2.32 

16.86 

502 

15.40 

3.14 

18.54 

503 

14.63 

2.32 

16.94 

504 

22.12 

1.43 

23.55 

701 

14.55 

182.88 

197.43 

702 

22.07 

73.35 

95.42 

703 

14.63 

182.81 

197.43 

704 

22.10 

74.80 

96.90 

901 

14.65 

1.41 

16.07 

902 

22.14 

0.95 

23.10 

903 

14.65 

1.41 

16.07 

904 

22.14 

0.95 

23.10 

1101 

0.00 

48.31 

48.31 

1102 

0.00 

48.31 

48.31 

1103 

0.00 

48.31 

48.31 

1104 

0.00 

48.31 

48.31 

1301 

3.51 

22.85 

26.36 

1302 

3.55 

22.92 

26.46 

1303 

3.55 

22.57 

26.12 

1304 

3.55 

22.65 

26.20 

1401 

95.22 

24.64 

119.86 

1402 

0.00 

132.44 

132.44 

1403 

93.98 

8.09 

102.07 
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TIME= 


1404 

0.00 

132.44 

132.44 

1405 

92.57 

72.55 

165.11 

1406 

0.00 

132.44 

132.44 

1407 

93.78 

6.71 

100.48 

1408 

0.00 

132.44 

132.44 

2201 

0.31 

0.12 

0.43 

2202 

2.51 

0.26 

2.77 

2203 

0.25 

1.80 

2.04 

2204 

2.51 

0.23 

2.74 

2205 

0.10 

3.59 

3.68 

2206 

2.51 

0.20 

2.72 

2207 

0.05 

2.71 

2.76 

2208 

2.51 

0.21 

2.73 

2209 

0.06 

0.15 

0.22 

2210 

2.51 

0.25 

2.76 

2211 

0.16 

0.07 

0.23 

2212 

2.51 

0.27 

2.78 

2213 

0.64 

0.17 

0.82 

2214 

2.90 

0.33 

3.23 

2215 

0.51 

2.92 

3.43 

2216 

2.90 

0.26 

3.16 

2217 

0.26 

7.03 

7.30 

2218 

2.90 

0.21 

3.11 

2219 

0.22 

5.00 

5.22 

2220 

2.90 

0.24 

3.13 

2221 

0.23 

0.04 

0.27 

2222 

2.90 

0.31 

3.20 

2223 

0.35 

0.06 

0.42 

2224 

2.90 

0.35 

3.25 


63.18 NODE IR_POWER UV_POWER TOT_POWER 


101 

48.13 

105.52 

153.65 

102 

48.13 

8.79 

56.91 

103 

409.41 

28.87 

438.28 

104 

47.18 

11.91 

59.09 

105 

47.17 

524.33 

571.51 

106 

6.47 

75.25 

81.73 

301 

14.65 

37.37 

52.02 

302 

15.45 

31.76 

47.21 

303 

14.65 

35.74 

50.39 

304 

22.16 

8.19 

30.34 
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501 

14.55 

4.48 

19.02 

502 

15.40 

5.93 

21.33 

503 

14.63 

4.47 

19.10 

504 

22.12 

2.71 

24.83 

701 

14.54 

184.68 

199.22 

702 

22.06 

64.22 

86.28 

703 

14.62 

184.61 

199.23 

704 

22.09 

65.71 

87.80 

901 

14.65 

3.33 

17.98 

902 

22.13 

2.11 

24.25 

903 

14.65 

3.33 

17.98 

904 

22.13 

2.11 

24.25 

1101 

0.00 

83.65 

83.66 

1102 

0.00 

83.65 

83.66 

1103 

0.00 

83.65 

83.66 

1104 

0.00 

83.65 

83.66 

1301 

3.51 

12.72 

16.23 

1302 

3.54 

12.76 

16.30 

1303 

3.54 

12.44 

15.98 

1304 

3.54 

12.49 

16.03 

1401 

95.21 

22.05 

117.26 

1402 

0.00 

229.35 

229.35 

1403 

93.99 

14.48 

108.47 

1404 

0.00 

229.35 

229.35 

1405 

92.64 

71.70 

164.34 

1406 

0.00 

229.35 

229.35 

1407 

93.86 

12.90 

106.76 

1408 

0.00 

229.35 

229.35 

2201 

0.31 

0.20 

0.51 

2202 

2.51 

0.44 

2.95 

2203 

0.25 

2.18 

2.43 

2204 

2.51 

0.41 

2.92 

2205 

0.10 

4.08 

4.18 

2206 

2.51 

0.38 

2.89 

2207 

0.05 

3.51 

3.56 

2208 

2.51 

0.39 

2.90 

2209 

0.06 

1.16 

1.22 

2210 

2.51 

0.42 

2.93 

2211 

0.16 

0.11 

0.27 

2212 

2.51 

0.44 

2.96 
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2213 

0.64 

0.30 

0.95 

2214 

2.90 

0.54 

3.44 

2215 

0.51 

2.89 

3.40 

2216 

2.90 

0.48 

3.37 

2217 

0.26 

7.31 

7.57 

2218 

2.90 

0.42 

3.32 

2219 

0.22 

6.02 

6.24 

2220 

2.90 

0.44 

3.33 

2221 

0.23 

0.26 

0.49 

2222 

2.90 

0.50 

3.40 

2223 

0.35 

0.12 

0.47 

2224 

2.90 

0.56 

3.45 


TIME= 71.08 NODE IR_POWER UV.POWER TOT_POWER 


101 

47.87 

14.57 

62.44 

102 

48.09 

10.73 

58.82 

103 

409.43 

33.09 

442.52 

104 

46.97 

14.62 

61.59 

105 

47.15 

526.16 

573.31 

106 

6.44 

68.86 

75.30 

301 

14.57 

5.43 

20.01 

302 

15.38 

7.15 

22.53 

303 

14.57 

5.43 

20.01 

304 

22.05 

3.27 

25.32 

501 

14.48 

5.47 

19.95 

502 

15.33 

7.19 

22.52 

503 

14.55 

5.45 

20.00 

504 

22.02 

3.28 

25.30 

701 

14.53 

185.32 

199.85 

702 

22.05 

60.77 

82.82 

703 

14.61 

185.25 

199.86 

704 

22.08 

60.76 

82.84 

901 

14.64 

4.06 

18.69 

902 

22.12 

2.55 

24.67 

903 

14.64 

4.06 

18.69 

904 

22.12 

2.55 

24.67 

1101 

0.00 

96.59 

96.59 

1102 

0.00 

96.59 

96.59 

1103 

0.00 

96.59 

96.59 

1104 

0.00 

96.59 

96.59 

1301 

3.50 

11.88 

15.37 
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1302 

3.53 

11.86 

15.38 

1303 

3.53 

11.59 

15.12 

1304 

3.53 

11.59 

15.11 

1401 

95.64 

17.06 

112.69 

1402 

0.00 

264.81 

264.81 

1403 

94.37 

17.23 

111.60 

1404 

0.00 

264.81 

264.81 

1405 

92.45 

70.76 

163.20 

1406 

0.00 

264.81 

264.81 

1407 

93.64 

15.16 

108.80 

1408 

0.00 

264.81 

264.81 

2201 

0.31 

0.23 

0.54 

2202 

2.51 

0.51 

3.02 

2203 

0.25 

2.01 

2.26 

2204 

2.51 

0.48 

2.99 

2205 

0.10 

4.11 

4.21 

2206 

2.51 

0.45 

2.96 

2207 

0.05 

3.96 

4.01 

2208 

2.51 

0.45 

2.96 

2209 

0.06 

1.82 

1.88 

2210 

2.51 

0.48 

2.99 

2211 

0.16 

0.13 

0.29 

2212 

2.51 

0.51 

3.02 

2213 

0.64 

0.35 

0.99 

2214 

2.90 

0.63 

3.52 

2215 

0.51 

2.13 

2.64 

2216 

2.90 

0.57 

3.46 

2217 

0.26 

7.03 

7.30 

2218 

2.90 

0.50 

3.40 

2219 

0.22 

6.77 

6.99 

2220 

2.90 

0.50 

3.40 

2221 

0.23 

1.81 

2.04 

2222 

2.90 

0.57 

3.46 

2223 

0.35 

0.14 

0.49 

2224 

2.90 

0.63 

3.52 


TIME= 78.97 NODE IR_POWER UV_POWER TOT_POWER 


101 

48.13 

11.97 

60.10 

102 

48.13 

8.78 

56.91 

103 

409.41 

28.44 

437.86 

104 

47.19 

105.66 

152.84 
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105 

47.17 

106 

6.47 

301 

14.65 

302 

15.45 

303 

14.65 

304 

22.16 

501 

14.55 

502 

15.39 

503 

14.62 

504 

22.11 

701 

14.54 

702 

22.05 

703 

14.62 

704 

22.10 

901 

14.65 

902 

22.14 

903 

14.65 

904 

22.14 

1101 

0.00 

1102 

0.00 

1103 

0.00 

1104 

0.00 

1301 

3.51 

1302 

3.54 

1303 

3.54 

1304 

3.54 

1401 

95.21 

1402 

0.00 

1403 

94.01 

1404 

0.00 

1405 

92.64 

1406 

0.00 

1407 

93.86 

1408 

0.00 

2201 

0.31 

2202 

2.51 

2203 

0.25 

2204 

2.51 

2205 

0.10 

2206 

2.51 


524.15 

571.33 

75.26 

81.73 

4.47 

19.12 

5.92 

21.37 

4.47 

19.12 

2.70 

24.86 

37.41 

51.96 

31.81 

47.21 

35.77 

50.39 

8.20 

30.32 

184.67 

199.20 

65.71 

87.76 

184.60 

199.22 

64.21 

86.31 

3.33 

17.97 

2.11 

24.25 

3.33 

17.97 

2.11 

24.25 

83.64 

83.64 

83.64 

83.64 

83.64 

83.64 

83.64 

83.64 

12.77 

16.28 

12.71 

16.25 

12.49 

16.03 

12.44 

15.98 

14.36 

109.56 

229.31 

229.31 

22.28 

116.29 

229.31 

229.31 

71.64 

164.28 

229.31 

229.31 

12.90 

106.75 

229.31 

229.31 

0.20 

0.51 

0.44 

2.96 

1.33 

1.57 

0.42 

2.93 

3.66 

3.76 

0.38 

2.90 
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2207 

0.05 

3.95 

4.00 

2208 

2.51 

0.38 

2.89 

2209 

0.06 

2.05 

2.11 

2210 

2.51 

0.41 

2.92 

2211 

0.16 

0.15 

0.31 

2212 

2.51 

0.44 

2.95 

2213 

0.64 

0.30 

0.94 

2214 

2.90 

0.56 

3.45 

2215 

0.51 

0.55 

1.07 

2216 

2.90 

0.50 

3.40 

2217 

0.26 

6.28 

6.54 

2218 

2.90 

0.43 

3.33 

2219 

0.22 

7.06 

7.28 

2220 

2.90 

0.42 

3.32 

2221 

0.23 

2.63 

2.86 

2222 

2.90 

0.48 

3.37 

2223 

0.35 

0.17 

0.53 

2224 

2.90 

0.54 

3.44 


TIME= 86.87 NODE IR_POWER UV.POWER TOT_POWER 


101 

48.14 

6.18 

54.32 

102 

48.14 

3.66 

51.81 

103 

409.41 

16.46 

425.87 

104 

47.20 

168.42 

215.62 

105 

47.18 

518.84 

566.03 

106 

6.47 

142.50 

148.97 

301 

14.65 

2.32 

16.97 

302 

15.45 

3.12 

18.57 

303 

14.65 

2.32 

16.97 

304 

22.16 

1.43 

23.59 

501 

14.54 

59.42 

73.96 

502 

15.40 

54.04 

69.43 

503 

14.63 

59.39 

74.01 

504 

22.13 

18.23 

40.35 

701 

14.55 

182.86 

197.41 

702 

22.06 

74.80 

96.86 

703 

14.63 

182.79 

197.42 

704 

22.10 

73.34 

95.45 

901 

14.65 

1.41 

16.07 

902 

22.14 

0.95 

23.10 

903 

14.65 

1.41 

16.07 
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904 

22.14 

1101 

0.00 

1102 

0.00 

1103 

0.00 

1104 

0.00 

1301 

3.51 

1302 

3.55 

1303 

3.55 

1304 

3.55 

1401 

95.22 

1402 

0.00 

1403 

94.01 

1404 

0.00 

1405 

92.60 

1406 

0.00 

1407 

93.78 

1408 

0.00 

2201 

0.31 

2202 

2.51 

2203 

0.25 

2204 

2.51 

2205 

0.10 

2206 

2.51 

2207 

0.05 

2208 

2.51 

2209 

0.06 

2210 

2.51 

2211 

0.16 

2212 

2.51 

2213 

0.64 

2214 

2.90 

2215 

0.51 

2216 

2.90 

2217 

0.26 

2218 

2.90 

2219 

0.22 

2220 

2.90 

2221 

0.23 

2222 

2.90 

2223 

0.35 


0.95 

23.10 

48.28 

48.28 

48.28 

48.28 

48.28 

48.28 

48.28 

48.28 

22.94 

26.45 

22.84 

26.38 

22.66 

26.20 

22.58 

26.13 

7.98 

103.20 

132.37 

132.38 

24.90 

118.91 

132.37 

132.38 

72.46 

165.06 

132.37 

132.38 

6.70 

100.48 

132.37 

132.38 

0.12 

0.43 

0.27 

2.78 

0.28 

0.52 

0.25 

2.76 

2.85 

2.95 

0.21 

2.73 

3.47 

3.52 

0.20 

2.72 

1.75 

1.81 

0.23 

2.74 

0.14 

0.30 

0.26 

2.77 

0.17 

0.81 

0.35 

3.25 

0.27 

0.78 

0.31 

3.20 

5.25 

5.51 

0.24 

3.13 

6.79 

7.01 

0.21 

3.11 

2.76 

2.99 

0.26 

3.16 

0.17 

0.52 
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2224 2.90 0.33 3.23 

TIME= 94.77 NODE IR_POWER UV.POWER TOT_POWER 


101 

48.12 

0.71 

48.83 

102 

48.12 

0.03 

48.15 

103 

409.42 

2.83 

412.25 

104 

47.18 

186.84 

234.02 

105 

47.17 

511.68 

558.84 

106 

6.47 

0.99 

7.46 

301 

14.65 

0.27 

14.92 

302 

15.45 

0.38 

15.83 

303 

14.65 

0.27 

14.92 

304 

22.15 

0.17 

22.33 

501 

14.54 

65.86 

80.40 

502 

15.39 

68.38 

83.77 

503 

14.62 

65.83 

80.45 

504 

22.12 

31.25 

53.37 

701 

14.54 

180.39 

194.93 

702 

22.05 

85.81 

107.87 

703 

14.62 

180.33 

194.94 

704 

22.09 

85.80 

107.89 

901 

14.65 

0.02 

14.66 

902 

22.13 

0.03 

22.16 

903 

14.65 

0.02 

14.66 

904 

22.13 

0.03 

22.16 

1101 

0.00 

0.00 

0.00 

1102 

0.00 

0.00 

0.00 

1103 

0.00 

0.00 

0.00 

1104 

0.00 

0.00 

0.00 

1301 

3.51 

0.40 

3.90 

1302 

3.54 

0.33 

3.87 

1303 

3.54 

0.24 

3.78 

1304 

3.54 

0.19 

3.74 

1401 

95.30 

1.12 

96.41 

1402 

0.00 

0.01 

0.01 

1403 

94.07 

25.27 

119.34 

1404 

0.00 

0.01 

0.01 

1405 

92.59 

73.41 

165.99 

1406 

0.00 

0.01 

0.01 

1407 

93.78 

0.38 

94.16 

1408 

0.00 

0.01 

0.01 
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2201 

0.31 

0.03 

0.34 

2202 

2.51 

0.06 

2.57 

2203 

0.25 

0.12 

0.37 

2204 

2.51 

0.04 

2.55 

2205 

0.10 

1.91 

2.01 

2206 

2.51 

0.02 

2.54 

2207 

0.05 

2.66 

2.71 

2208 

2.51 

0.02 

2.53 

2209 

0.06 

0.99 

1.05 

2210 

2.51 

0.03 

2.54 

2211 

0.16 

0.10 

0.26 

2212 

2.51 

0.05 

2.56 

2213 

0.64 

0.03 

0.68 

2214 

2.90 

0.09 

2.99 

2215 

0.51 

0.17 

0.68 

2216 

2.90 

0.06 

2.96 

2217 

0.26 

4.23 

4.49 

2218 

2.90 

0.02 

2.92 

2219 

0.22 

6.05 

6.27 

2220 

2.90 

0.02 

2.91 

2221 

0.23 

2.16 

2.39 

2222 

2.90 

0.03 

2.93 

2223 

0.35 

0.13 

0.49 

2224 

2.90 

0.07 

2.97 


TIME= -94.77 NODE IR.POWER UV_POWER TOT_POWER < 
Orbital Averages > 


101 

48.09 

62.76 

110.85 

102 

48.13 

3.00 

51.12 

103 

409.41 

50.35 

459.76 

104 

47.15 

62.36 

109.51 

105 

47.17 

515.25 

562.42 

106 

6.47 

67.75 

74.22 

301 

14.64 

22.69 

37.33 

302 

15.44 

25.09 

40.53 

303 

14.64 

22.56 

37.19 

304 

22.14 

10.40 

32.54 

501 

14.53 

22.71 

37.25 

502 

15.39 

25.12 

40.51 

503 

14.61 

22.55 

37.17 

504 

22.10 

10.42 

32.52 
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701 

14.54 

182.15 

196.69 

702 

22.06 

79.68 

101.74 

703 

14.62 

182.09 

196.71 

704 

22.09 

79.68 

101.77 

901 

14.65 

1.15 

15.80 

902 

22.13 

0.77 

22.91 

903 

14.65 

1.15 

15.80 

904 

22.13 

0.77 

22.91 

1101 

0.00 

30.04 

30.04 

1102 

0.00 

30.04 

30.04 

1103 

0.00 

30.04 

30.04 

1104 

0.00 

30.04 

30.04 

1301 

3.51 

19.98 

23.49 

1302 

3.54 

19.97 

23.51 

1303 

3.54 

14.07 

17.62 

1304 

3.54 

14.07 

17.62 

1401 

95.30 

36.92 

132.22 

1402 

0.00 

82.37 

82.37 

1403 

94.07 

37.09 

131.16 

1404 

0.00 

82.37 

82.37 

1405 

92.58 

100.72 

193.30 

1406 

0.00 

82.37 

82.37 

1407 

93.78 

5.33 

99.11 

1408 

0.00 

82.37 

82.37 

2201 

0.31 

0.35 

0.66 

2202 

2.51 

0.72 

3.23 

2203 

0.25 

0.99 

1.24 

2204 

2.51 

0.31 

2.82 

2205 

0.10 

2.40 

2.50 

2206 

2.51 

0.24 

2.75 

2207 

0.05 

2.24 

2.29 

2208 

2.51 

0.24 

2.75 

2209 

0.06 

0.75 

0.82 

2210 

2.51 

0.31 

2.82 

2211 

0.16 

0.21 

0.37 

2212 

2.51 

0.72 

3.23 

2213 

0.64 

0.50 

1.15 

2214 

2.90 

2.53 

5.43 

2215 

0.51 

1.34 

1.85 

2216 

2.90 

0.66 

3.56 
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2217 

0.26 

5.31 

5.57 

2218 

2.90 

0.30 

3.19 

2219 

0.22 

5.04 

5.26 

2220 

2.90 

0.30 

3.19 

2221 

0.23 

0.93 

1.16 

2222 

2.90 

0.66 

3.56 

2223 

0.35 

0.21 

0.57 

2224 

2.90 

2.53 

5.43 


ITAS ABSORB SUB PROGRAM FINISHED NORMALLY 
*************************************** 

ABSORBED HEAT CALC CPU TIME (second) = 605.660 

*************************************** 


ITAS THERMAL ANALYZER: 

vvvvvvvvvvvvvvvvvvvv ADDING USER NODES vvvvvvvvvvv 


IPTR_N= 

2522 





RECORD 

2522=199 

-5420 0 

2299 

-298 0 

RECORD 

2523=20301 

-596 

5 

20302 

-596 0 

RECORD 

2524=20303 

-596 

5 

20304 

-298 54 

RECORD 

2525=20501 

-596 

5 

20502 

-596 0 

RECORD 

2526=20503 

-596 

5 

20504 

-298 54 

RECORD 

2527=20701 

-596 

10 

20702 

-596 35 

RECORD 

2528=20703 

-596 

10 

20704 

-596 35 

RECORD 

2529=20901 

-596 

10 

20902 

-596 35 

RECORD 

2530=20903 

-596 

10 

20904 

-596 35 

RECORD 

2531=1199 

-298 

0 

1299 

-894 20 

RECORD 

2532=1399 

-298 

0 




/wwwwwAvwwvwg^j) OFUSER NODES ^wwwwsa 
A fWVVWV ^^^ 


/WWVVw ^^ 


NODE 

LIST 

101 (REL NODE 

1 ) IS 

NODE 

LIST 

102 (REL NODE 

2) IS 

NODE 

LIST 

103 (REL NODE 

3) IS 

NODE 

LIST 

104 (REL NODE 

4) IS 

NODE 

LIST 

105 (REL NODE 

5) IS 


BEING ADDED TO THE CURRENT 
BEING ADDED TO THE CURRENT 
BEING ADDED TO THE CURRENT 
BEING ADDED TO THE CURRENT 
BEING ADDED TO THE CURRENT 
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APPENDIX J - STRUCTURAL DESIGN CALCULATIONS 


A. CENTRAL CYLINDER 

Assuming the payload mass to be concentrated at its center and the spacecraft 
bus mass to be concentrated at its center as shown in Figure J.l, the axial load, P, 
and the bending moment, M, at the bottom of the cylinder are determined as 
follows: 

P = (363+900)(7.2)(9.81)(1.5) = 1.34X10 5 N 
M = 363(2. 5)(9.81)(1.5)(0.762) = 1.02xl0 5 Nm 

The critical axial buckling load and critical bending moment are: 

P„=1.2-|Et 2 

M cr = 0.6yiErt 2 
where y= 1 - 0.9(1 - e -4 * 1 ) 

Using the radius of the cylinder, r = 0.254 m, Young’s Modulus for 6061-T6 

Aluminum, E = 67 x 10 9 N/m 2 , and assuming a thickness, t = 2 mm, the critical 

buckling load and bending moment are 

P = 1.75xl0 5 N 

CT 

M„ r = 6.99x1 0 4 Nm 

cr 

The margin of safety is thus: 

MS = (P/P CT + M/M cr ) 1 - 1 = 0.100 

B. ADAPTER CONE 

The axial load, P, and the bending moment, M, at the bottom of the adapter 

cone are determined as follows: 

P = (363+900)(7.2)(9.8 1 )( 1 .5) = 1.34X10 5 N 
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M = 363(2.5)9.81(1.5)1.27 + 900(2.5)9.81(1.5)0.508 = 3.38xl0 1 * * 4 Nm 

The critical buckling load and bending moment are given by 

P CT = 2 n r x t o c cos a 

M ct = 7i rj 2 t ^ cos a 

where 

a =0.19972 E-J- 

C r e 

a,, = 0.24814 E -f- 

r e 

r e = r 2 / cos a 

Using a = 41.85 ° (Figure J.l) , the smaller radius, ri = 0.254 m, the larger radius, 
T 2 = 0.709 m, and assuming the shell thickness, t = 5 mm,the critical buckling load 
and bending moment are 

P cr = 4.18xl0 5 N 
M cr = 6.59x1 0 4 Nm 

The margin of safety is therefore: MS = 0.201. 

C. HONEYCOMB PANELS 


1. Payload Panel 

This panel supports the load of the payload (363 kg). It transmits this load 

to the central cylinder by being clamped to the top of the cylinder. To determine 
the dimensions of the panel, it was assumed to be made up of two panels, each 
clamped on one side and simply supported on the other three as shown in Figure 
J.2. Each of these panels supports half the payload or approximately 185 kg. The 
mass per unit area, y, is 


y = 185 

1 (0.762)(1.524) 


159.3 — 
m 2 
m 
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Assuming a natural frequency of 50 Hz to avoid coupling with the launch vehicle, 
the panel stiffness, D, can be found from 


f= ^r Pi V 

50 = ^- 17.33 



V 


D 

159.3(0.762) 4 


D = 17650 Nm 

Assuming that the skin thickness, t, of the honeycomb is 1.0 mm, the heighth of the 


honeycomb can 


be found from the following equation: 


D = 


E t h 


2 ( 1 - 0 ) 

67xl0 9 (0.001) h 2 
2(.89) 


h = 2.34 cm 


Assuming a uniform dynamic acceleration of 36 g's across the panel, the maximum 
stress in the face skin at the center of the panel is given by: 


P 2 yg (dyn accel) a 

°max= 67h 

0.6102(159.3)(9.81)(36)(.762 2 ) 

6(.001)(0.0234) 

= 141.9 N/mm 2 

This maximum stress is below the allowable yield stress of aluminum which is 240 
N/mm 2 . Hence, the panel design is acceptable. 


2. Equipment Panels 

There are two equipment panels which must support the loads from the 
various subsystems which total approximately 600 kg. It was assumed that 300 kg 
must be supported by each equipment panel. The top panel was modeled 
identically to the payload panel. The mass per unit area is 129.2 kg/m 2 . Following 
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a procedure similar to that above, the maximum stress is 115.1 N/mm 2 . This is 
again below the allowable yields stress. 

The bottom equipment panel was assumed to be made of four panels, each 
clamped on one end and simply supported on the other three as shown in Figure 
J.3. Each of these panels supports 1/4 of the 300 kg. Using the same honeycomb 
as above, the maximum stress was determined to be 54.2 N/mm 2 , which is again 
below the yield stress considerably. 

3. Side Panels 

Each of the side panels are manufactured from the same honeycomb as 
the other panels. The panels are simply supported on all sides and are attached to 
the frame's square tubing. The panels adequately support any components 
attached to them. 
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pi = 17.33 
p2 = 0.6102 


Figure J.2 Side Pannel Diagram 


^-0.762-^ 

S 

pi =23.65 
p2 = 0.2874 

S 

Figure J.3 Top Pannel Diagram 


♦ I 

0.762 S 

I L 
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APPENDIX K - COST ANALYSIS 

The following are a few of the definintions used in cost analysis 

(A) Component: A component is a functional unit that is viewed as an entity 
for purposes of analysis, manufacturing, maintenance, or record keeping. 
Examples are hydraulic actuators, valves, batteries, electrical harnesses, and 
individual electronics packages such as transmitters, receivers, or multiplexers. 

(B) Subsystem: A subsystem is an assembly of two or more components 
including the supporting structure to which they are mounted, and any 
interconnecting cables or tubing. A subsystem is composed of functionally related 
components that perform one or more prescribed functions. Typical space vehicle 
subsystems are electrical power, attitude control, thermal control, telemetry, 
structure, and propulsion. 

(C) Space Vehicle: A space vehicle is a complete, integrated set of 

subsystems and components capable of supporting an operational role in space. A 
space vehicle may be an orbiting vehicle, a major portion of an orbiting vehicle, or 
a payload which performs its mission while attached to a recoverable launch 
vehicle. 
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